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^ ' Abstract 

o ; 

We review the present status in the theoretical and phenomenological understanding of char- 
monium and bottomonium production in heavy-ion coUisions. We start by recapitulating the basic 
notion of "anomalous quarkonium suppression" in heavy-ion collisions and its recent amendments 
^ ' involving regeneration reactions. We then survey in some detail concepts and ingredients needed 

for a comprehensive approach to utilize heavy quarkonia as a probe of hot and dense matter. 
The theoretical discussion encompasses recent lattice QCD computations of quarkonium proper- 
ties in the Quark-Gluon Plasma, their interpretations using effective potential models, inelastic 
rate calculations and insights from analyses of electromagnetic plasmas. We illustrate the pow- 
erful techniques of thermodynamic Green functions (T-matrices) to provide a general framework 
for implementing microscopic properties of heavy quarkonia into a kinetic theory of suppression 
and regeneration reactions. The theoretical concepts are tested in applications to heavy-ion reac- 
tions at SPS, RHIC and LHC. We outline perspectives for future experiments on charmonium and 
bottomonium production in heavy-ion collisions over a large range in energy (FAIR, RHIC-II and 
LHC). These are expected to provide key insights into hadronic matter under extreme conditions 
using quarkonium observables. 
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1 Introduction 

The properties of heavy quarkonium states (charnionium and bottomonium) in a hot and dense QCEf^ 
medium have been intensely studied for over 20 years now, both experimentally and theoretically. 
Current and upcoming heavy-ion collision experiments at the Relativistic Heavy-Ion Collider (RHIC, 
at Brookhaven National Laboratory (BNL) in New York), at the Large Hadron Collider (LHC, at the 
European Organization for Nuclear Research (CERN) in Geneva) and at the Facility for Antiproton 
and Ion Research (FAIR, at the Helmholtz Center for Heavy-Ion Research (GSI) in Darmstadt) put a 
large emphasis on heavy quarkonium programs in their campaigns. The interest in heavy quarkonia in 
medium is motivated by their unique role in the diagnostics of the highly excited medium created in 
ultrarelativistic heavy- ion collisions (URHICs). Early on, J/ip suppression in URHICs was suggested 
as a signal of the formation of a Quark-Gluon Plasma (QGP) pLj. This idea was instrumental in 
triggering a corresponding experimental program at the CERN-SPS. The experimental results have been 
accompanied and pushed forward by a broad spectrum of theoretical work (see, e.g., Refs. |21 El IH E] 
for various reviews). After many years of analysis and interpretation of the SPS data [B], with a first 
round of RHIC results completed and with new insights from the theoretical side (including thermal 
lattice QCD, effective models and phenomenology), it is timely to assess the current state of affairs to 
help facilitate the next stage of developments. In the remainder of this introduction we will give an 
initial view of the physics of quarkonia in a hot and dense medium, illustrating some of the difficulties 
in the interpretation of (charmomium) observables in URHICs at SPS and RHIC. 

As a starting point, we collect in Tabs. [T] and [2] basic properties of the bound-state spectrum of a 
heavy quark {Q=c,b) and its antiquark (Q) in the vacuum (note that the lifetime of the top quark is 
too short to allow for developing a ti bound-state spectrum). These spectra can be well understood in 
terms of a potential-model approach, where the underlying potential is of the so-called Cornell-type [8] , 
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consisting of a (color-) Coulomb term dominant at small Q-Q separation, r, and a linearly rising 



^Quantum Chromodynamics, the theory of the strong force and part of the Standard Model of elementary particle 
physics. 
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Table 1: Selected properties of cc mesons (generically referred to as \E') in the vacuum, as extracted 
from the recent review of particle physics [7J . The particle name in the first column is supplemented by 
a nonrelativistic classification of principal quantum number (n) and orbital angular momentum (L), the 
second column gives the total spin, parity and charge-conjugation quantum numbers, the third column 
the meson's mass (errors below 1 MeV are not quoted), the fourth column the total decay width, and 
the last column the mass difference to the open-charm threshold (taken as 2m£)0 ^ 3730 MeV). 
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Table 2: Same as Table [H but for low-lying bb mesons (generically referred to as T) and an open-bottom 
threshold of 2mB± — 10558 MeV (mass differences have been rounded to steps in 5 MeV). 



("confining") term at large r. The potential description is now understood as a low-energy effective 
theory of Quantum Chromo dynamics (QCD) utilizing an expansion in the inverse heavy-quark mass 
(l/mg) [lOl [9]. Moreover, the pertinent vacuum potential has been computed in lattice QCD [11] 
and found to agree well with the phenomenological Cornell potential as deduced from applications to 
quarkonium spectroscopy. 

The good understanding of the vacuum properties of heavy quarkonia in a relatively simple frame- 
work is one of the reasons why they are believed to be a good probe of medium effects. The latter can 
be roughly categorized into screening effects in the two-body potential and dissociation reactions with 
constituents of the heat bath. Since the heavy-quark mass is large compared to the typical temperatures 
realized in a heavy-ion reaction, a further practical benefit emerges: heavy-quark production is believed 
to be largely restricted to the earliest phase of the collision, i.e., in primordial "hard" (high-momentum 
transfer) collisions of the incoming nucleons. On the one hand, this implies a separation of the (hard) 
production process from the subsequent (soft) medium effects, and, on the other hand, it provides a 
baseline to determine the initial abundance prior to the formation of the medium. For total charm 
production, this picture is consistent with current experimental information [121 [13] and also supported 
by theoretical estimates [2]. 

The above setup defines the basic framework to analyze modifications in the production of heavy 
quarkonia in URHICs, due to final-state interactions induced by the surrounding hot/dense medium. 
The first and most widely discussed probe in this context is the J /ip. Somewhat contrary to the expecta- 
tions at the time, which considered the production of charmonium as a primordial plasma probe P^ [TB] , 
the suggestion by Matsui and Satz [T] asserted that with increasing centrality in URHICs a suppres- 
sion of the J /if) peak in the dimuon invariant- mass spectrum should occur. As a suitable reference 
Drell-Yan dileptons at high mass were suggested, as a well-established "hard" process. The underlying 
mechanism for the dissociation of the charmonium bound-state in a dense medium was associated with 
the "Debye screening" of the binding potential, largely driven by a deconfinement of color charges and 
thus intimately connected to the formation of a QGP. This phenomenon is rather general, very simi- 
lar to the dissolution of atomic bound states in electromagnetic plasmas, or to the Mott transition in 
semiconductors pTf [18] where under high pressure electrons become delocalized and a conduction band 
emerges, signaling the plasma state. 

Shortly after the prediction of J /ip suppression, the NA38 experiment at the CERN-SPS found 
evidence for this effect in collisions of 200 AGeV oxygen (O) projectiles with uranium (U) nuclei [TH], see 
Fig. [TJ The interpretation of this observation in terms of QGP formation was immediately challenged by 
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Figure 1: First observation of the J/ip suppression effect in 0(200 AGeV)-U collisions in the NA38 
experiment at CERN-SPS. When comparing the invariant-mass spectrum of muon pairs produced in 
peripheral collisions (characterized by a small transverse energy, Et < 34 GeV; left panel) with that in 
central collisions (at high transverse energy, Et > 85 GeV; right panel), a reduction of the J ftp signal 
relative to the Drell-Yan continuum is apparent (figures from Ref. [T9]). 



more conventional explanations in terms of inelastic J/^p scattering on "primordial" target and projectile 
nucleons |20]. and/or on secondary produced hadronic "comovers" [2T|[22]. The systematic experimental 
analysis of the nuclear mass-number dependence of J/ijj production in proton-nucleus (p-A) collisions 
at CERN (NA3 [23], NA38 (IHlElj) and FNAL (E772 [25]), as well as in collision systems with light-ion 
projectiles (0-U, S-U, 0-Cu) at CERN (NA38 [26j), indeed suggests that the J/ip suppression in all these 
experiments can be understood in a unified way by primordial "nuclear absorption" of the charmonium 
state (with no further reinteractions before being detected via its dimuon decay). The magnitude 
of this suppression has been characterized by an empirical mean path length, L, of the charmonium 
traveling through cold nuclear matter at normal density, Uq = 0.16 fm~^. The suppression systematics 
could then be quantitatively described following the simple absorption law, S'nucL_= exp (— o'abs'^'O-^) 
|27j . with the nuclear absorption cross section (Tabs extracted from p-A collisionqj. In fact, nuclear 
absorption systematics could not only account for p-A but also for A-B collision systems, see Fig. |2l 
where L = La + Lb is a schematic measure for the combined path length through projectile and target 
nuclei. The QGP signal had seemingly vanished p7]! 

The situation changed when the Pb beam became available at the SPS. In central Pb(158 AGeV)-Pb 
collisions the successor experiment of NA38, NA50, found a stronger J/ip suppression than predicted 
by the extrapolation of the nuclear absorption law (based on the Glauber model) [2S1 [22], see Fig. El 
This deviation was dubbed "anomalous J/ip suppression", as opposed to the "normal suppression" 
caused by nuclear absorption (following the simple exponential dependence on the L- variable). A first 
interpretation of the anomalous suppression as a signal for QGP formation was given by Blaizot and 



^A more detailed account underlying such analysis, based on the Glauber model including realistic nuclear density 
profiles, will be given in Sec. 13.21 of the main text. 
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Figure 2: Early comparisons of nuclear-absorption calculations for J ftp suppression to data in p-A and 
A-B collisions with light-ion projectiles p7]. Recent data and analyses are discussed in the main text. 



OUitrault [30] within a "threshold-suppression" scenario. It was assumed that a charm-anticharm quark 
pair cannot hadronize into a charmonium state when it is produced in a region where the transverse 
density of participant nucleons exceeds a certain "critical" value, n^"*, see Fig. HI This scenario provided 



a consistent explanation of the available data in the sense that the value for n^'''* (characterizing the 
maximal energy density of the subsequently formed medium) required to describe the suppression in 
Pb-Pb was larger than the maximum value Up reached with light-ion projectiles (see also Ref. [3T]). 
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Figure 3: Anomalous J/ if) suppression in Pb(158 AGeV)-Pb collisions at the CERN-SPS measured by 
the NA50 experiment as a function of the Glauber variable L (left panel, from Ref. j28j) and as a 
function of the number of charged particles N^h (right panel, from Ref. [32] ) . 
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Figure 4: The transverse density of participants, np, in S-U (left) and Pb-Pb (right) colhsions along 
the direction of the impact parameter vector, b, for different values of 6 = 0,2,4, ...fm. The dashed 
line indicates the adopted threshold value rip"* = 3.3 fm~^ for the onset of anomalous J/ip suppression. 
Figure taken from Ref. j30] . 

This initial interpretation was followed by a wide variety of theoretical investigations including 
hadronic comover dissociation, parton-induced break-up reactions in the QGP, pre-equilibrium effects, 
and combinations thereof. The NA50 collaboration was led to conclude that the J/ip suppression pattern 
provides evidence for the deconfinement of quarks and gluons [35] Q One of the controversially discussed 
issues was (and still is) whether the data provide evidence for a "threshold" or "onset" behavior in terms 
of step-like patterns in the centrality dependence. Further clarification was hoped to be gained from 
studying intermediate size collision systems, as recently done by the NA60 collaboration [331 El]- In 
In(158 AGeV)-In collisions, the "onset" of anomalous J/ip suppression does not seem to follow a scaling 
with the Glauber variable L, but rather with the number of participants, Npart, in the collision. Note 
that Npart, contrary to L, is a quantity closely related to the density of secondary particles, potentially 
forming a QGP, see Fig. [51 

In the year 2000 the experimental program at RHIC commenced. First J/ip data for Au-Au colli- 
sions became available in 2003 [36], while recent ones may be found in Ref. ^7]. Figure El shows the 
experimental data for J/ip production in terms of the so-called nuclear modification factor, Raa, as a 
function of centrality. It is defined as the yield observed in heavy-ion reactions relative to the one in 
p-p collisions scaled by the number, iVcou, of binary N-N collisions. 
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In the absence of medium effects one expects Raa=^- An anomalous suppression of charmonium pro- 
duction was confirmed in semi-/central Au-Au collisions, with a magnitude similar to the one observed 
at the SPS. This seems rather surprising in view of the factor of 10 higher collision energy at RHIC, 
inducing higher (initial) energy densities. However, it had been predicted [SI HSl SHI HI] that at RHIC 
and LHC the copious production of open charm leads to an additional source of J/ip production, through 
the recombination of charm quarks (or charmed mesons) in the hot and dense medium. Calculations 



^For the final NA50 results concerning J/ip suppression, see Ref. 
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Figure 5: "Anomalous" J/ip production (i.e., relative to expectations from normal nuclear absorption) 
in In(158 AGeV)-In collisions measured by the NA60 experiment compared to other projectile-target 
combinations at the CERN-SPS. The "threshold" for anomalous suppression appears to scale with the 
number of participants (right panel) rather than with the Glauber variable L (left panel); figures taken 
from Ref . 



including both suppression and regeneration mechanisms 01] had anticipated that a stronger "anoma- 
lous" suppression is largely compensated by secondary production via charm-anticharm coalescence. 
This is illustrated in the left panel of Fig. |6] where a comparison with reaction kinetic models for J/ip 
production [3S1, fSU] shows that a naive extrapolation of "anomalous" J/ip suppression in a QGP environ- 
ment from SPS to RHIC fails. Additional support for a secondary component of the J/ip yield at RHIC 
follows from the observation that J/ip suppression is less pronounced at central rapidity relative to 
forward rapidity, since at central rapidity the regeneration effect should be larger due to a larger charm 
content in the medium. On the contrary, dissociation mechanisms should increase the suppression due 
to larger energy densities at midrapidity. 

The above discussion has given a first indication of the complexity in describing and understanding 
the various facets figuring into quarkonium observables in heavy-ion collisions, let alone their interpre- 
tation in terms of properties of the Quark-Gluon Plasma. It is thus mandatory to identify key concepts 
and combine them into a comprehensive approach that allows for a systematic interpretation of ex- 
perimental data under a broad range of conditions (encompassing, e.g., all available collision energies, 
from FAIR via SPS and RHIC to LHC). This will be the main objective of this article, by means of 
a critical review of existing approaches and their comparison to data. An essential ingredient to such 
approaches are controlled reference points, e.g., equilibrium limits for in-medium properties, or p-A 
collisions to separate "normal" nuclear effects from hot/dense matter effects. Our article is roughly 
organized according to the following decomposition of the problem: 

(1) Equilibrium properties of quarkonia (Sec. [2]): these provide the main link between properties 
of the QGP and the medium created in heavy-ion collisions. The basic quantity to be determined is the 
thermodynamic quarkonium spectral function: its pole mass (binding energy) largely determines the 
equilibrium abundance (in connections with the properties of open charm/bottom), while in- /elastic 
reaction rates determine thermal and chemical relaxation times and are encoded in the spectral widths. 

(2) Quarkonium transport (Sec. 13.11) . which is required to evaluate the coupling of the quarkonia to 
the medium by computing the evolution of their phase space distribution, starting from realistic initial 
conditions. In principle, this evolution progresses toward equilibration, e.g., as a function of trans- 
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Figure 6: J/ip production, relative to the binary-collision scaled yield in p-p collisions, measured in 
the PHENIX experiment in y/sj^ = 200 GeV Au-Au collisions at RHIC. Left panel: comparison with 
kinetic models [381 ES] of charmonium production at central rapidity, illustrating the importance of 
secondary J/ip production via charm-anticharm recombination (figure from Ref. [IQ]). Right panel: in- 
crease of J/ijj suppression at forward relative to central rapidity, suggestive for regeneration mechanisms 
(from Ref. [37]). 



verse momentum. The method of choice is a transport treatment incorporating regeneration reactions 
to enable a relaxation toward thermal and chemical equilibrium. Since regeneration processes involve 
open-charm/-bottom spectra, a reliable assessment of the latter becomes mandatory. 

(3) Initial spectra and pre-equilibrium interactions (Sec. 13.2]) : these determine the initial condi- 
tion of the transport equation in item (i); they include modifications relative to spectra in p-p collisions 
{Pt broadening via Cronin effect, shadowing), as well as absorption and rescattering on primordial 
nucleons and secondary particles in the pre-equilibrium evolution of the medium. 

(4) Observables (Sec. 13. 3p : these need to be computed by combining all of the above. An additional 
input is a realistic modeling of the bulk medium evolution, ideally in terms of a locally thermalized 
medium (if applicable), e.g., within a hydrodynamic evolution, or in terms of transport models in 
regimes of incomplete thermalization. 

In Sec. m we give an outlook on the capabilities of, and physics questions to be addressed with, future 
experiments at LHC, RHIC-II and FAIR, while Sec. contains a brief summary. 



2 Equilibrium Properties 



As emphasized in the Introduction, the basic quantity encoding the equilibrium properties of a quarko- 
nium state of given quantum number a, is its spectral function, p^. Roughly speaking, Pa{uj) character- 
izes the spectral distribution of strength as a function of energy, u. Bound or resonance states manifest 
themselves as peaks with well defined mass, rria, and spectral width, F^ (in principle, more than one 
state is possible for a given a, e.g., in vacuum J/ip, ip' etc. for J^*^ = 1 , recall Tab. [1]). Note, 
however, that in the context of charmonium production in heavy-ion collisions, it is unlikely that spec- 
tral modifications can be directly detected in the dilepton invariant- mass spectrum of (vector-meson) 



quarkonium decaycl. Thus, relevant connections of the in-medium spectral properties of quarkonia to 
URHIC phenomenology are: (a) their masses, rriq,, which determine the equilibrium abundances, (b) 
their inelastic widths, V^"^^, which determine formation and destruction rates (or chemical equilibration 
times), (c) their dissolution temperatures, T^^^^, which determine the absence of formation processes for 
T > T^^^^, and, (d) their elastic widths, F^, which affect momentum spectra (and determine the kinetic 
equilibration times). 

First-principle information on quarkonium properties can be obtained from numerical lattice-discre- 
tized computations in QCD at finite temperature. We start our discussion from this perspective by 
reviewing the current status of this approach in Sec. l2.1[ As we will see, lattice-QCD (IQCD) simulations 
do not directly provide the physical spectral function, nor do they readily allow for insights into the 
mechanisms underlying "observable" medium effects. Effective models thus play an important role in 
the interpretation of IQCD results, as well as in furnishing quantitative input for applications in heavy- 
ion reactions. A natural starting point is the extension of the potential model, which works well for 
quarkonium spectroscopy in vacuum, to finite temperatures; this is discussed in Sec. 12. 2[ with emphasis 
on recent developments incorporating information and constraints from IQCD. Potential models are 
typically used to assess how medium effects in the 2-body interaction affect the bound-state spectrum 
(i.e., the location of the peaks in the spectral function). The thermodynamic T-matrix approach enables 
to assess both bound and continuum states. Of equal importance for phenomenological applications is 
the determination of inelastic reaction rates, both in the QGP and hadron gas (HG), which is elaborated 
in Sees. 12.31 and 12.41 respectively. The problem of in-medium quarkonia shares several features with 
the well-studied fermionic two-body problem in electromagnetic plasmas. Techniques developed for the 
latter may thus provide useful insights which we attempt to exhibit in Sec. 12.51 Regarding notation, 
we adopt 3 variants of the single heavy-quark (HQ) mass with the following meaning: rriQ: bare mass, 
uiq: constant effective mass, mq: temperature-dependent effective mass. 

2.1 Correlation and Spectral Functions from Lattice QCD 

The properties of the QQ interaction in a given hadronic (color- singlet) channel with quantum numbers 
a are encoded in the two-point correlation function of the pertinent current, ja- In coordinate space 
the correlation function (or correlator) is defined by the amplitude of creating the current at the origin, 
propagating it to point {T,r) and absorbing it, 

G„(r,f) = ((j„(r,f)jt(0,0))). (3) 

In a heat bath of temperature T the amplitude corresponds to a thermal average over the partition 
function of the system. For meson correlators, the QCD currents are given by quark bilinears, ja = 
QTaQ, where Fq, specifies the spin-flavor channel. In momentum space the imaginary part of the 
(retarded) correlation function, G^{u,'p), is commonly referred to as the spectral function, 

p^{oo,p) = -2lmG^{u,p) , (4) 

where u and p denote the energy and 3-momentum modulus, respectively (in the literature, the spectral 
function is also denoted by a^ = Pa/'^T^', to avoid redundancy in notation with the string tension and 
cross sections, we will not use this notation in the present article). In the timelike regime, uP' —]? > 0, 
Pa characterizes the physical excitation spectrum in the channel a, which, in principle, can be measured 
in experiment. 



''E.g., the J/tp lifetime in free space is rj/^ = l/F °, '^'^ ~ 2000 fm/c, compared to a typical fireball lifetime of ~10 fm/c. 
This means that the ratio of J/tp decays inside to outside the medium is about 1/200. If the average in-medium J/ip width 
is ^200 MeV, a detection of this effect would require a mass resolution of about IMeV, corresponding to an unrealistic 
0.03% at the J/tp mass. 



with the finite-T kernel 



In field theory, the implementation of temperature into the partition function involves the transfor- 
mation of the vacuum transition rate in real time to imaginary time. This amounts to replacing time 
evolutions with thermal averages. The temporal direction is thereby restricted to the interval [0, /3] with 
the upper limit being identified with the inverse temperature of the system, /3 = 1/T. The evaluation of 
thermal expectation values of correlation functions, Eq. (j3]), is routinely performed in IQCD. At large r, 
an exponential decay of the spatial correlator renders the lowest mass state dominant, Ga{r) oc e"'""'"^'', 
which can be used to extract hadronic screening masses, m^^, at finite temperature [l5l |16]. For tempo- 
ral correlators, usually projected onto a fixed 3- momentum, G{t,p), the limited extent of the r interval 
renders an extraction of physical ground-state masses much more difficult. For a better comparison of 
quenched (A^/=0) and full (A^/=2,3) lattice QCD, it is customary to discuss temperature dependencies 
in units of Tc in the respective simulation. This accounts at least qualitatively for the difference in the 
absolute values of critical temperature by comparing the results at roughly equal parton densities (which 
is often the relevant quantity in the discussion of medium effects). For example, for j'qi'^nch^ 1.5T^^^^, 
the factor of ~3 larger gluon density in quenched QCD is roughly compensated by the extra quarks in 
Nf=3 QCD. In addition, the quantitative values for the (pseudo-) critical temperatures computed in 
unquenched IQCD, for a given number of quark fiavors, are currently beset with a systematic uncer- 
tainty of ca. ±10-15% corresponding to T^ ~ 180±30MeV |228j (the uncertainty is smaller in quenched 
QCD with Tc~ 270 MeV). 

As mentioned above, the information on the excitation spectrum in a given hadronic channel is 
encoded in the spectral function. The temporal correlator is related to the spectral function via 

G^iT,p;T) = I ^ p^icu,p;T) K{lo,t;T) (5) 



^ ' ' ^ sinh[u;/2T] ^ ' 

The extraction of the spectral function thus requires an inverse integral transform, which, on a finite 
number of lattice points, is, in principle, not possible in a unique way. In Ref. |17] it has been suggested 
to employ the so-called Maximum Entropy Method (MEM) to perform a probabilistic reconstruction of 
the most likely spectral function. Subsequently, this method has been widely used to extract hadronic 
spectral functions from IQCD correlators ^421 SH EHl EB IS2] • A few examples in various charmonium 
channels are collected in Fig. [71 The general picture emerging from these calculations is that the 
ground state S-wave charmonia, i.e., the rjc (J^=0~) and J/ip (J^=l^), are surviving in the QGP 
up to temperatures of ~2Tc, while the P-wave states, Xc (-'^'^=0^,1^,2+), dissolve at temperatures 
slightly above the critical one, ~1.2Tc or so. The first excited S-wave state {ip'), cannot be resolved 
numerically, but presumably dissolves earlier than the stronger bound Xc states, possibly even below Tc. 
The ip' should therefore be quite sensitive to modifications in the hadronic phase of a heavy-ion collision. 
Ground-state bottomonia seem to survive to even higher temperatures, 2.5-3 Tc or more [511 [53] . 

It is instructive to examine the underlying temporal correlators, examples of which are displayed in 
Fig. M (at vanishing 3-momentum, p=0). To facilitate their interpretation, they have been normalized 
to a so-called reconstructed correlator, which is evaluated with the kernel at temperature T, 

G^ecir; T)=l£ Paico; T*) K{uj, r; T) , (7) 



but employs a spectral function at a low temperature, T*, where no significant medium effects are 
expected. The correlator ratio, 

R^a{r,T) = Gc.{r-T)/Gl^{r-T), (8) 
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Figure 7: Charmonium spectral functions as a function of energy (u) and at vanishing 3-monientuni 
(p=0) as evaluated from IQCD correlation functions using the maximum entropy method (MEM). 
Upper panels: rjc (left) and Xc (right) in quenched QCD (with Tc=295MeV), compared to their vacuum 
spectral functions (dotted lines) [51]; the 5- wave (rjc) signal at 1.5 Tc is essentially unmodified compared 
to the vacuum, while the P-wave (xc) signal at 1.16 Tc is largely distorted (the "data" points in the 
upper left panel, as well as the two solid lines in the right panel, indicate the uncertainty in the 
MEM reconstruction procedure). Lower panels: J/ip (left) and Xc (right) computed in full QCD for 
2 light flavors [52] (the y scale is in arbitrary units); the critical temperature in these simulations is 
Tc ~ 210 MeV (slightly larger than typical values in the range of 160-190 MeV); a J/ip signal around 
a;~3 GeV survives up to temperatures of at least 392 MeV~1.9 Tc while the Xc peak is essentially smeared 
out at T=252MeV~1.2Tc (the quark masses, m, quoted in the lower panels are unrenormalized input 
values; after renormalization, they give rise to physical zero-temperature charmonium masses within 
10-15 % of the experimental values). 



is then an indicator of medium effects in (^^(r; T) through deviations from 1 (note that a normalization 
to Ga{T]T*) is not meaningful due to different r ranges for T* and T). The stability of the temporal 
correlator ratios, i.e., their small deviations from one, in the vector and pseudoscalar channels has been 
associated with the survival of the ground-state S-wave charmonia. The variations for the Xc (which are 
P waves) are much larger and set in much closer to Tc (e.g., up to 40-50% at T=1.2 Tc). The opposite 
trend of the J/ip correlators in the upper left and lower left panel could be due to the different vector 
components considered (Z]/^=o-3^y' i^ the upper panel vs. I]j=i_3 Gy in the lower panel). 

Another quantity that can be computed with good accuracy in IQCD, and which encodes information 
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Figure 8: Temporal charmonium correlators at zero 3-monientum computed for VC=J /ip^ PS=t]c 
(left panels) and SC=Xco, AX=Xci (right panels) in quenched (upper panels) |19] and Nf=2 (lower 
panels) [52] lattice QCD. The critical temperature in these simulations is T^ ~ 270 MeV and 210 MeV, 
respectively. Note the comparatively small variations in the S*- waves (left panels), and the rather 
large increase in the P-waves (right panels) indicative for additional low-energy strength relative to the 
reconstructed correlator. For clarity, the PS and AX correlators in the upper panels have been offset 
by the indicated constant. 



about charmonium properties at zero and finite T, is the free energy, FQQ{r), of a static pair of heavy 
quark and antiquark (more precisely, IQCD obtains FqqIt) as the difference between free energies for 
a thermal system with and without the QQ pair). In gluo-dynamics the large distance limit of the free 
energy, F^q = Fqq^r — )• oo), can be related to the expectation value of the Polyakov loop [53] . 

e-^QV^=|(L)r; (9) 

L characterizes a heavy-quark (HQ) source whose expectation value is zero (finite) in the (de-) confined 
phase, and thus serves as an order parameter of deconfinement. At zero temperature, the free energy 
computed in the color-singlet channel (which we denote by a superscript "(1)") closely follows the form 
expected for a phenomenological Cornell- type potential [TOl El ITT] - 

^^^\r;T = 0) = —— + ar , (10) 



QQ^ 



3 r 



recall Eq. (|T]o. This finding, in connection with the development of HQ effective theories, has confirmed 



^Strictly speaking, the projection of the color-averaged free energy as computed in lattice QCD, onto color-singlet 
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Figure 9: Left panel: free energy of a color-singlet heavy quark-antiquark pair as a function of its 
size for various temperatures as computed in Nf=2 lattice QCD [llj; the vertical bar in the upper 
left part indicates the estimated string-breaking energy, V^Mrsh) — 1.1 GeV, in vacuum. Middle panel: 
temperature dependence of the asymptotic value of the color- averaged free energy, F^q(T), in quenched 
(A^/=0) [56], Nf=2 [m and A^/=3 [57] IQCD; the pertinent critical temperatures are Tc=270, 200 
and 193 MeV, respectively. Right panel: entropy contribution to the color-averaged free energy at 
asymptotically large quark-antiquark separation corresponding to the free energies of the middle panel. 
The band enclosed by the dotted lines for T > l.STc indicates a perturbative estimate of TSqq{T). 



a posteriori the suitability of potential models to quantitatively compute charmonium and bottomonium 
spectroscopy in free space based on first principles |9]. In the presence of light quarks, string breaking 
occurs, being characterized by FJj^{r;T = 0) leveling off at a finite value. This may be interpreted as 



QQ^ 



an effective quark (or D-meson) mass according to nic = m^ + FqI^°° /2 (with a bare quark mass of 
m° ~ 1.2-1.3 GeV). Numerically, with a typical string-breaking separation of r^h — 1.2 fm, one finds 
Fii(rsb;0) ~ 1.1 GeV, which roughly recovers the empirical D-meson mass. At finite temperature, 

IQCD computations for FqUt^T) find the expected color-Debye screening, which gradually penetrates 
to smaller distances as the temperature increases, cf. left panel of Fig. |9l The implications for HQ bound 
states (or spectral functions) are, however, quite subtle. First, at finite T the free energy receives an 
extra contribution due to an entropy term. 



FQQ(r;T) = f/QQ(r;T)-T5QQ(r;T) 



(11) 



It is currently not clear whether the internal energy, UQQ(r]T), or the free one, FqqIv^T), is a more 
appropriate quantity to be identified with a finite-temperature HQ potential. Second, the temperature 
introduces an extra scale into the problem which needs to be properly implemented into the construction 
of a HQ effective theory [5HI EH] . Third, the presence of the entropy term, which does not vanish at large 
separation r, renders the identification of an effective HQ mass more problematic, especially close to 
Tc, where TS^q becomes very large, see right panel of Fig. [91 Nevertheless, the application of potential 
models at finite temperature has provided valuable insights into heavy quarkonium properties at finite 
T. In particular, the synthesis of independent information on correlators and free energies from IQCD 
has enabled interesting insights, as we discuss in the following Section. 



and -octet channels is not gauge invariant i55, . However, they can be constrained at sufficiently short distances where 
temperature effects are not operative, or guided by Casimir scaling relations inferred from perturbative QCD. 
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Figure 10: Nonperturbative extraction of the Debye mass {mD=fJ'D) based on exponential fits of a 
screened Coulombic term to the in-medium color-singlet HQ free energy in Nf=2 (solid circles) [11] 
and quenched {Nf=0, open squares) [61] IQCD. The bands are analytic fits using the perturbative 
form of the Debye mass, fio/T = A{1 + A^//6)^/^5'2-ioop with a 2-loop finite-T running coupling and a 
multiplicative constant 74~1.5 to account for nonperturbative effects. Absolute units are obtained by 
using the pertinent critical temperatures of Tc ~ 200(270) MeV for iVj-=2(0). 



2.2 Potential Models in the QGP 

Early applications [60j of in-medium heavy-quark potentials have employed a phenomenological ansatz 
to implement color-screening effects in a deconfined plasma into the Cornell potential, Eq. ([T]), via 



Vooir^T) 



a 



f^D{T) 



(1 



-MD(r)r 



4:as 
3r 



-Mo(T)r 



(12) 



The entire temperature dependence is encoded in the Debye mass, fioiT) (the vacuum potential is 
recovered for fi£, — )■ 0). Defining the Hamiltonian as 



H{r;T) = 2mQ + VQQ{r;T) 



the Schrodinger equation for the bound state problem. 



H (pniir; fin) = ^uAi^d) 4>ni{r] fin) 



(13) 



(14) 



has been solved for the eigen-energies, Uni, as a function of the Debye mass. A state [n, I) was considered 
to be dissolved at a "critical" Debye mass, fijf^, if the "dissociation" energy. 



.diss 



a 



co^r = 2mQ H iOni 

fJ'D 



(15) 



reaches zero. Note that the term -^- corresponds to the large-separation limit of the in-medium potential, 
Eq. 012p . cf. also the discussion in the text following Eq. (II Op in the previous Section. The critical Debye 
masses for dissolution were found to be /i^*'^(n = 1) ~ 0.7(1.6) GeV for ground state (n=l) charmonium 

2) ~ 0.35(0.6) GeV for the first excited states, tjj' and Xc (T' and Xb) [60] . 



(bottomonium) , and ijLJf^{n 
Recent IQCD computations of the in-medium Debye mass based on the Coulombic term in the color- 
singlet free energy are displayed in Fig. [10] [11]. Comparing these to the values for /ijf^ extracted 
from the screened Cornell potential, Eq. flT^ . suggests dissociation temperatures of all charmonia at or 
slightly above T^. 



The next step became possible with the availabihty of quantitative IQCD results for heavy-quark 
free energies. In Ref. [62], the color-singlet free energy, FI^q was used in a Schrodinger equation, and 
J/ip dissolution occurred at about 1.1 Tc, while Xc (V'') were found to melt (well) below T^.. The recent 
IQCD spectral function results, indicative for J/ip states surviving well above T^., appear to challenge 
this conclusion (recall Fig. [7]). In Ref. [63] it has been suggested to employ the internal, rather than 
the free, energy as an interquark potential. It was argued that Uqq is the appropriate quantity in the 
limit where the typical bound state (formation or life) time is much smaller than the typical time for 
heat exchange with the thermal environment, Tbound "^ Theat- In this case, the entropy term is not 
active and should be removed from the free energy, leading to the internal energy. While, according 
to Eq. flTTl) . one has Uqq > Fqq, the actually used potentials become deeper upon subtraction of the 
asymptotic values of the internal energy (recall that the force is given by the derivative of the potential). 
As noted in connection with Fig. [9l the subtraction procedure is problematic in the immediate vicinity 
of Tc. This problem is currently the weakest link in converting IQCD results to effective heavy-quark 
potentials. From a pragmatic point of view, the use of Uqq in a Schrodinger equation [631 EH 165] 
improved the qualitative agreement of the potential approach with the IQCD spectral functions in that 
J/ip bound states are supported up to temperatures of ~2-2.5 T^, while the Xc dissolves at around 1.1- 
1.2 Tc. In Ref. [BB], the internal energies have been implemented into a T-matrix approach for Q-Q (and 
q-q) interactions, which has the important advantage over the Schrodinger equation that it provides a 
unified description of bound and scattering states (i.e., above the nominal Q-Q threshold). For S-wave 
charmonia, it was found that the lowest bound state moves into the continuum at ~2 Tc, after which it 
rapidly dissolves. However, it was also noted that the extractions of the internal energy from different 
IQCD calculations are not unique, thus adding to the uncertainty in the potential definition. For 
example, in perturbative approaches [Hll |82] the free energy has been identified with the heavy-quark 
potential. In the vicinity of a phase transition and/or in the presence of a confining term, this notion 
is less obvious. In Ref. [67] it has been pointed out that, in the static limit, the thermal expectation 
value of the interaction part of the QQ Hamiltonian is precisely the internal energy of the QQ pair and 
thus the appropriate quantity to serve as in-medium two-body potential. Another alternative has been 
suggested in Ref. [64j in terms of a hnear combination of F and U, based on the idea of subtracting 
the internal energy of the induced gluon cloud. Clearly, the question of extracting a two-body potential 
from the IQCD free energy deserves further study. 

Model comparisons to lattice QCD "data" can be made more quantitative by computing the perti- 
nent spectral functions, Pa{^), within a given approach (see, e.g., Ref. [68] for an early application in 
the light quark sector): the integral in Eq. ([5]) can then be easily carried out and compared to IQCD 
results for Euclidean correlation functions. In the following paragraphs we discuss such line of work 
within potential models. One idea is that the different binding properties resulting from assuming either 
F or f/ as potential can be distinguished by such comparisons. 

The first applications of a finite-T potential model for charmonia to IQCD correlators has been 
conducted in Ref. [69j. For the bound-state part a Schrodinger equation was solved using a screened 
Cornell potential, Eq. ( TT2l) . as well as IQCD internal energies, with fixed heavy-quark masses. The 
continuum was modeled as non-interacting with a temperature dependent threshold sq = so(T), leading 
to a spectral function 

Pc^iu; T) = 47r 5^ M,F^6{u^ - Ml) + ^u^e{u - so)Uuj, so) (16) 

(the kinematic coefficients fa depend on the meson channel a, Fi are wave-function overlaps at r=0, 
and the sum over i includes all bound states for a fixed a). While some trends of the lattice correlator 
ratios were reproduced, a comprehensive agreement could not be achieved; similar results have been 
obtained in Ref. [70] . 
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Figure 11: Diagrammatic representation of the Q-Q correlation function as obtained from resumming 
the HQ potential (upper panel) resulting in a folding with the T-matrix (lower panel) corresponding to 
Eq. ( 1T8|) . The solid dots represent operators projecting into different mesonic channels a (giving rise 
to the coefficients a/ in Eqs. (fT9|) and ( 120|) ) and the 2 connecting lines represent the non- interacting QQ 
propagator Gqq; figure taken from Ref. [71j. 



In Ref. [7T] the T-matrix approach introduced in Ref. [50] has been further developed for a sys- 
tematic analysis of heavy-quarkonium correlator ratios. Starting from the covariant Bethe-Salpeter 
equation for two-body scattering, one can apply standard reduction schemes to obtain a 3-dimensional 
Lippmann-Schwinger equation. After a partial-wave expansion, the latter becomes a 1-dimensional 
integral equation for the Q-Q T-matrix which for vanishing total 3-momentum (j9=0) takes the form 



Ti{u-q',q) = Vi{q',q) + - dk P Vi{q',k) Gqq{u- k)Ti{u; k,q) f'^'^M 

TC Jo 



(17) 



Gqq{u!, k) denotes the intermediate heavy quark-antiquark propagator, fQQ=l—fQ~fQ with f'^''^ HQ 
Fermi distributions, and q, q' and k are the (off-shell) 3-momenta of the quarks in the initial, final 
and intermediate state, respectively. The in-medium T-matrix equation can be derived from a finite- 
temperature Green's function approach [72] and constitutes a consistent many-body framework to 
compute in-medium 1- and 2-particle correlations (it has been widely applied, e.g., in nuclear many- 
body theory [73], as well as in the analysis of electromagnetic plasmas, as will be discussed in Sec. 12. 5p . 
It is particularly suited for the problem at hand since the potentials can be directly identified with the 
ones extracted from the lattice (after Fourier transformation and partial-wave expansion with angular- 
momentum quantum number /). It treats bound and scattering states on the same footing, which is 
mandatory for situations where bound states gradually dissolve into a continuum. At the same time, 
rescattering effects in the continuum, including possible resonance formation, are accounted for by a full 
resummation of the potential, without the need for matching procedures or K factors. The 2-particle 
propagator allows for the implementation of medium effects on the heavy quarks via (complex) single- 
particle self-energies encoding mass changes (real part) and finite widths (imaginary part). Charmonia 
widths are important for phenomenological applications in heavy-ion reactions and thus should be 
accounted for in correlator analyses. The T-matrix is directly related to the correlation function in 
momentum space by a double folding over external vertices, schematically given by 

G = G° + G° T G° = G^ + AG , (18) 

and illustrated in Fig. [TTl Performing a partial wave expansion and keeping the leading terms in l/mq 
the explicit form of the correlator becomes 



G'A^) 



2N, 

IT 



k'dkai{k,k)GQQ{uj-k) f'^'^M 



(19) 



7]-o 



AGi{u) = ^k'dkGQQ{u-k)f'^'^{ujk)k''dk'GQQ{uj-k')f'^'^{ujk')ai{k,k')Ti{uj,k,k'), (20) 
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Figure 12: Charmonium spectral functions (upper panels) and pertinent Euclidean-time (r) correlator 
ratios (lower panels), Eq. (|8]), within the T-matrix approach of Ref. [71] using the internal energy 
extracted from Nf=3 lattice QCD [57] as input potential. Note that the energy range in the upper 
panels is smaller than in Figs. [7] and [TU to better resolve the bound-state and threshold regime. The 
correlator ratios are normalized using the calculated vacuum spectral function. The left panels represent 
S'-wave states {J/i', rjc) while the right ones correspond to P-waves {Xco, Xci)- The charm-quark mass 
has been fixed at mc=1.7 GeV and a small width of 20MeV has been introduced in the Q and Q 
propagator for better display of the bound-state regime. 



with coefficients ao,i(fc, k')=2, 2kk' /rriQ in the S- and P-wave channels, respectively. The imaginary part 
of the correlator is just the spectral function (up to a constant factor). To leading order in l/mg, HQ 
spin symmetry implies degeneracy of different spin states within a partial wave (i.e., for S*- waves r]c and 
J/ip, as well as for the P-wave Xc states). Results of the T-matrix approach for spectral and correlation 
functions are shown in Fig. [12] for a constant effective c-quark mass {mc=1.7 GeV) and a small width 
(~20 MeV) in the quark propagators (for better resolution we here focus on the energy regime around 
the open-charm threshold, Uthr = ^rric)- For the potential the (subtracted) Nf=3 internal energy of 
Ref. [57] has been employed. The spectral functions displayed in the upper panels of Fig. [12] indicate 
that the J/ip (or rjc) dissolves around ~3Tc (~2.5Tc when using the internal energy of Ref. [H]), while 
the Xc states melt at lower temperatures, ~1.3Tc. These values should be considered as upper limits, 
since larger widths will lead to a melting at smaller temperatures. The continuum part of the spectral 
functions in both S- and P-waves exhibits a large enhancement over the non-interacting case; this 
effect cannot be neglected in quantitative applications to correlation functions. The pertinent temporal 
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Figure 13: Two effects at the single HQ level relevant to quarkonium correlators. Left panel: in- medium 
charm-quark mass following from the asymptotic value of the internal energy (upper curve) [571 E] 
and as inferred from a "fit" to the ?7c correlator [75] as shown in the upper right of Fig. [151 Right 
panel: contributions to temporal correlation functions computed in finite-temperature lattice QCD [76] . 
In addition to the standard mesonic point-to-point correlator (upper panel), the periodic boundary 
conditions in r direction allow for a quark scattering diagram (lower panel). 



correlator ratios, i?g (recall Eq.(|8l)), are plotted in the lower panels of Fig. [121 In the S'-wave they 
are qualitatively similar to the ones in IQCD, while in the P- waves they are quite different in both 
magnitude and temperature dependence, recall Fig. [H] In particular, the decreasing trend at large r 
indicates that the model spectral functions are missing low-energy strength with increasing temperature. 
In the ?7c channel, this problem is largely overcome when the in-medium charm-quark mass, as following 
from the asymptotic value of the internal energy, is included, i.e.. 



m:(r) = m^ + [/^^'°°(r)/2 



(l),oo. 



(21) 



displayed in the left panel of Fig. [131 This stabilizes the T-dependence of the r/c correlator ratio [71] . 
but it is not enough to quantitatively improve on the agreement with the IQCD results for J jil) and Xc 



An essential part of this puzzle has been identified in Ref. [76], in terms of an energy-independent 
(but T-dependent) contribution to the correlators which arises due to the periodic boundary conditions 
in the (finite) r direction at finite temperature (see right panel of Fig. [T3|) . This term may be inter- 
preted as a scattering contribution of a single heavy quark in the medium which generates a low-energy 
peak in the quarkonium spectral functions that can be related to a contribution from a heavy-quark 
susceptibility, Xa; to the lattice QCD correlators [77] . These contributions have been readily imple- 
mented into potential model analyses and are found to largely resolve the discrepancies with the IQCD 
correlator ratios [751 [THl [79] . While the pseudoscalar (?7c) correlator appears to be free of the zero- mode 
contribution, the P-wave correlators in the QGP are enhanced appreciably, in particular at large r (the 
J/?/' is mildly affected). This is due to the additional very- low energy strength in the spectral functions 
which, in quasiparticle approximation, is generated by extra terms ~ XabJb{iS). When implementing 
these into the correlation functions and evaluating the susceptibilities, Xa? for a free gas, the description 
of the IQCD P-wave correlator ratios (Fig. [8]) is much improved [751 [751 [79], see lower panels in Figs. [H 
and [151 However, the agreement of the potential approaches with IQCD correlator ratios still allows for 
a significant redundancy in the underlying spectral functions, leading to rather different conclusions on 
the "melting" temperatures of the various charmonia, as we will now discuss. 
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Figure 14: Quarkonium spectral functions and correlator ratios in a Gluon Plasma (quenched QCD) 
within a nonrelativistic Green's function approach employing a screened Cornell potential [78], sup- 
plemented with a perturbative QCD continuum including a K factor at high energies. Upper panels: 
spectral functions and correlator ratios (insets) for 5'-wave charmonium (77^, left) and bottomonium (77^, 
right) at various temperatures (the data points in the inset are IQCD results with the same color code 
as for the calculations represented by the lines); the lattice results are for quenched QCD [51], with 
the data points in the spectral functions indicating the statistical uncertainty (vertical error bar) of the 
MEM over the energy interval indicated by the horizontal bar (same as in upper left panel of Fig. [7]). 
Lower panels: Euclidean correlator ratios for P-wave charmonia (left) and bottomonia (right) at 1.5 T^; 
closed and open symbols are quenched IQCD results using isotropic [IHIED] and anisotropic [51] lattices. 



In Ref. [7S], a screened Cornell potential has been employed within a nonrelativistic Green's function 
approach to describe the low-energy part of the spectral function. The latter has been combined with 
a high-energy part which has been approximated with a perturbative continuum multiplied by a i^ 
factor (simulating radiative corrections) to match the large-r behavior of the IQCD data. Above T^ 
the S-wave charmonium and bottomonium spectral functions exhibit a moderate variation with T close 
to threshold, cf. upper panels in Fig. [TH The rjc state melts at a temperature as low as T ~ 1.3Tc, 
but the corresponding S'-wave correlator changes little with temperature, roughly consistent with the 
small variations found in the IQCD correlator ratios (cf. the insets in the upper left panels). It was 
therefore concluded that the IQCD correlators are consistent with a melting of all charmonium states 
at temperatures below 1.3 T^. Essential to this argument is a reduction of the in- medium charm- 
anticharm quark threshold to 2?7i* ~ 2.7 GeV, together with a nonperturbative threshold enhancement 
in the spectral function, as first emphasized in the T- matrix approach of Ref. [7T] . 

In Ref. [79] the in-medium HQ potential was taken as a linear combination of Fii and t/iA (as 
advocated in Ref. [M]), in connection with a constant HQ mass and a noninteracting perturbative con- 
tinuum, corresponding to a spectral function as in Eq. f lTB]) . Here, the (semi-) quantitative description 
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Figure 15: Charmonium correlator ratios, Rg = G/G^ec as in Eq. (IHl), including zero-mode contributions. 
Left panels: 5*- and P- waves (upper and lower left, respectively) in the Sclirodinger-equation approach 
of Ref. [79] where the HQ potential is based on a combination [6l] of free and internal energy extracted 
from unquenched lattice QCD. Right panels: rjc and Xa (upper and lower right, respectively) in the 
T-matrix approach of Refs. [Til EH] (employing Nf=3 internal energies from IQCD) including an in- 
medium charm-quark mass according to the left panel of Fig. [131 



of the IQCD correlator ratios (see left panels of Fig. [T5|) goes along with a J/ip melting at ~1.6 Tg. 

The right panels of Fig. [151 depict S- and P-wave charmonium correlator ratios within the T-matrix 
approach based on the spectral functions in the upper panels of Fig. [121 where the HQ potential is taken 
as the (subtracted) internal energy. The in-medium charm-quark mass has been calculated including 
the asymptotic value of the internal energy, m* = m° + f/iA'°°(T)/2 (except close to Tc where the large 



■,(l),oo 



entropy contribution, TS^k"", is problematic, cf. Fig. [TS]). As in Refs. [THl [79], the IQCD correlator 
ratios can be approximately reproduced, only that in this approach the dissolution temperature of the 
J/ip is around ~2.5Tc. 

In a slightly different line of work, potential models have recently been developed utilizing perturba- 
tive techniques, starting with a suitable identification of the static finite-T HQ potential [581 ESI El 112] ■ 
It has been found that, even at leading order in as, the potential develops an imaginary part related 
to Landau damping in t-channel gluon exchange (cf. Ref. ^83j for an analysis using classical IQCD). 
When applied within a Schrodinger equation in the bottomonium sector, bound-state solutions emerge 



which gradually melt with temperature [SI]; pertinent spectral functions exhibit a large threshold en- 
hancement, not unlike the ones displayed in Fig. [12] |7l]. In Ref. [85] it has been shown that, in weak 
coupling, quarkonium resonances survive up to temperatures parametrically given by T ~ g^rnq, but 
have melted at T ~ gruQ (the melting occurs at T ~ g^'^mg [58]). 

Let us briefly reflect on the main points of Sees. 12.1] and 12.21 The main strength (and best accuracy) 
of IQCD lies in the computation of euclidean correlators, which, however, are difficult to transform 
unambiguously into spectral functions containing the information on the physical excitation spectrum 
(such as the existence, masses and widths of the quarkonium states). Potential models, on the other 
hand, are less general in principle (e.g., based on a static approximation), but have the advantage of 
enabling direct calculations of spectral functions and providing insights into the underlying mechanisms 
(e.g., color screening and quark- mass dependencies). Moreover, the spectral functions can be readily 
used to calculate pertinent euclidean correlators. It is thus the synergy of IQCD and effective models 
which provides the best prospects to better understand the properties of heavy quarkonia in the QGP. 
More concretely, IQCD computations of Euclidean correlator ratios and associated spectral functions 
(Figs. Ul [H]) exhibit a weak temperature dependence for 5- wave states {rjc, J/ip), while P-wave (xc) 
correlators increase substantially. Potential models employing IQCD-motivated input potentials (or 
screened Cornell potentials) can essentially reproduce the IQCD results (Figs. [TH [T5|) . but the inter- 
pretation is not conclusive: the required low-energy strength in the S-wave spectral functions can be 
provided by either bound states surviving up to ~2Tc, or by a rather strong reduction in the open-charm 
threshold in connection with a nonperturbative rescattering enhancement. The increase in the P-wave 
correlators can be accounted for by a low-energy transport peak which is not directly related to bound- 
state physics (but may be sensitive to the charm-quark mass [86]). A quantitative determination of the 
in-medium HQ mass will therefore be important to make further progress. Another ingredient that has 
not received much attention in the context of correlators is the (inelastic) width of the charmonium 
spectral function, which is a very relevant quantity for heavy-ion phenomenology. This is the topic of 
the following Section. 

2.3 Quarkonium Dissociation Reactions in the QGP 

The survival of a charmonium state in the QGP does not imply that it is protected from suppression, 
since it can be dissociated by inelastic collisions with surrounding partonqj. By detailed balance the 
inverse of these reactions leads to charmonium formation, and therefore the pertinent inelastic width, 
r^ = (r*)~^, is directly related to the relaxation time, r^, of the charmonium abundance toward its 
thermal equilibrium number, N^. This is immediately borne out of the pertinent kinetic rate equation 
for the time evolution of the charmonium number A'^g,(r), which for a spatially homogeneous system 
takes the simple form 

^ = -r^(iV^-iV;'i) (22) 

(its solutions in the context of heavy- ion collisions will be discussed in Sec. 13.11 below). As usual, such 
a rate equation is applicable if a well-defined quasiparticle state \& exists and if the deviation from 
thermal equilibrium is not too large so that the relaxation time approximation is valid. In general, the 
dissociation of a composite quasi-particle in a medium can be classified into the following processes: 

(1) direct decays, \E'— J-c + c (or \E'—)-c + c + (y'); 

(2) inelastic collisions with partons, most notably 

^In the following, we refer to ^ = rjc, J/ipj4'' iXc, ■■■ as a generic charmonium state, and Y = r^f,, T, T', x^, ... as a 
bottomonium state. To streamline the discussion, we concentrate on charmonia. Most of the arguments also apply to 
bottomonia although regeneration effects are expected to be less important than for charmonia. 





Figure 16: Charmoniuin dissociation reactions via parton impact; left panel: gluo-dissociation [881 189] 
representing the leading-order QCD process; right panel: quasifree dissociation [H] representing a 
next-to-leading order process. 



(2i) -^ + g ^ c + c (cf. left panel in Fig. [16]) 

(2ii) \E' + 2— >c + c + i, with i = g,q,q (cf. right panel in Fig. [16 



Note that the distinction between a bound and a resonance state from the point of view of the rate 
equation ( [22]) is immaterial; it matters, however, for quantitative assessments of the dominant contribu- 
tions to the inelastic width. E.g., if a \l/ state moves into the continuum, i.e., above the cc threshold, at 
a certain temperature, the opening of the direct-decay channel usually renders it the dominant process. 
In this sense the dissolution of a state due to screening corresponds to the F^ — )■ oo limit of direct 
decays. The potential-model calculations discussed in the previous section indeed suggest that once 
a charmonium state becomes unbound, resonance states can no longer be supported. On the other 
hand, for binding energies of the order of 100 MeV, inelastic widths of a similar magnitude could induce 
a "premature" melting of the state (the dissociation temperatures quoted above are therefore to be 
considered as upper limits). 

A standard approach to compute inelastic reaction rates in a thermal medium is to calculate the 
cross section for a given dissociation reaction and fold it over the thermal distribution of the (on-shell) 
medium particles. 
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(frei: relative velocity of the incoming particles in the center of mass of the collision; s={p+k)'^: center-of- 
mass energy squared of the parton- \l/ collision). The cross section for gluo-dissociation, g + J/ip — )■ c-(-c, 
has been evaluated in Refs. [SHI [SSI 1148] , using hydrogen-like wave functions. The latter are expected 
to be a good approximation as long as the Coulombic part of the Cornell potential dominates the 
Q-Q bound-state properties. As noted in Refs. [SSI [SH] this approximation is not really (marginally) 
satisfied for the J/ip (T) with vacuum binding energies of ~0. 6(1.0) GeV. Within these limitations, the 
gluo-dissociation cross section for massless gluons of energy ko incident on a ground-state quarkonium 
reads 
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where Eb denotes the quarkonium binding energy and tuq the effective heavy-quark mass (as used, e.g., 
in fits to vacuum quarkonium spectra). The second line has been obtained by using the Coulomb- 



potential expressions for binding energy, Eb 



7 a.. 



rnq, and Bohr radius, Oq = 4/(3as'^Q)- It 



illustrates the leading-order (LO) character {0{as)) of the cross section, being proportional to the 
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Figure 17: Total cross sections (left panel) and pertinent dissociation rates (right panel) for J /ip disso- 
ciation reactions in the QGP: gluo-dissociation, g + J /tp -^ c + c, according to Eq. (!24|) with vacuum 
(eB=640MeV, dashed line) and in-medium (e^— 30MeV, dotted line) binding energies (the pertinent 
phase distribution, p'^f^{p; T), for massless gluons is illustrated by the dash-dotted line); parton- induced 
quasifree dissociation, i + J/ijj -^ c + c + i, using an in-medium binding energy, /i^i = gT and thermal 
parton masses m, oc gT 



"geometric size" of the bound state. The cross section has a threshold at k^^ = eb and acquires its 
maximum at /c™^^ = y^s, cf. left panel of Fig. [T71 For QGP temperatures of 200-300 MeV, typical 
for SPS and RHIC energies, the average thermal gluon energy of 3T amounts to 600-900 MeV, provid- 
ing good overlap with the vacuum gluo-dissociation cross section and resulting in dissociation rates of 
Tj/^ ~ 30 — 200 MeV, see right panel of Fig. [T71 However, if the binding energy decreases with tempera- 
ture, the available phase space of gluo-dissociation shrinks substantially. E.g., in the screened potential 
model of Ref. [60], corresponding to Eq. ( IT2l) . a moderate screening mass of fiD — gT with as ^0.25 
((7~1.75), yields a J/tp binding energy of only about 200(30) MeV at T=180(300) MeV, reaching zero at 
around 360 MeV. Under these conditions the gluo-dissociation rate decreases with increasing tempera- 
ture for T > 250 MeV, cf. dotted line in the right panel of Fig. [T71 This signals that other dissociation 
processes become important (see also Refs. |M1 EOl E]). For the same reason, the gluo-dissociation 
process cannot provide a realistic description for the dissociation of excited charmonia whose binding 
energies are small already in the vacuum. 

An early nonperturbative approach to compute heavy quarkonium break-up from quark impact 
(corresponding to the category (2ii) of the above classification) has been suggested in Ref. [93] in 
analogy to electron-induced atom ionization; the pertinent "thermal activation" cross section. 
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is characterized by a temperature dependent dissociation energy, k^ = Eb, while the coefficient has 
been estimated by the "geometric" transverse size of the bound state, Trr^^ (cf. Ref. [91] for an alterna- 
tive derivation of a similar result, which was shown to apply to sufficiently large dissociation energies, 
^diss ^ 2"^_ To obtain more quantitative results, microscopic approaches are required. 

In Refs. [H| 112] the dissociation reactions of weakly bound charmonia have been treated within a 
quasifree approximation for parton impact on either one of the heavy quarks within the bound state, 
2 + \E' — 7- c + c + i. The basic process is (quasi-) elastic i + c — )■ z + c scattering [98] , including the proper 
dissociation kinematics (i.e., finite Eb and 4-momentum conservation), as well as thermal parton and 
Debye masses (the latter are essential to render the dominant t-channel gluon-exchange contributions 
finite). The key quantities controlling the rate are a^, /i_D and Eb, while the dependence on rric is weak. 
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Figure 18: Parton-induced J/ip dissociation rates in the QGP. Left panels: NLO results using in- 
medium binding energies ranging from ~40— J-O MeV for T=190— ;'280 MeV (upper left: T dependence of 
rjM(p=0) for (a) quark- plus antiquark- and (b) gluon-induced reactions for constant thermal parton 



masses 



lower left: 



ptot _ 



+Tj/^g, as a function of the J/i/j 3-velocity, 



J/Vg+g-^J- J/V9' 



p/up, 



relative 



to the heat bath for different temperatures with constant thermal parton masses of 600 MeV [^S] and 
Tc=170MeV). Right panels: comparison of the 3-momentum dependence of J/ip (upper right) and Xc 
(lower right) dissociation rates at 3 different temperatures using either gluo-dissociation with vacuum 
binding energies [881 IM] (dashed lines) or quasifree destruction with in-medium binding energies |Hl [97] 
(solid lines). 



Formally, these reactions are of next-to-leading order (NLO) compared to LO gluo-dissociation, but 
the available phase space in the cross section is substantially increased (see solid line in the left panel 
of Fig. [TT]) . leading to significantly larger dissociation rates at temperatures where the charmonium 
binding is weak (see solid line in the right panel of Fig. [TT]) . In fact, the infrared singularity in the 
t-channel exchange of a (Debye-screened) gluon for inelastic quarkonium scattering off thermal partons 
renders this process of order g'^. 

A complete NLO calculation has recently been performed in Ref. [HS] confirming the prevalence 
of i + J/ip — )■ c + c + i over g + J/ip -^ c + c processes for weak binding. More quantitatively, for 
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Figure 19: Bottomonium dissociation rates in the QGP. Left panel: quasifree calculation [99j, T + i — >^ 
6 + 6 + z (i=g, g, g), using a^— 0.25 with vacuum (lower curves) and in-medium (upper curves) binding 
energies for T, %{, and T' (the in-medium binding energies are based on Ref. [00], with £:^~550(200) MeV 
at T=180(350) MeV). Right panels: NLO calculations for quark- and gluon-induced break-up of the 
ground-state T [95] (sub-panels (a) and (b), respectively) including T-|-(yf — ;■ 6-1-6 and T-|-i — ;■ h^h^i 
with in-medium binding energies, as=0.5 and two different values for thermal parton masses. 



as=0.5, T=250MeV and thermal parton (Debye) masses viiq^g = 400 MeV, the total dissociation width 
of a weakly bound J/ip has been found to be F j/^ ~ 350 MeV. This is in reasonable agreement with 
Refs. [HEg where, using a,~0.25, Fj/^(T=250MeV)~80MeV, since the NLO rate is proportional 



to aj (modulo screening- mass corrections in t-channel gluon exchange diagrams for forward scattering; 
e.g., at T=250MeV, the Debye mass in Refs. [^ 192] is /i£)=440MeV). For a constant Debye mass, 

[95] exhibit a rather strong increase of Fj/^ with temperature (upper left 



the NLO results of Refs. 

panel in Fig. [18]), which is, however, tamed if a perturbative temperature dependence, fin oc g{T)T, 
were included (since the rates decrease with increasing /xd). The dependence of the NLO rate on the 
3-momentum, p, of the J ftp is rather flat, see lower left panel in Fig. [18] [96]. This is quite different from 
the LO gluo-dissociation process whose rate drops substantially with p (even more so for weak binding 
for which the cross section is concentrated at small gluon energies), see dashed lines in the right panels 
of Fig. [18] The p dependence following from the quasifree dissociation rate [97j (solid lines in the right 
panels of Fig. [T8]l is quite similar to the full NLO calculation (recall that q;^ differs by a factor of ~2 in 
these two calculations). 

Inelastic collision rates of the ground-state T are compiled in Fig. [19] Analogous to charmonia, 
in-medium reduced binding energies render the quasifree process dominant over gluo-dissociation, see 
Ref. [ni] for a detailed discussioiu; this has also been confirmed in a complete NLO calculation 



Since the vacuum binding energy of the T is quite large (e^— 1-1 GeV) compared to typical temperatures 
at RHIC (T~0.3GeV), its dissociation rates are rather sensitive to color-screening (which has important 
consequences for T suppression in heavy-ion collisions), e.g., ~7MeV for £:y~1.1 GeV vs. ~80MeV for 
Ey^'^ ^0.25 GeV at T=0.3GeV. For comparable values of the parameters (strong coupling constant, 
binding energy as well as Debye and thermal masses), the widths within the quasifree calculation [99] 
{eB{T)-cuTves in the left panel of Fig. [T^ tend to be larger than within the NLO analysis [SS] (right 
panel of Fig. [T^ when accounting for the factor 2 difference in a^. This discrepancy could be due a 
reduced accuracy of the quasifree approximation for more tightly bound systems [HI [82] . As discussed 
in the latter paper, in perturbation theory, the imaginary part of the finite-temperature heavy-fermion 



^For the vacuum binding energy, the ghio-dissociation rate of T is comparable to the quasifree one. 




Figure 20: S'-wave chamionium correlator ratios obtained within the T-matrix approach [TT] using 2 
different values for the in-medium charm-quark width corresponding to total r]c (or J/ip) widths of 
40MeV (lower squares) and 100 MeV (diamonds); the reconstructed correlator used for normalization 
is the one at 1.5 T^, and the charm-quark mass has been held fixed at mc=1.7GeV (no zero-mode 
contributions are included). 



potential referred to in the previous section is closely related to quasifree dissociation. In particular, it 
has been shown that in a QED plasma the infinite-distance limit of the imaginary part, ImVQQ(r — )■ 
oo) = —q'^T/Att, precisely corresponds to (twice) the damping factor of a single heavy fermion in the 
static limit, Tq = —g'^T/Sn (in pQCD, an additional Casimir factor, Ci?=4/3, appears). For small 
distances (i.e., for tightly bound states), interferences occur which suppress this quarkonium width 
from the quasifree limit. 

Recent calculations have investigated the effects of anisotropic momentum distributions in the QGP 
on the quarkonium width |100l llOli 1102] . The emerging consensus seems to be that anisotropics do not 
have a large impact on the quarkonium dissociation temperatures and widths as long as the thermal 
parton density (or entropy density) is kept constant. This conclusion may depend on the microscopic 
mechanism underlying the dissociation. For example, if the dissociation rate increases (decreases) with 
3-momentum, harder momentum spectra of the medium partons may lead to an increase (decrease) in 
width over the isotropic thermal case. 

Finally, let us address the impact of finite- width effects on the charmonium correlator ratios discussed 
in the previous Section. In the T-matrix approach of Ref. [71], a J/ip width has been implemented 
via an elastic charm-quark width, Fc, in the two-particle propagator, Gqq, figuring into the Lippmann- 
Schwinger Eq. ( TT71) (see also Ref. |103] for a more recent calculation). As discussed above, this generates 
a charmonium width in the spirit of the quasifree dissociation mechanism where the "elastic" scattering 
of thermal partons essentially occurs with one of the charm quarks within the bound state. Using 
Fc = Fg = 50 MeV, the J/ip spectral function acquires a width of about 100 MeV, and the corresponding 
correlator ratio with this width is compared to the calculations shown in Fig. [12] (corresponding to a 
total J/ip width of ~40MeV), using the correlator at 1.5 Tc as the reconstructed one. One finds that 
the moderate increase of the J/ip width increases the large-r regime of the correlator ratio by a few 
percent, due to the shght increase of strength in the low-energy part of the spectral function caused 
by the broadening. More systematic investigations of width effects, especially close to the "melting" 
temperature, remain to be carried out, e.g., by implementing microscopic dissociation width calculations 
into the T-matrix approach. Generally, finite-width effects tend to accelerate the dissociation of the 
bound state, i.e., lower its dissociation temperature jl03] . 
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Figure 21: Diagrammatic representation of the J /ip self-energy due to interactions with pions from the 
hadronic heat bath, with effective vertices in a hadronic basis (left panel) and resolved into a quark 
basis with microscopic pion (G^g), -D-meson [Gqq) and J jip {Gqq) propagators as well quark-exchange 
interactions (f/ex (right panel) [91]. The J/ifj dissociation width follows from the imaginary part of 
the self-energy, r,i, = — ImS^/m<i»; the underlying process may be visualized by cutting the left-hand 
diagram through the middle (note that this cannot be done at the quark level since confinement forbids 
external quark lines in the hadronic phase). Diagrams of this type naturally appear in the cluster 
expansion for two-particle properties, see Fig. [29] in Sec. 12.51 



2.4 Quarkonium Dissociation in Hadronic Matter 



The theoretical description of the J /iIj in medium is not complete without an understanding of its 
properties in hadronic matter. Since the latter is characterized by confined color charges, color screening 
of a tightly bound, spatially compact cc state is beheved to be negligible below Tc (on the order of 10 MeV 
or so) [104J. Thus, most of the calculations of J /ip properties in hadronic matter have focused on the role 
of dissociation reactions of the type J /ip + M ^ D + D where M denotes a meson from the heat bath 
(e.g., 71, p) and D generically denotes a D-meson (e.g., D(1870) or L'*(2010)); for reaction with baryons 
[B) from the medium one has J /ip + B ^ Ac, Ec + D. On the phenomenological side, a quantitative 
assessment of the inelastic reaction rates in the hadronic phase of a heavy- ion reaction is indispensable for 
disentangling the "anomalous suppression" in the QGP. Here, we focus on the most recent developments 
over the last 5 years or so (for earlier reviews, see, e.g., Refs. |105[ I106[ I173[ I107[ 1108] ). 

The "natural" framework to evaluate hadronic J/ip dissociation are effective hadronic models to 
compute the pertinent inelastic cross section for the reactions of the type listed above. However, even 
after the exploitation of standard symmetry in hadronic physics (such as vector-current conservation, 
vector-meson dominance |113] and chiral [111] 1112] or broken flavor- S'f/(4) |113] 1114] symmetry), a 
quantitative control over the dissociation rates could not be achieved due to uncertainties in hadronic 
formfactors. The latter are required to account for the finite size of the vertices in the effective theory. 
When varying the pertinent cutoff parameters over a typical hadronic scale, say, A=l-2GeV, the J/ip 
dissociation cross sections can vary by about 2 orders of magnitude (from sub-mb to tens of mb). 
Alternatively, parton-based calculations have been performed which, in principle, could provide a more 
microscopic description of the effective vertices. Combining the LO gluo-dissociation cross section (as 
discussed in the previous section for deconfined gluons in the QGP) [88] with the gluon distributions 
within a hadron, the resulting hadron-J/ip break-up cross section was found to be small (sub-mb) [145] , 
primarily due to the soft gluon distribution function within the hadrons (also, the thermal motion of 
hadrons below Tc is suppressed compared to deconfined gluons above T^, medium effects in the bound 
state have been studied in Ref. [90], and NLO calculations have been carried out for the T |115] but not 
for the J/ip). However, the hadronic dissociation reactions of charmonia are presumably dominated by 
other processes, specifically quark exchanges as illustrated in Fig. ^T\ which have been addressed within 
effective quark models |117[ 1116^ 1120^ 1106^ 1107^ 1112] . A realistic description of such processes requires 
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Figure 22: Dissociation cross sections of the J/ip induced by pions (left panel) and p-mesons (right 
panel) as a function of center of mass energy of the reaction. In each panel, results for a chiral 
(Nambu-Jona-Lasinio) quark model (labeled "MI") |112j are compared to those for a chiral effective 
hadronic model [109] in which the vertex formfactor cutoffs (quoted in parenthesis) have been adjusted to 
approximately match the strength in the quark model calculations. The low-energy thresholds in the left 
panel correspond to (from left to right) 2miy ~ 3.75 GeV, mjy + m£,* ~ 3.88 GeV and 2m£,* ~ 4.02 GeV, 
and in the right panel to (from left to right) rrij/^ + rUp ~ 3.87 GeV and 2m/). ~ 4.02 GeV. Figures are 
taken from Ref. |112] . 



non-perturbative input in terms of both the effective quark propagators (especially for their momentum- 
dependent mass functions in the chirally broken and confined hadronic phase) and the interaction 
vertices (which are usually resummed in ladder approximation). In principle, these quantities can be 
constrained by vacuum phenomenology (hadron spectra, decay constants and quark/gluon condensates), 
but their implementation in a hadronic medium poses formidable challenges, including a fully relativistic 
treatment, cancellations due to constraints from chiral symmetry, etc. Note that the quarkonium 
propagator, Gqq, figuring into the right- hand- side of Fig. |2T1 is closely related to the potential models 
discussed in Sec. 12.21 For example, the DD state appropriate for hadronic matter could be implemented 
into the charmonium T-matrix by coupling an additional DD channel, Gj^q, into Eq. (IT7|) using a 
suitably constructed transition potential, [/ex = Vcc^dd- 

A recent example for the calibration of hadronic vertices evaluated in chiral effective theory |109] 
using a chiral quark model |112] is shown in Fig. [22] for the processes vr + J/ip -^ D + D (left panel) and 
p + J/ijj -^ D + D (right panel; see also Ref. [121] for a similar analysis). These plots illustrate several 
features: (a) The energy dependence of the cross sections for hadronic and quark models does not match 
well; most notably, the quark-model results exhibit threshold peaks which do not appear in the hadronic 
calculations (whether this is a realistic feature remains to be seen); nevertheless, with formfactor cutoff 
values in a typical hadronic range, At,j/^=1 GeV and Apj/^=1 GeV, the overall strengths are comparable; 
(b) the cross sections for p-induced dissociation are significantly larger than the pion-induced ones, and 
both agree within a factor of ~2 with their counterparts in quark-model calculations [105[ I106[ 1107] . 
This is quite encouraging in view of the much larger uncertainties alluded to above; similar results are 
obtained in Ref. [121] . A recent comparison for nucleon-induced dissociation can be found in Refs. [118[ 
1110] (note that the low-energy "^-N break-up cross section, which is relevant for J/ip absorption in an 
equilibrated hadron gas, is different from (and probably much smaller than) the high-energy absorption 
cross section figuring into "nuclear absorption" in heavy-ion collisions discussed in Sec. 
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Figure 23: Left panel: energy (s) and temperature (T) dependence of the effective cross section (a*) 
for J/ip breakup by p-meson impact. We display p^a*{s\T) {p: initial-state hadron momentum) for 
better visibility of the effective lowering of the breakup threshold when T exceeds the D-meson Mott 
temperature T^°** ^ 172 MeV; right panel: T dependence of the thermally averaged J/ip breakup cross 
section in a n-p gas; the calculation with vacuum D-mesons (dashed line) is compared to one with in- 
medium broadened D-meson spectral functions (due to the Mott effect at the chiral phase transition) 
causing a step-like enhancement (solid line) caused by the reduced breakup threshold; from Ref. [IPS] . 



The energy-dependent cross sections for quarkonium break-up by hadron impact enable to calculate 
the temperature- (and density-) dependent dissociation rates in hadronic matter. If one assumes the 
short-distance vertex functions to be unaffected by the surrounding medium, the issue remains whether 
mass and width of the final-state open-charm hadrons change with temperature and density. These 
medium effects, in particular, imply modifications of the thresholds for the breakup processes [38| 
I123[ I124[ I125[ I126[ I127[ 1128] . The theoretical basis for the discussion of quark-exchange effects in 
the self-energy of heavy quarkonia in strongly correlated quark matter is given by systematic cluster 
expansion techniques developed in the context of plasma physics (more details are contained in the 
next subsection in connection with Figs. [261 123 and l29l) . Fig. [23] illustrates the effect of the spectral 
broadening of D-mesons at the chiral phase transition due to the opening of the decay channel into 
their quark constituents (Mott effect) for temperatures exceeding the D-meson Mott temperature, 
j^Mott ^ 172 MeV. The temperature-dependent lowering of the J/ip breakup threshold in the thermally 
averaged J ftp cross section in a vr-p meson gas induces a step-like enhancement [105]. This effect has 
been discussed as a possible mechanism underlying the threshold-like anomalous suppression pattern of 
J/f/^'s observed by the NA50 experiment |125] 1129] . 

We finally make a rough estimate of the magnitude of the dissociation widths of the J/ip in hadronic 
matter, based on the most recent evaluations discussed above, and compare it to the values found in 
the QGP as discussed in the previous Section. In analogy to Eq. ( 123|) . the J/ip dissociation rate in 
the hadron gas can be written as a convolution of the h-J/ip break-up cross section with the thermal 
distributions of all hadrons h in the heat bath, 
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Let us approximate this expression by T,^ = {a'^^^^ n)!^^ Vrei) with a total density n)^^ = Skq = 0.8 fm~^ 
for a hadron-resonance gas at T=175 MeV. From Figs. [22] and [23] we estimate (cr'^'^^ f rei) — 0.5 mb as an 
estimate for the average over all hadrons (note that mass thresholds, the relative velocity fj-ei < 1 and 
an average over the thermal motion significantly reduce the average value for a'^^^^ fj-eb compared to the 
relatively narrow maxima in the cross section plots), yielding F^^ ~ 8 MeV. This is significantly smaller 



than the estimates for the QGP-induced suppression shown in the right panel of Fig. [T71 even close to 
Tc. We can therefore expect that the anomalous suppression of the J/ip in heavy-ion collisions is very 
small (and even less for the T). However, the situation may be quite different for excited charmonia, tp' 
and Xc, which contribute "indirectly" to J ftp production via late decays after freezeout ("feed-down"). 
Unfortunately, much less is known about hadronic dissociation of ip' and Xc states. 

2.5 Insights from Plasma Physics 

In developing a theoretical approach to heavy quarkonia as messengers of the deconfinement and/or 
hadronization transition of a quark-gluon plasma as formed in a heavy-ion collision, one should aim 
at a unifying description where hadrons appear as bound states (clusters) of quarks and gluons. The 
situation is reminiscent of the problem of two-particle states in QED plasmas where a well-developed 
theory in the framework of the Green function technique exists. These methods have been widely 
elaborated for the case of the hydrogen plasma, where the electrons and protons as the elementary 
constituents can form hydrogen atoms as bound states of the attractive Coulomb interaction. The 
problem is tractable analytically for the isolated two-particle system, with a discrete energy spectrum 
of bound states and a continuous spectrum of scattering states. More complex bound states, such as 
molecular hydrogen can also be formed. 

In a many-particle system, the problem of bound state formation needs to account for medium 
effects. They give contributions to a plasma Hamiltonian, 

^pl _ ^Hartree , ^Fock , ^Pauli , jjMW , ^Debye , ^pp , jjvdW , j'2y\ 

where the first three terms - the Hartree and Fock energies of one-particle states and the Pauli blocking 
for the two-particle states - are of first order in the interaction and represent the mean-field approxima- 
tion. The following two terms of the plasma Hamiltonian are the MontroU-Ward (MW) term accounting 
for the dynamical screening of the interaction in the self-energy, and the dynamical screening (Debye) of 
the interaction between the bound particles. These contributions are related to the polarization function 
and are of particular interest for plasmas due to the long-range character of the Coulomb interaction. 
In a consistent description, both terms should be treated simultaneously. The last two contributions to 
the plasma Hamiltonian are of second order in the fugacity (particle density): the polarization potential 
(pp), describing the interaction of a bound state with free charge carriers, and the van der Waals (vdW) 
interaction, accounting for the impact of correlations (including bound states) in the medium on the 
two-particle system under consideration [TSl [72] . 

Approximations to medium effects in the self-energy and the effective interaction kernel have to 
be made in a consistent way, resulting in predictions for the modification of one- and two-particle 
states. On this basis, the kinetics of bound-state formation and breakup processes can be described, 
establishing the ionization equilibrium under given thermodynamic conditions |13U] . see, e.g., Eq. (122]) . 
Coulomb systems similar to the hydrogen plasma are electron-hole plasmas in semiconductors |131j . 
where excitons and bi-excitons play the role of the atoms and molecules. Other systems which have 
been widely studied are expanded fluids like alkali plasmas or noble gas plasmas [18|. Applications of 
the plasma physics concepts for cluster formation and Mott effect to the rather short-ranged strong 
interactions have been given, e.g., in Refs. |132[I133] for nuclear matter and in Refs. [13411135] for quark 
matter. 

In this Section, we would like to discuss basic insights from the investigations of bound-state for- 
mation in electromagnetic plasmas, as far as they might concern our discussion of heavy-quarkonia 
formation in hot and dense matter. Before going more into the details, let us give a brief overview. The 
bound-state energy remains rather inert to changes of the medium since the self-energy and interaction 
effects partially compensate each other to lowest order in density. The small size of the bound states 
plays an important role in this respect. The compensation is not operative for continuum states which 



are influenced by self-energy effects only, usually resulting in a lowering of the in-medium ionization 
threshold. This often leads to a strong enhancement of the transition rates from bound to free states, 
inducing a sequential "melting" of different bound-state excitation levels into the continuum of scat- 
tering states for certain critical plasma parameters (Mott effect [I7j), until eventually the ground state 
becomes unbound. 

The theory of plasma correlations has also been developed for strongly non-ideal systems, where 
the formation of clusters in the medium needs to be taken into account. This situation is reminiscent 
of a hadronizing quark-gluon plasma; we will therefore refer to cluster expansion techniques as the 
theoretical basis. 

Bethe-Salpeter equation and plasma Hamiltonian. A systematic approach to the description 
of bound states in plasmas starts from the Bethe-Salpeter equation (BSE) for the thermodynamic 
(Matsubara-) two-particle Green function of particles a and h (in analogy to Eq. ( TT8l) in Sec. 12. 2p . 

Gab = G^f^ + G^f, Kab Gab = ^afc + ^afc "^ab G^f, , (28) 

which is equivalent to the use of the two-particle T- matrix Tab- The 2-body equation has to be solved 
in conjunction with the Dyson equation for the full one-particle Green function, 

Ga = G^ + G^TjaGa , (29) 

defined by the dynamical self-energy Sa(p, w) and the free one-particle Green function G^{p,uj) = 
[u — Eaip)]^^ for a particle of species a with the dispersion relation ea{p) = Jp^ + w^a ~ raa + p^ / {2nia) , 
see Fig. El 
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Figure 24: The coupled one- and two-particle problem in the medium. Dyson equation (left) and Bethe- 
Salpeter equation (right) need to be solved in consistent (conserving) approximations for self-energy 
(E) and interaction kernel {K). 

The BSE contains the information about the spectrum of two-particle bound as well as scattering 
states in the plasma. A proper formulation of the plasma effects on the two-particle spectrum is essential 
to understand why bound and scattering states are influenced in a different way by the surrounding 
medium, leading to the Mott-effect for bound states. We here give the essence of a detailed discussion 
as presented in Ref. [72] . 

The homogeneous BSE associated with Eq. fl28l) can be cast into a form of an effective Schrodinger 
equation for the wave function 4'ab{pi,P2, z) of two-particle states in the medium |131] . 

H {[^a(pi) + £fe(p2) - Z] (5g,o - Vab{q)] ^Pab{pi + Q, P2 - Q, z) = J^ Hal{PuP2, Q, z)^ab{pi + ?, P2 - Q, z),{30) 

where a, b denote a pair of particles with 3-momenta pi and p2 which transfer a 3-momentum 2g in their 
free-space interaction Vab{q), and z is a complex two-particle energy variable. The in-medium effects 
described by Eq. (l30|) have been singled out in the definition of a plasma Hamiltonian, containing all 
modifications beyond the two-body problem in free space (721 1131] , 

Hiip,,p2,q,z) = Vab{q)[Nab{pi,P2)-l] -T.Vabiq')[NabiPi + q',P2-q')-l]5,,o 

' ^ ' q' 

(i) Pauli blocking ^ ^ 



(ii) Exchange self— energy 

+ AVabipi, P2, q, z)Nab{pi, P2) -J^ ^^'^b{pi, P2, q' , z)Nab{pi + q , P2 -q')Sq.o .(31) 

^ • ' q' 

(iii) Dynamically screened potential ^ y 



(iv) Dynamical self— energy 



Here, AVab{pi,P2,fl,z) = Kab{pi,P2,(l,z) — Vab{q) stands for the in-medium modification of the bare 
interaction potential to a dynamically screened interaction kernel Kab{pi,P2, Qi z). The effects of phase 
space occupation are encoded in the function Nab{pi,P2), which for the case of an uncorrelated fermionic 
medium takes the form of the Pauli blocking factor Nab{pi,P2) = 1 — fa{pi) — fb{p2), where fa{p) = 
{exp[(£a(p) — Ha) IT] + 1}~^ is the Fermi distribution and Ha the chemical potential of species a; cf. the 
factors f^^ in Eq. ( TT71) . Eq. ( l30l) is a generalization of the two-particle Schrodinger equation, where the 
left-hand side describes the isolated two-particle problem while many-body effects due to the surround- 
ing medium are given on the right-hand side. The in-medium effects denoted in the plasma Hamiltonian 
(!27ll can be obtained from the one derived in the Bethe-Salpeter approach (I3T1) upon proper choice of 
the interaction kernel Kab so that Eq. (1271) appears as a special case of Eq. (l3Tll . 

The influence of the plasma Hamiltonian on the spectrum of bound and scattering states can be 
qualitatively discussed in perturbation theory. Since bound states are localized in coordinate space, 
their momentum-space wave functions extend over a flnite range A and we may assume them to be 
g- independent, ipabipi + (1,P2 — q, z = Eni) ~ 4'ab{pi,P2, z = Eni), for small momentum transfer, g < A, 
and to vanish otherwise [Eni energy of the bound state with quantum numbers n and /). Assuming 
further a flat momentum dependence of the Pauli blocking factors Nabipi + <1,P2 — l) ~ Nab{pi,P2) for 
small q where the interaction is strong, we obtain a cancellation of the Pauli blocking term (i) by the 
exchange self-energy (ii), and of the dynamically screened potential (iii) by the dynamical self-energy 
(iv). Therefore, the bound-state energy remains largely unmodifled by medium effects. For scattering 
states which are extended in coordinate space and can be represented by a 5-function in momentum 
space, the above cancellations do not apply and a shift of the two-particle continuum threshold results. 
For this mechanism to work it is important that approximation schemes for the self-energy and the 
interaction kernel are consistent as, e.g., in the conserving scheme of ^-derivable theories |136j . 

Summarizing the discussion of the plasma Hamiltonian, we state that bound-state energies remain 
unshifted to lowest order in the charge carrier density while the threshold for the continuum of scattering 
states is lowered. This formally implies a reduction of the binding energy if the latter is identifled with 
the difference between the 2-particle threshold and the bound-state energy, Eb = ^rria + m^ — Eni. The 
intersection points of bound state energies and continuum threshold deflne the Mott densities (and 
temperatures) for bound-state dissociation. 

When applying this approach to heavy quarkonia in a medium where heavy quarks (either free or 
bound in heavy hadrons) are rare, then Nab = 1 so that both, (i) and (ii), can be safely neglected. 
The effects (iii) and (iv) originate from the dynamical coupling of the two-particle state to collective 
excitations (plasmons) in the medium. In the screened-potential approximation, the interaction kernel 
is represented by 

Vab{PlP2,q,Uj) = V^f,{q,Uj)6p,p,+pj2g,p^-p2 

Vf,{q,uj) = Vab{q) + Vab{q)Ilab{q,oo)V^b{q,oo) = Vab{q)[l-Ilab{q,oo)Vab{q)]-\ (32) 

with the total momentum P, momentum transfer 2g and energy transfer u in the two-particle system. 
The most frequently used approximation for the polarization function, Ilab{q,u}), or for the equivalent 
dielectric function, eab{q,<-^) = 1 — ^ab{q,^)Vab{q), is the random phase approximation (RPA). In the 
next two paragraphs we discuss the static, long wavelength limit of the RPA and its generalization in 
a clustered medium. 

Example 1: Statically screened Coulomb potential. The systematic account of the modiflcation 
of the interaction potential between charged particles a and b through polarization of the medium is 
taken into account in the dynamical polarization function, Ilab{q,(^), which in RPA reads 



'''' ^^' ' ' (27r)^ E^-E^_g-u ^ ' 



For the Coulomb interaction, corresponding to the exchange of a massless vector boson, the potential is 
obtained from the longitudinal propagator in Coulomb gauge, Vab{q) = eae^/q^. For a recent discussion 
in the context of heavy-quarkonium correlators and potentials see, e.g., Refs. |811 [821 159]. Due to the 
large masses of the constituents in the heavy-quarkonium system, one may use a Born-Oppenheimer 
expansion to replace the dynamically screened interaction by its static (w = 0) and long-wavelength 
(g — !■ 0) limit. For nondegenerate systems the distribution functions are Boltzmann distributions and 
their difference can be expanded as 

fa{E;) - UE;_^) = e-^.VT (i _ ,-iE;_,-E;)/T^ ^ -fa{E;){E; - E;__^)/T , (34) 

so that the energy denominator gets compensated and the polarization function becomes 



nsr"(g, 0) = -2^ J j^JaiE;) = -Sa,^ . (35) 

The corresponding dielectric function e^^^{q,ijj) takes the form 

lim£^P^(g,0) = l + 4, ^^l = ^J:^lMT). (36) 

The screened Coulomb potential in this approximation is therefore V^Kg) = Kb(g)/£^^'^(g, 0) = eaeb/{q^-\- 
lji\)). In this "classical" example of the statically screened Coulomb interaction, the contribution to the 
plasma Hamiltonian is real and in coordinate representation given by 

AV;fe(r) = --(e-'^-^ - 1) ^ a^iD - ^^^lr , (37) 

r 2 

where a = e^/(47r) is the fine structure constant. For the change in the Hartree self-energy of one- 
particle states due to Debye screening we can perform an estimate in momentum space. 
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which compensates to lowest order in the density the shift of the bound-state energy levels due to the 
screening of the interaction ( 137]) , 



{aBfl'- Bohr radius in vacuum) so that [721 1137] 

AE^i ^ --ixl{r)ni = 0{nal,) . (40) 

The Debye mass, fi£,, is equivalent to the inverse of the Debye radius r^ characterizing the effective 
range of the interaction, and depends on the square root of the density n{T) of charge carriers. It is this 
different response of bound states and scattering continuum to an increase of density and temperature 
in the medium which leads to the Mott effect for electrons in an insulator pTl IT8] : bound states of 
the Debye potential can only exist when the Debye radius is larger than rD,Mott = 0.84 ae.o |138] . 
This entails that above a certain density even the ground-state electrons become unbound and form 
a conduction band, resulting in an insulator-metal transition (also called Mott-transition). Further 
details concerning this example can be found in Ref. |139j . 

In complete analogy to the electronic Mott effect it is expected that in hadronic matter at high 
density the hadrons - as bound states of quarks - undergo a Mott transition which results in a 



phase transition from the color-insulating phase (hadronic matter) to a color- conducting or even color- 
superconducting phase (deconfined quark matter). This applies to light hadrons as well as to heavy 
quarkonia, whereby due to the different scales of Bohr radii the Mott dissociation of heavy quarkonia 
is expected to occur at higher densities than for light hadrons. In some approaches quark self-energy 
effects are neglected and one is only left with the medium effect due to a statically screened potential. 
Consequently, in such a picture the continuum edge of the scattering states remains unshifted and, due 
to the lack of a compensating effect, the effective interaction entails an appreciable medium dependence 
of the bound-state energies (masses). For the electron- hole plasma in highly excited semiconductors 
it could be shown experimentally, however, that the compensation picture is correct and the bound 
state energies remain almost unshifted |140] . In the context of heavy quarkonia in the QGP this could, 
at least qualitatively, provide a natural explanation of the approximate constancy of the Euclidean 
correlator ratios found in thermal lattice QCD. 

One may of course absorb the self-energy effects into a redefinition of the effective interaction, by 
adding a homogeneous mean-field contribution. This is equivalent to the use of the Ecker-Weitzel 
potential |141] . 

^Ecker-Weitzel (?") = ^~^°'' " tt/^D • (41) 

r 
Recent investigations of the screening problem in the context of Debye-Hiickel theory |142 j and QQ 
correlators [59| 182] have obtained this continuum shift {—afio) as a homogeneous background field 
contribution. According to the above lesson from plasma physics, however, this contribution should be 
attributed to the self-energy of the constituents rather than to the interaction kernel, since it determines 
the shift of the continuum edge (cf. Eq. f l2Tl) ). 



Example 2: Heavy quarkonia in a relativistic quark plasma. We consider the problem of a 
heavy quark-antiquark pair interacting via (an ansatz for) the statically screened Cornell potential ( !T2l) 
in a relativistic quark plasma. The thermodynamics of this medium is described in good agreement 
with lattice QCD data by a Nambu-Jona-Lasinio model coupled to a Polyakov-loop potential (PNJL 
model |143] which allows us to estimate the Debye mass, ficiT), by evaluating the RPA polarization 
function, Eq. fl33|) . for A^^ colors and Nf flavors of massless quarks with E^ = \p\ as 

TtRPA^ ^nn\ OA f J^^VaU^ f°° dp p (5^6 2 f.o\ 

n- (^^o,o) = 25.,y^^^^ = -25.,y^ -^MP) = -^,m)T . (42) 

Here, J($) = (12/7r2) J^ dxxU{x) and /$(x) = [$(1 + 2e-^)e-^ + e-3^]/[l + 3$(1 + e-^)e-^ + e'^^'] is 
the generalized quark distribution function |144j which in the case of a deconfined medium ($ = 1) is 
just the Fermi function, yielding /(I) = 1. In the deconfinement transition region, where < $ < 1, 
quark excitations are strongly suppressed, e.g., /(O) = 1/9. Taking as the bare potential a color singlet 
one-gluon exchange V{q) = —ATcas/q"^, a^ = fi'^/(4vr), the Fourier transform of the Debye potential, 
V^{r) = —aexp{—iij:){T)r)/r, results as a statically screened potential with the Debye mass f^oiT) = 
4:'n'asI{^)T'^- The Hartree self-energy for heavy quarks in a PNJL quark plasma is then given by 
Eq. dMD. 

In the spirit of Eq. ( |38|) . one can evaluate the Hartree shift due to the screened confinement part, 
V^o„f(r) = (cr//iD)(l — exp(— /Z£)r)), of the potential ( 1T^ and show that it vanishes, Sconf = 0. 

In Fig. 125] we show the temperature dependence of the two-particle energies for heavy quarkonia as 
a solution of the effective Schrodinger equation 

H^\r; T)(t>ni{r- T) = Eni{T)<f)ni{r; T) (43) 

for the plasma Hamiltonian 

m\r- T) = 2mQ - a^D{T) - — + ^/^^(r; T) , (44) 
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Figure 25: Two-particle energies of charmonia (left panel) and bottomonia (right panel) in a statically 
screened potential, from Ref. |145j . The full results of the PNJL model (solid lines) are compared to 
results without the coupling to the Polyakov loop potential (dashed lines). 



as evaluated recently in Ref. |145] . 

We can draw three conclusions from this example: (1) Due to the suppression of light-quark excita- 
tions by the confining effects of the Polyakov loop the lowering of the continuum edge and the merging 
of two-particle bound states with it (Mott effect) occur at higher temperatures when compared to the 
simple NJL model case; (2) The inclusion of the Hartree selfenergy for heavy quarks in the plasma 
Hamiltonian, Eq. f l4^ . leads to a compensation and even overcompensation of the upwards shift of 
quarkonia masses obtained from a solution of the Schrodinger equation with the naive Hamiltonian, 
Eq. ( fT3l) : see, e.g., upper left panel of Fig. [12] or upper right panel Fig. [14] for the effect of neglecting 
the downwards shift of the continuum due to the heavy quark self-energy on the temperature depen- 
dence of the quarkonium spectrum. The compensation of a decreasing m* and decreasing eb indeed 
play an important role in understanding the rather constant euclidean correlator ratios from lattice 
QCD. (3) From Fig. [25] one can also read off the in-medium lowering of the dissociation threshold 
k^^^ = Eb, which is the energy difference between the considered bound-state level and the continuum 
edge. This quantity is indifferent to the inclusion of a rigid self-energy shift. The lowering of fcg''^'' 
with increasing density and/or temperature leads to a strong increase in the quarkonium breakup cross 
section, Eq. f[25|) . by thermal impact, and to bound-state dissociation even before the binding energies 
vanish at the critical Mott densities/temperatures for the corresponding states. To estimate this effect 
we show in the lower panels of the graphs in Fig. [25] the thermal energy Eth = T together with the 
in-medium binding energy Eb = —eb- The temperatures for the onset of bound-state dissociation by 
thermal collisions with medium particles are obtained from the crossing points of these curves which 
are considerably lower than the naive Mott temperatures corresponding to vanishing binding energies. 

For the development of a comprehensive approach to heavy quarkonia in hadronizing hot/dense 
QCD matter another insight from plasma physics may be of relevance and will be discussed next: the 
effect of strong correlations (bound states) in the bulk of the medium. To this end, the bound states will 
be treated like a new species occurring in the system. Accordingly, additional diagrams have to be taken 
into account which arise from a cluster expansion of the interaction kernel Kat and the corresponding 
self-energy Eq, see Figs. [26l[28l In the plasma Hamiltonian, if p\ this leads to a generalization of the 
self-energy contributions (cluster Hartree- Fock approximation), the distribution functions in the Pauli- 
blocking factors and the dynamical screening (cluster- RPA). The van-der-Waals interaction in Eq. f l27|) 
appears naturally as a contribution to the cluster expansion, describing polarization effects due to bound 
states in the medium. 
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Figure 26: Cluster expansion for the interaction kernel for the two-particle problem in a strongly 
correlated medium (upper equation) and the corresponding self-energy (lower left equation) with a 
dipole ansatz for the vertex (lower right equation). 



K 



A— 1 A 

X ' — X 



+ 



-^ je- 



f. 



+ 




CV^'V^v/) 




+ ... 




I\^y 



+ 



L +t 



'v/V' 



Figure 27: Alternative way of drawing the diagrams for the cluster expansion of the interaction kernel 
and the corresponding self-energy of Fig. |2S]in a form familiar in plasma and nuclear physics; T2 denotes 
the 2-body T-matrix computed in ladder approximation. 



Two-particle states in the medium: Cluster expansion. In the vicinity of the plasma phase 
transition, correlations play an important role, and their proper accounting requires rather sophisti- 
cated theoretical methods such as cluster expansion techniques. For the problem of charmonium in 
dense hadronic matter at the deconfinement transition, i.e., in the strong-coupling case, we suggest a 
systematic Born series expansion of collisions with free and bound states in the surrounding matter so 
that all terms linear in the density of free particles and bound states are taken into account. 

We here describe the cluster expansion in terms of its diagrammatic expressions for the interaction 
kernel and the corresponding self-energy. The l*** Born approximation diagrams of this expansion are 
given in Fig. |26l see also the monograph [72] . The wavy lines denote the dynamically screened interaction 
V„|, which in a strongly correlated plasma receives contributions from the polarization of the medium 
beyond the RPA, denoted as generalized (cluster-) RPA in Fig. [2S] |146] . Bound and scattering states 
are described consistently in the two-particle T-matrices. For a generalization to higher n-particle 
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Figure 28: Left panel: The dynamically screened interaction potential, V"J(w, q) (wavy line), determined 
by the bare potential (dashed line) and the polarization function, liahi'^A)- Right panel: Cluster 
expansion for the generalized RPA, where in addition to free particles (RPA) also two-particle states 
(cluster- RPA) contribute to the polarizability of the medium |l46j . 
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Figure 29: Contributions to the dynamical self-energy of a two-particle system in a correlated medium at 
2nd -Qqtj-t^ order. Left panel: impact by two-particle states without constituent exchange (van-der-Waals 
or dipole-dipole interaction). Right panel: constituent-rearrangement collisions (string-flip process), 
from Ref. |150] . see also Ref. |151] . 



correlations, see Refs. |133[ I147[ [72] . The diagrams containing T- matrices do not contribute to the 
charmonium spectrum as long as the densities of the charmed quarks and of charmed hadrons in the 
medium are negligible. This is the situation expected for FAIR, SPS and RHIC energies, but for the 
expected copious charm production at LHC these terms may become significant. 

At the 2"^^ Born order, we distinguish two classes of collisions with light clusters (hadrons) that can 
give rise to spectral broadening of the charmonia. The first class concerns hadron impact without quark 
rearrangement inducing transitions to excited states, shown in the left panel of Fig. [291 These processes 
have been considered for charmonium-hadron interactions within the operator product expansion tech- 
niques following Peskin and Bhanot [HHl [89] , see Refs. |148l 1149] . The result is a deformation of the 
charmonium spectrum under conservation of the spectral weight integrated over all charmonia states. 
In the second class are quark rearrangement (string-flip) processes, as indicated in the right panel of 
Fig. [29l They induce transitions to open-charm hadrons which are also responsible for charmonium 
dissociation in hadronic matter, cf. Sec. 12. 4[ 



Example 3: String-flip model of charmonium dissociation. As a third example for utilizing 
insights from plasma physics we discuss charmonium dissociation within the string-flip model of quark 
matter |134[ 11351 11521 1153] . In this model string- type color interactions between quarks are saturated 
within the sphere of nearest neighbors so that in a dense system of overlapping quark-antiquark pairs 
frequent string-flip processes occur to drive the system into its minimal energy configuration at any time. 
The microscopic evaluation of those processes requires diagrams as the one shown in the right panel of 
Fig. [29] involving at least four quarks. Evaluating the imaginary part of this diagram, in connection with 
the optical theorem {a^iss oc Im T), gives access to the heavy-quarkonium dissociation rate in hadronic 
matter, cf. Eq. f[2^ in the previous section [531 When considering a heavy quark-antiquark pair in dense 



matter with negligible heavy-flavor fraction, the Pauli blocking and exchange self-energy contributions 
are negligible, but the strong correlations with light quarks of complementary color within the nearest 
neighbor sphere will result in a mean- field self-energy shift (Hartree shift A^) for all quarks |154j which 
determines the shift of the continuum edge. Because of the compensation in the Bethe-Salpeter kernel 
between the effects of screening of the interaction and self-energy shifts calculated with it (see discussion 
above), it is suggested that to lowest order the bound-state energies remain unshifted when increasing 
the temperature and/or density of the medium. In contrast to the first example of Debye screening of 
long-range Coulombic interactions, the screening mechanism in the string flip model is color saturation 
within nearest neighbors, applicable for strong, short-range interactions as appropriate for the case 
of the sQGP at RHIC or dense systems at CBM-FAIR. The resulting two-particle energy spectrum 
for charmonium and bottomonium states is discussed in Ref. [93], where the static screening picture 
is compared to the string-flip picture. Within the latter the in-medium lowering of the dissociation 
threshold k^^^^ is directly given by the behavior of the continuum edge, i.e. by the Hartree shift, as 
known for rate coefficients in strongly coupled plasmas [130j . 

3 Quarkonia in Heavy-Ion Collisions 

The utilization of heavy quarkonia as a probe of the medium created in heavy-ion reactions has to rely 
on the analysis of their final number and momentum spectra. Since the equilibration of quarkonia is 
generally not warranted, the interface between their in-medium properties and observables is usually 
based on transport approaches (Sec. 13. ip . If the bulk medium is in local thermal equilibrium, the 
in-medium properties (as elaborated in the previous Section) enter in the form of relaxation times 
(dissociation rates) and asymptotic limits for the yields (equilibrium abundances). With realistic initial 
conditions, usually inferred from p-A collisions (Sec. 13. 2p . phenomenological studies over a wide range 
of energies can be performed (Sec. 13. 3p . 

3.1 Charmonium Transport 

The evolution of the abundance and spectrum of a quarkonium state, \1/, in a hot and dense fireball can 
be tracked via its phase space distribution, f^, satisfying a transport equation of the form 

P'^d^hir, t;p) = -Up T^{f, t; p) f^{f, t;p)+Up l3^{f, t; p) . (45) 

Here, V^, denotes the dissociation width discussed in Sees. 12.31 and 12. 4[ /3^ the gain term responsible 
for regeneration (discussed in more detail below), and pQ=ujp={m% +p^)^/^ the on-shell energy of state 
\1/ implying that a well-defined quasiparticle state exists (mean-field terms have been neglected which 
is justified if mass corrections are small). The right-hand-side of this equation is the collision term 
which accounts for both dissociation and formation reactions (possibly including transitions between 
different charmonium states). Detailed balance (and the correct equilibrium limit) requires that both 
terms include the same reaction mechanisms. 

Let us first discuss the situation where the gain term in the transport equation can neglected; this 
can be realized in scenarios (or evolution phases) where either no charmonium states can be formed 
(are not supported in the medium) and/or the actual charmonium number is much larger than the 
equilibrium abundance in connection with relatively small reaction rates. Under these circumstances, 
the evolution of the charmonium number is given by a sequence of suppression stages for which the 
solution of the rate equation may be schematically written in terms of a (3-momentum dependent) 
survival probability (or suppression factor) 

S^i = 5'hg 'S'qgp Snuc 




hg' '" 



( ^^ dT 



(T)^ exp - JVqgAT)^\ exp (-nTv^abs^) • (46) 

/ \ ^0 / 

In this approximation the key ingredients are the inelastic reaction rates in hadron gas (HG) and QGP 
phases, as well as primordial nuclear absorption (which will be discussed in more detail in the following 
Section. 13.21) . If the latter can be accurately determined from p-A data, and if the hadronic suppression 
is small, S^, would provide direct information on the (temperature-dependent) inelastic reaction rate 
in the QGP (recall that the complete melting of a charmonium state (e.g. due to color screening) may 
be considered as the limit F^ — )■ oo). As discussed in the introduction, the suppression approach has 
been widely applied in the interpretation of SPS data on charmonium production in light- and heavy- 
ion reactions, encompassing a wide range of models for both the charmonium dissociation mechanisms 
(hadronic comovers [211 |22l 11551 I156[ 1157] . screening-based models [U [301 1158] . as well as thermal 
dissociation-rate calculations) and the evolution of the medium (transport, fireball and hydrodynamic 
simulations). An updated account of the state-of-affairs in this enterprise will be given in Sec. l3.3l below. 
The suppression factor, Eq. (H6|) . encodes no notion of the charmonium masses, of their equilibrium 
number, nor of any coupling to the open-charm content in the system. These aspects, however, become 
essential upon inclusion of regeneration processes, i.e., the gain term in Eq. (HSj) . as we will now discuss. 
Early evaluations of secondary charmonium production in heavy-ion collisions were focusing on 
hadronic regeneration reactions in the hadronic |159] and mixed phases |160] under LHC conditions. 
In Refs. [m HHl 1161] (see, e.g., Ref. |162] for a recent update) the statistical hadronization model has 
been invoked assuming that charmonia thermally equilibrate (and diffuse into a homogeneous distri- 
bution) upon completion of the hadronization phase transition (on the hadronic side), without further 
modifications in their abundances in the subsequent hadronic phase. In terms of the rate equation ( H5|) 
this implies that the thermal relaxation time (or inelastic width) of charmonia changes rather rapidly 
across the phase transition, which is not inconsistent with the microscopic calculations discussed in the 
previous two Sections (at least for the J/ip)- It also implies that the open-charm states forming the 
charmonia have kinetically equilibrated. This assumption is less certain, with current estimates of the 
thermal relaxation time of charm quarks of around r^*^ ~ 5-7 fm/c |163] . However, it is presently an open 
question by how much the J /ip regeneration yield is affected by deviations of charm-quark distributions 
from thermal equilibrium. E.g., when comparing the limiting cases of either primordial spectra (i.e., no 
reinteractions) or a thermal source (-|-flow), a factor of ~3 difference in the J /ip yield in central Au-Au at 
RHIC has been found in the coalescence model of Ref. |164] , while the variations in the transport model 
of Ref. |165] are smaller. In Ref. (161j it was additionally assumed that charm equilibrates chemically 
with the surrounding medium, so that the charmonium density is determined solely by its mass and 
and the ambient temperature. With a hadronization temperature of Tc = 175 MeV, a surprisingly good 
description of the J/ip-over-pion ratio at the SPS was found. In Ref. [H] the assumption of chemical 
equilibration of charm-quarks was relaxed based on theoretical expectations [HI I166[ 1160] that cc pair 
production is restricted to primordial N-N collisions and frozen thereafter (available experimental data 
at RHIC [12] and SPS [13] support this assertion). Under these conditions, charmonium production at 
the hadronization transition corresponds to a coalescence of available charm and anticharm quarks. In 
thermal equilibrium, this leads to a charmonium density 

where 7c=7c denote a (anti-) charm-quark fugacity (particle densities and numbers are denoted by n 
and A^, respectively). In relative chemical equilibrium, the fugacity is determined by matching the 
number of charm-quark carriers to the total number, Ncc-, of cc pairs in the system, 

^- = ^^°^77§4 + ^^^ ^ <^^)' ^48) 

2 W'^op) ,,e,J/^,... 



n^{mm] T, 7c) = dq,jl J 7^^/*(^p; ^) > (47) 



where Nop denotes the number of all open-charm states in a fireball of volume VpB- In a QGP, where 
the charm-quark number is carried by deconfined c-quarks, one has Nop=VFBlc'^nl^{m*,T), whereas in 
the HG, Nop=VpBlc J2a ^"^aiT, hb) with a running over all known charmed hadrons (D, Z), Ac, . . .; note 
that for a finite baryon chemical potential, /i^ > 0, A^a^ 7^ A^c implying Nd 7^ A^/5, etc.). The ratio 
of modified Bessel functions, h/Io, in Eq- ( l48l) approaches one for large Nop, but goes to O.SA'op for 
Nop -C 1. In the latter limit, /1//0 can be interpreted as the probability of producing an extra (anti-) 
charm state, i.e., it enforces exact charm conservation in the canonical ensemble (c and c can only 
be produced together) |13]. Note that the charmonium equilibrium limit does not depend on nuclear 
absorption or pre-equilibrium effects. The QGP evolution prior to hadronization is, however, imprinted 
on the collective properties of the statistically produced charmonia, i.e., their p^ spectra and elliptic 
flow |164[ I165[ 1167] . If c-c coalescence occurs off equilibrium, the pertinent charmonium spectra reflect 
on the collective properties of the charm quarks, which, in turn, may serve as an additional means to 
identify a coalescence component in the spectra |164[ 1168] . In Refs. PH [92] , a 2-component model has 
been proposed where statistical production at hadronization is combined with primordially produced 
charmonia subject to suppression in the QGP, followed by hadronic suppression of both components. 

Solutions of the transport equation (145 1) for J/ijj production in heavy-ion collisions under inclusion 
of the gain have first been evaluated in Ref. 02]; neglecting spatial dependencies in the temperature 
profile as well as in the J/ip and open-charm densities, a simplified rate equation of the form 



dNj/^ 

— = AFiyc'll-c - ADiyj/i^Ug , 



XpNcU-c - XoNji^Ug , (49) 



has been employed, very similar to Eq. (I22l) . The reactivity Xo = {vreico) is given in terms of the J/ip 
dissociation cross section, for which the gluo-dissociation process, Eq. (1M1) . has been used (and detailed 
balance for the formation reactivity, Xp). As emphasized in Ref. |12] (as well as in the statistical models 
with charm conservation [HIHS]), the J/ip formation rate depends quadratically on the number, N^c, 
of (anti-) charm quarks in the system, rendering an accurate knowledge of the charm-production cross 
section an essential ingredient for reliable predictions of charmonium regeneration. In the equilibrium 
limit, the quadratic dependence is encoded in the squared fugacity factor in the charmonium densities, 
Eq. P7|) . As elaborated in Refs. [381 1173[ 1174] . the equilibrium limit of the J/ip abundance is sensitive 
to the masses of both open- and hidden-charm states. E.g., for fixed charmonium masses, a reduced 
charm-quark mass in the QGP (or reduced charm-hadron masses in the HG) implies a reduced 7c for 
fixed A'cc, which in turn leads to a smaller A^^*^; in other words, reduced (increased) open-charm masses 
make it thermodynamically more (less) favorable to allocate a given number of cc pairs into open-charm 
states rather than into charmonia, see Fig. [301 This interplay has been quantitatively worked out in 
Ref. [38] based on the rate equation f[22|) with in-medium charm-quark and -hadron masses in QGP 
and HG phases, respectively. In particular, m*(Tc) ^ 1.6-1.7 GeV on the QGP side has been chosen to 
provide a continuous transition of jc into the hadronic phase at the (pseudo-) critical temperature (this 
range of c-quark masses is somewhat below the value of m*{Tc) ~ 1.8 GeV as inferred from the lattice 
correlator analysis using the internal energy as the heavy-quark potential, recall Fig. [131). ^^ discussed 
in Sec. 12.21 the screening of lattice-based potentials leads to reduced binding energies, while the in- 
medium change in the J/ip mass is expected to be rather small. The latter is rendered approximately 
constant by a decrease in the in- medium charm-quark mass m*, as seems to be required in potential- 
model analyses of lattice QCD correlators (recall also the paragraph prior to Eq. (l4ip ). In Ref. [97] , 
the rate-equation approach of Ref. [38j, which was restricted to inclusive charmonium production (i.e., 
for p=0), has been extended to finite 3-momentum. The suppression part has been treated explicitly 
(including spatial dependencies, such as leakage effects) in an expanding fireball background, whereas 
the momentum dependence of the regeneration part was approximated with a blast-wave model (within 
the same fireball). The latter implies the underlying charm-quark spectra to be thermalized, at least in 
the relevant regime of low momentum [1631 1175] . 
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Figure 30: Equilibrium abundance of J /ip mesons as a function of temperature in an isotropic fireball 
representative for central Pb-Pb and Au-Au collisions at SPS and RHIC, respectively [38]. While the 
variation in temperature is comparatively moderate, a large sensitivity of the equilibrium numbers with 
respect to the charm-quark (charm-hadron) masses in the QGP (HG) is found. The absolute numbers 
are based on open-charm cross sections in p-p collisions of cr^p = 570 /zb and 5.5 /xb at RHIC and SPS, 
respectively (these values are roughly consistent with PHENIX measurements [121 1169] and compilations 
for fixed target energies |171l 1172] ). extrapolated to central A-A using binary-collision scaling, and for 
a rapidity interval of 3.6 units around midrapidity (corresponding to 2 thermal fireballs). The arrows 
to the left indicate the J /if) abundance from primordial (hard) production. 



A more advanced treatment of the gain term in the rate equation (H5ll within a hydrodynamically 
evolving background medium was carried out in Ref . [165] (see also Ref . |176] ) . Assuming the prevalence 
of the 2-(-;'2 dissociation rate from gluo-dissociation, ^f -|- \& -f-)- c -|- c, the gain term takes the form 



P^iPT^n^r) 



d^k 



d^Pc 



d^Pc 



2ujp J {2n)^2oJk (27r)32wp^ (27r)32a;p 



-riPc, n, t) f{p,; n, t) W^^{s) Q{Tin, r) - T,) 

x{27i)W^'\p + k-p,-p,) ,{50) 



where the charmonium phase space has been restricted to the transverse plane at central rapidity 
{pz=0=y). With W^^{s) = a'^^^^Vj.ei{2cL!k2ujp) representing the "^ dissociation probabihty, the forma- 



tion probability, Wcc (s), is inferred from detailed balance 

energies of gluon, charmonium, charm and anticharm quark, respectively, with charm(onium) masses 



cjfc, ujp, ujp^ and Up^ denote the on-shell 



of mc=1.87GeV and mj/^=3.1GeV (corresponding to a "vacuum" binding energy of e^ =0.64 GeV). 
The 9(T — Tc) function in Eq. (!50|) restricts the inelastic processes to the QGP phase, with rc=165 MeV 
in Ref. |165] . Note that an explicit treatment of NLO dissociation processes in the regeneration term 
requires the evaluation of a 3-body initial state. The initial conditions for the charmonium distribu- 
tion function, /i'(p;r, tq), for its evolution in the thermally evolving medium are typically taken from 
a Glauber collision profile for production in primordial N-N collisions, plus the effects of nuclear ab- 
sorption. In the "pre-equilibrium" phase, i.e., for times earlier than the thermalization time, tq, of the 
fireball, the gain term is switched off. While the production of a cc pair is expected to occur on a 
short time scale, Tcc ~ l/2m^ ^ O.lfm/c, the development of the charmonium wave function in its rest 
frame requires a time duration on the order of its inverse binding energy, t^ ~ ^/sb ~ 0.3fm/c. For a 
slowly moving J/ip, this time is comparable or smaller to commonly assumed thermalization times of 
To = 0.3-1 fm/c. However, for J/tp^s at large transverse momentum, the formation time in the fireball 



frame is time-dilated by the (transverse) Lorentz factor 77- = rriT/Tn, and therefore may well reach 
into the QGP phase, presumably reducing the suppression effects |157[ 11771 1178[ I179[ llSOj 1181] (in a 
naive picture, the cc wave package is of smaller size than a fully formed bound state, thus reducing 
dissociation cross sections with thermal partons). The interplay of charmonium-formation and fireball- 
thermalization time scales will be further pronounced for in-medium reduced binding energies, as well 
as for excited charmonia (which even in vacuum have rather small binding energies). 

An alternative approach to thermal descriptions of the bulk medium are microscopic transport mod- 
els where the evolution of essentially all particle distributions is treated by numerical simulations of 
the Boltzmann equation. On the one hand, transport models have the important advantage that no 
explicit assumption on equilibration is required, which in particular eliminates the introduction of ther- 
malization and freezeout times (or conditionsjj. On the other hand, the treatment of phase transitions 
is more problematic than in thermal models, where pertinent aspects are readily encoded in the equa- 
tion of state (EoS). The main input for charmonium interactions are 202 dissociation reactions (plus 
detailed balance), i.e., gluo-dissociation {g + '^ -^ c + c) in the partonic phase and meson dissociation 
(M + \E' -^ D + D with M=7c, K, p, . . .) in the hadronic phase (sometimes augmented with elastic cross 
sections). Especially hadronic cross sections at low energy, as relevant for the typical reactions in the 
hadronic phase (corresponding to temperatures T < 200 MeV), are currently beset with rather large 
theoretical uncertainties (and, in the case of meson dissociation, not accessible experimentally), recall 
Sec. 12.41 This problem applies, of course, to both transport and thermal descriptions of the expanding 
medium. Systematic comparisons of charmonium observables within complementary descriptions of the 
bulk medium are thus valuable to extract robust information on charmonium properties in heavy-ion 
collisions. Recent calculations of charmonium production within transport models, including regenera- 
tion reactions, can be found in Ref. |182l 1183^ |TB4j, and will be reiterated in the context of heavy-ion 
data in Sec. 13. 3[ E.g., in Ref. |182j . it was found that an increase in the masses of open-charm states 
(c-quarks and D-mesons in the partonic and hadronic phase, respectively) entails an enhanced J/ijj 
production, which is consistent with the systematics of the thermal model as illustrated in Fig. |30l In 
Ref. [185J . the uncertainty in meson- induced dissociation reactions was circumvented by introducing a 
single matrix element for all dissociation reactions, augmented by spin-isospin factors and phase space 
factors (subsequently used in Refs. |183l 1184] ). In a similar spirit, the hadronic comover model jl55] 
has recently been extended to include regeneration reactions |186j . 

3.2 Nuclear Absorption and Initial-State Effects 

A commonly employed strategy to isolate cold-nuclear-matter (CNM) effects on the production of heavy 
quarkonia in heavy-ion collisions, encoded in the suppression factor S'nuc in Eq- P6|) . is to study its A- 
dependence in p-A collisions relative to p-p collisions. The pertinent effects on quarkonium production 
may be classified into three different stages: (1) initial state, prior to QQ production by parton fusion, 
due to the nuclear modification of parton distribution functions (PDFs) and/or parton energy loss; 
(2) formation stage of the quarkonium states from the initial QQ pair; and (3) final-state inelastic 
interactions leading to quarkonia absorption in the target matter. It is at present an open problem how 
to disentangle these three CNM effects from the experimental results for the A-dependence measured 
in p-A collisions. To obtain a simplified baseline for predictions and applications in heavy-ion collisions, 
it is assumed that absorption is the dominant process in all three stages and may be analyzed using 
Glauber theory in terms of a survival probability for a J/ijj produced in a p-A collision [2] 
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^Much like for charmonia, primordial particle production in hadronic transport models account for a formation time, 
Tp ~ 1/Aqcd — 0.5-1 fm/c, in the rest frame of the hadron. 
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Figure 31: Left panel: J ftp cross section in p-A collisions from three different NA50 data sets |187] 
at 400 GeV and 450 GeV together with a Glauber fit of an effective nuclear-absorption cross section, 
<7abs = 4.2 ± 0.5 mb. Middle panel: Cross sections for J/tp production in p-A collisions normalized 
to p-Be for 158 GeV (squares) and 400 GeV (circles) from the NA60 experiment [188] . The lines are 
Glauber model fits with aabs(158 GeV) = 7.6 ± 0.7(stat) ± 0.6(syst) mb and (Tabs(400 GeV) = 4.3 ± 
0.8(stat) ±0.6(syst) mb. Right panel: A-dependence of the J/ip over Drell-Yan (DY) cross section ratio, 
compared to a Glauber model calculation taking into account feed-down from higher cc states |189j . 



where b is the impact parameter and nfyf{b,z) is the density profile of the target nucleus. As a result 
a "convoluted" J/ip absorption cross section, a^hs, is extracted from experimental data for the A- 
dependence of the survival probability, Eq. (15T]) . parameterizing the combined result of all CNM effects. 
A Glauber model fit to three different NA50 data sets |187j gives an effective nuclear absorption cross 
section of (Tabs = 4.2 ± 0.5 mb (see left panel of Fig. EI]) for both 400 GeV and 450 GeV incident 
proton energies. Recently, the NA60 collaboration has published results for the A-dependence of J/ip 
production in p-A collisions at 158 GeV and 400 GeV |188j which show a larger effect at lower incident 
energy: aabs(158 GeV) = 7.6±0.7(stat)±0.6(syst) mb and aabs(400 GeV) = 4.3±0.8(stat)±0.6(syst) mb, 
see middle panel of Fig. [3TJ The origin of the increase of (Tabs towards lower collision energies is presently 
under investigation. One may speculate that an energy dependence of the J/ip breakup cross section on 
nucleons could be responsible for such a behavior (note that the high-energy absorption cross section 
relevant in the present context might be quite different from the low-energy dissociation cross section 
discussed in Sec. 12.41 as relevant for the hadronic phase of heavy- ion collisions). Another option are 
formation time effects which we will return to below. 

In extracting absorption cross sections from experiment one has to keep in mind that only about 60% 
of the observed J/ip mesons are directly produced, while about 30% (10%) arise from Xc (^') feed-down. 
If one approximates a^N geometrically, as 7rr^, with rj/^ = 0.25 fm, r^/ = 2 x rj/^ and r^^ = 1.5 x rj/^, 
and repeats the Glauber-model calculation, an equally acceptable description of the A-dependence for 
400 GeV and 450 GeV is obtained |189j , see the right panel of Fig. [311 This reasoning would point 
to a smaller absorption cross section for the J/ip on nucleons, on the order of 2 mb. Formation-time 
effects have been evaluated in Ref. [190] within a quantum mechanical description of the formation 
of charmonia from an initially small-sized "pre-meson" state. An energy- and time- dependent cross 
section has been derived. At a collision energy of ^/s = 10 GeV, the initial pre-meson starts out with 
o"abs — 3mb, evolving into (Tj/^ n = 2.8 ± 0.3 mb and ct^'at = 10.5 ± 3.6 mb at asymptotic times 
(t —7- oo). In this framework, which is in accordance with geometric scaling, a good description of the 
(xi?- dependence of the) E866 data [200] can be achieved, and an apparent discrepancy between the 
o"abs values extracted from hadro- and photo-production data can be resolved [191j . A characteristic 



1.1 


- ■pTnT4°In ° 0, i 

1 Wit t«ot 

~ . . . 1 . . . 1 . . . 1 . . . 


a 

D 














1.05 

1 

0.95 

0.9 


5- 
-5-1.3- 

1.2- 

1.1- 

1- 

0.9- 

0.8- 


□ NA50 @ 400 GgV 
* E866@800GeV 
. HERA-B @ 920 GeV 










0.85 
0.8 

0-1% 






4 


-0.2 0.2 0.4 


0.6 


0.8 6 0^5 i lis 2 2^5 


3 3.5 


4 


4.5 5 
[GeV/c] 




. HERA-B @ 920 GeV 

1 E672/E706 (o) BOO GeV 

. E771 @ 300 GeV 

■ E789 @ 300 GeV 

1 NA50@400GeV 



Figure 32: The nuclear-dependence parameter a as a function of xp (left panel) and of pr (middle 
panel), for J/ip mesons measured in p-A collisions at three different projectile energies (the symbols 
in the left panel correspond to data from HERA-B (open triangles) [201] . E866 (open squares) |200] . 
NA50 at 450 GeV (open circles) [1871 [202], NA60 at 400 GeV and 158 GeV (closed circles and squares, 
respectively) |188] and NA3 (open stars) [23] ■ A measure for pj- broadening (Cronin effect), {p"^), 



increases with nuclear target size and with collision energy (right); compilation of data in the left panel 
taken from Ref. [188j and in the middle and right panels from Ref. [189] . 



feature of the quantum mechanical treatment of charmonium formation in p-A (also applicable to AA) 
collisions [1921 11931 11941 1195] is the possibility for an oscillatory behavior of the survival probability with 
time, most pronounced for the ijj'. This is different from the classical picture of a cc state expanding in 
time [196[ll97ill98[ll99"] and could be tested via detailed kinematic (e.g., xp) dependencies of production 
cross sections [195] . 

The kinematic dependencies of CNM effects on xp and pt have been measured for charmonium pro- 
duction [23l I200[ 11871 12021 I20H 1188] and are often parametrized with the nuclear-dependence parameter, 
a{xF,pT), defined via 

v4° . (52) 
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Results from fixed-target experiments at different center-of-mass (cm) energies are compared in Fig. [321 
The values of a at xf ~ (relevant for heavy-ion collision experiments) tend to increase with collision 
energy. Since the scaling of a with px (middle panel of Fig. [32]) appears to be independent of energy, 
the increase with cm energy might originate from the increase of the average px with collision energy, 
see right panel of Fig. [321 

The two parameterizations of CNM effects, Eqs. (1521) and (15T|) . can be related to each other [Sj [2]. 
Introducing the effective path length L of the J /ip trajectory inside the target nucleus with nuclear 
matter density wq, Eq. (jSTl) can be approximated as 
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(53) 



Comparing with Eq. f[52|) one finds 



CTabs = -^ In A . 



(54) 



If one furthermore assumes a uniform density profile, one has L = StqA^^^ /A, and with A^^^ ^ InA 
(valid for large A > 50), the relation f[5^ takes the form 
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Figure 33: Left panel: nuclear modification of gluon distribution functions in a Pb nucleus, according 
to the EKS98 [2031 1204] weight functions, from Ref. ^172j . Right panel: changes of the cc and 66 cross 
sections induced by the nuclear modification of the PDFs at mid-rapidity, from Ref. |172] . 



With this "pocket formula" a range of a-values, a = 0.95 ± 0.02, translates into a range of absorption 
cross sections of o"abs = 4.0 ± 1.6 mb. For a similar example, see Ref. [2]. 

Next we turn to CNM effects in nuclear PDFs (nPDFs), //^(x, Q"^), characterized by the ratio |205] 



Rtix,Q') = ffix,Q')/f!{x,Q' 



(56) 



The left panel of Fig. 133) shows this ratio for a lead nucleus within the EKS98 parametrization [20311204] . 
for two values of the momentum transfer, Q^ = 4(50) GeV^, appropriate for cc (bb) production. Relevant 
values for the Bjorken variable, x = Mj ^^ at central rapidity (y = 0) and for charm (bottom) 
production, M ~ vtiq = 1.5(5) GeV, are in the range of 0.03 < a; < 0.3 from SPS to HERA-B energies 
(200 < i^iab < 920 GeV, corresponding to y/s = 20-40 GeV). Of particular interest are gluons, which 
are expected to dominate the production cross section for QQ pairs via the two-gluon fusion process 
(except for close to the production threshold). Typical nuclear effects are: (a) shadowing at low x, 
where i?^(x,Q^) < 1, (b) anti-shadowing at intermediate x, where i?'^(x, Q^) > 1, (c) EMC effect 
at 0.25 < X < 0.8, where i?'^(x, Q^) < 1, and (d) Fermi motion at x ~ 1, where i?^(x,Q^) > 1. 



The right panel of Fig. [3S1 indicates a 10% cc production cross section enhancement at midrapidity 
for SPS energies, while there is essentially no effect at RHIC. On the other hand, for d-Au collisions 
at backward rapidities {y = —2.0) at RHIC, a 15% antishadowing is predicted, compared to 20% 
shadowing at forward rapidities oi y = 2.0 (in the direction of the Au beam). This trend for the 
y dependence can be seen in the recent PHENIX data at ^/sNN = 200 GeV [207] . see left panel of 
Fig. EH Based on EKS98 PDFs, a best fit to the measured nuclear modification factor, RdAn, results 
in a "residual" nuclear absorption cross section of a^hs = 2.8ti;4mb. For the nDSg PDF [206j a smaller 
cross section of o"abs = 2.2li;5mb has been extracted. This is an example of how, in principle, a "de- 
convolution" of initial-state effects due to the nuclear dependence of PDFs (in particular for gluons) 
may be performed. A recent systematic study of the "residual" J/ip-imdeon absorption cross section 
was performed in Ref. [208] . After factorization of the charmonium production process, a global fit 
involving data from SPS to RHIC energies resulted in cXabs = 3.4±0.2mb when neglecting uncertainties 
in the gluon distribution functions. It was furthermore concluded that present uncertainties in the gluon 
(anti-) shadowing do not allow for a precise determination of a^hs- On the one hand, the new EPS08 
nPDFs [209] , which, for the first time, included RHIC results for forward-rapidity hadron production in 
d-Au |210] into the analysis, show a strong enhancement of gluon anti-shadowing (shadowing) in the SPS 
(LHC) domain. On the other hand, the rather mild x-dependence of J/tp production in p-A collisions at 
^JWn = 38.8 GeV in E866 [200] disfavors the nDSg, EKS98 and EPS08 parameterizations [mj. This 
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Figure 34: Left panel: PHENIX data for the nuclear modification factor, -RdAu, as a function of rapidity, 
compared to theoretical curves using the EKS shadowing model [2031 I204J for different values of the 
nuclear absorption cross section (Tabs (=o"breakup)- The best fit to these data gives o"abs = 2.8li;4mb. 
Middle panel: Raa for Au-Au collisions at mid-rapidity compared to projections of CNM effects ex- 
trapolated from the d-Au data, for EKS98 and nDSg |2U6] shadowing. Right panel: Raa at forward 
rapidity; from Ref. |2U7] . 



might suggest that the suppression observed in the BRAHMS data |210j may not be entirely attributed 
to nPDFs alone, but also to other mechanisms such as initial-state parton energy loss |212] . 

Another aspect which complicates the "de-convolution" of CNM effects is the microscopic nature of 
the J/ip production process. Recent reanalyses of the two-gluon fusion process have been performed in 
Refs. |213[ I214| I215[ [2T6] (see also Ref. |217] ). These works suggest that production processes involving 
a hard gluon recoiling off the J/ip can account for the low- and intermediate-transverse momentum 
spectra from the Fermilab Tevatron (CDF) to RHIC (PHENIX), as opposed to the soft-gluon emission 
mechanism of the color evaporation model (CEM) at leading order (LO) or to the color-octet model 
(COM) of NRQCD at LO. In this context, one distinguishes the kinematical effect of an "extrinsic" 
source of Pt due to a recoiling gluon from a "primary" source of Pt due to an "intrinsic" Pt inherited 
from the initial gluon pair (as in the COM and CEM). In the "extrinsic" case, a shift of the rapidity 
dependence of the J/ip RdAu occurs, not inconsistent with PHENIX data [207| . In connection with the 
EKS98 nPDFs this requires a larger absorption cross section, o"abs = 4.2 mb, to be compatible with the 
data fITE\ . 

It is important to assess the limits of a classical description of CNM effects, i.e., whether a quantum 
mechanical approach to the propagation of a charm quark pair (pre-meson) in the nuclear environment 
is warranted. The pre-meson evolution is governed by two time scales: the proper time for the formation 
of a spectrum of charmonium eigenstates typically estimated as tj/^ ~ (m^/ —mj/^)~^ ^ 0.35 fm/c, and 
the so-called coherence time, Tc = 2Ej/^/m'^j,^, representing the lifetime of a cc fluctuation in the lab 
frame. When the formation length becomes larger than the size of the Lorentz contracted nucleus, the 
phenomenon of color transparency occurs |219[I220] . also interpreted geometrically by a small-size pre- 
meson state |221] . When the coherence length becomes comparable to this size (e.g., by Lorentz time 
dilation), the amplitudes for processes on different nucleons interfere destructively, producing additional 
suppression. A review of the quantum mechanical treatment including applications to heavy-flavor 
production in a nuclear environment can be found in Ref. 



An example for the extrapolation of the CNM effects from d-Au to Au-Au collisions is shown in 
Fig. |3l] for mid- (middle panel) and forward rapidity (right panel) in comparison to PHENIX data 
|207j . While at forward rapidity the onset of effects beyond CNM suppression appears at A^part ~ 70 
(which roughly coincides with the onset of anomalous suppression in the NA50 experiment), there is 
almost no deviation from the CNM baseline at midrapidity, even for central collisions. This counter- 



intuitive result points to the existence of an additional J /il) production process which is predominantly 
operative in the central-rapidity region. In the previous Section, charm recombination in the medium 
has been discussed as a gain process in the kinetic equation for charmonium production. Its role in 
the interpretation of present and future experiments will be discussed quantitatively in the following 
Section. 

3.3 Applications at SPS, RHIC and LHC 

In this section we attempt to convert the theoretical and phenomenological developments described 
above (charmonium equilibrium properties as well as implementations into transport approaches) into 
interpretations of available data from ultrarelativistic heavy-ion collisions (URHICs) at SPS and RHIC 
As a rough guideline, let us start by recalling some of the basic relations which have been proposed to 
associate different observables with key properties of charmonia in medium. 

(A) Centrality dependence of J /ip production, which is typically normalized to the expected yield in 
the absence of any medium effects, -Raa (A^part) ^ based on an increasing matter density with cen- 
trality (i.e., decreasing impact parameter), this is the classic observable to search for an anomalous 
suppression of charmonia possibly linked with the onset of QGP formation. 

(B) Transverse-momentum {pt) spectra, normalized to the expected spectra without medium effects 
(i.e., in p-p collisions), Raa{pt) ^ here, the original objective is the extraction of the pt dependence 
of the charmonium dissociation rates. In addition, py-spectra from regeneration are expected to 
be softer than the primordial power-law spectra, which provides a handle to disentangle secondary 
and (suppressed) primordial production. The discrimination power may be further augmented in 
the elliptic flow, V2{pt)- 

(C) Excitation function, Raa{\/s) ^ the suppression effect on primordial charmonia is expected to 
increase with increasing collision energy, providing complementary information on the temperature 
and density dependence of charmonium disintegration (in addition to centrality dependencies at 
fixed energy). However, the possibly largest virtue of an excitation function is the stronger 
increase in charm-quark production relative to light particles, leading to a large variation in the 
charm-quark densities in the fireball which (quadratically) enhances regeneration mechanisms 
with increasing energy. The excitation function thus promises to be a prime observable to study 
the interplay of suppression and regeneration mechanisms. 

(D) Rapidity dependence, RAA{y) ^ the idea behind this observable is reminiscent of (C) in that light- 
and charm-quark production are characterized by different rapidity distributions (narrower for 
heavy quarks); regeneration, being proportional to N^^ (grand-canonical limit), should thus lead 
to a narrowing of the charmonium distributions compared to p-p collisions. 

(E) Excited charmonia {xc, V^') ^ different binding energies, dissociation rates and dissolution tem- 
peratures are hoped to provide a systematic suppression/production pattern that can serve as 
a "thermometer" upon varying control parameters such as centrality, system size and collision 
energy (all of the J/ip observables listed under (A)-(D) apply here). 

(F) Bottomonia - due to stronger binding energies, bottomonium suppression could provide a more 
sensitive measure of color screening via an associated increase in dissociation rates. This notion is 
corroborated by the small open-bottom cross sections which presumably suppresses bottomonium 
regeneration. 

We recall that feed-down contributions via the decays of excited charmonia into lower lying ones need 
to be considered. For the J/ip about 30% of its inclusive production yield in p-p collisions arises 



from decays of Xc states |223] . and ~10% from ifj' \12A\ 1225] . Such fractions cannot be neglected 
in quantitative calculations of J /ip observables in heavy-ion reactions, rendering a simultaneous and 
consistent treatment of medium effects on Xc and ip' production mandatory. At high pt>5 GeV, relevant 
for collider energies, the inclusive J/ip yield is expected to receive additional feed-down contributions 
from 5-meson decays, on the order of 20-40% |226l 1225] . 

In the remainder of this Section we discuss the 6 classes of observables in the order given above, 
focusing on the recent status and developments, and with preference for approaches which have been 
applied at both SPS and RHIC energies. 

3.3.1 Centrality Dependence 

We start our discussion with the centrality dependence of J/ip production at SPS energies, specifically 
with implementations of the J/ijj transport equation into a thermally expanding background medium. 

In Ref. |227j the gluo-dissociation rates for J/ip and Xc (using vacuum binding energies) have been 
folded over a 2-1-1 hydrodynamic evolution of Pb(158 AGeV)-Pb collisions, employing an EoS consisting 
of a 2-|-l-flavor QGP and a hadron resonance gas (connected via a mixed phase at Tc ~ 165 Me V). 
Hadronic suppression of the J/ip and regeneration effects (expected to be small at SPS energies due 
to small open charm production) have been neglected. The primordial nuclear absorption cross section 
has been fixed at o"abs = 4.3 mb. The overall magnitude of the suppression in the NA50 data has been 
reproduced by introducing a minimal initial entropy density of s™™ ^ 32fm~^, below which the matter 
is assumed to decouple instantaneously (i.e., not to become part of the thermally evolving medium); 
with a thermalization time of tq = 0.8fm/c (the moment when the hydrodynamic evolution is assumed 
to start), the centrality dependence of the J/ip yield (upper left panel of Fig. [35]), as well as its average 
momentum squared (lower left panel), are well reproduced. 

Within a similar spirit, but adopting a slightly different approach, solutions of the rate equation ( 122|) 
have been evaluated in Refs. [381 [97]. The thermal evolution of the medium has been treated in a more 
simplistic isotropic fireball expansion, parametrized to resemble hydrodynamic models, with an EoS in- 
cluding QGP, mixed and hadronic phases and T^ = 170(180) MeV at SPS (RHIC)|!|, similar to Ref. [227] . 
The inelastic charmonium reactions in the QGP include an in-medium reduction of the binding energies, 
incorporated using the quasifree dissociation process [44j (rather than gluo-dissociation, recall Fig. [T6]) 
corresponding to the rates displayed in the right panel of Fig. [TT] The dissociation rates in the hadronic 
phase have been estimated based on SU{A) effective theory, leading to rather small effects on the J/ip 
(somewhat more significant for Xc and ip'). Fixing cXabs = 4.4 mb and Tq = Ifni/c (translating into an 
average initial temperature of Tq ~ 210 MeV in central Pb-Pb), the coupling constant in the quasifree 
process has been adjusted to as—0.25 to fit the suppression in central Pb-Pb. Also here the centrality 
dependence of inclusive J/ip production and its mean (p^) are well described. The contribution from 
regeneration turns out to be small (with a total open-charm cross section of o"*g*=5.5/ib distributed over 
2 thermal fireballs, i.e., ~3.5 units in rapidity); it is further reduced by a schematic implementation 
of incomplete charm-quark thermalization via a relaxation-time factor TZ{t) (with a smaller relaxation 
time, T^'^, leading to a larger regeneration yield), see below. 

Hadronic comover-interaction [155] and transport [229] models have also been successful in describing 
the J/ip suppression pattern observed by NA50. With a typical comover dissociation cross section of 
0.65-1 mb, the required comover energy densities in central Pb-Pb are rather large, (well) above the 
typical value of 1 GeV/fm^ associated with the phase transition. The comover interaction is thus to 
be understood as at least partially of partonic origin, consistent with the rather small value of the 
dissociation cross section. In the hadronic transport calculations of Ref. [229] , the hadronic comover 



^A smaller value of Tc at SPS energies is expected due to its decrease with jib, while Tc—l80MeV at /iB=0 is within 
the current range of lattice results for full QCD [228) . 
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Figure 35: J/^jj production in Pb-Pb collisions at SPS employing the rate equation (H5l) within a 
hydrodynamic evolution with gluo-dissociation in the QGP (left panels, gain term neglected) [227] , 
and within an expanding thermal fireball with quasifree dissociation in the QGP and meson-induced 
break-up in the hadronic phase (right panels) [38l [97] . Upper panels: nuclear modification factor as a 
function of centrality quantified by the measured transverse energy, Et, of produced particles. Lower 
panels: centrality dependence of the average transverse momentum squared of the J/ip spectra. 



interactions are significantly larger (a few millibarns for baryon-induced dissociation, and 2/3 of that 
for mesons), but "pre-hadronic" degrees of freedom still play a role. 

An often discussed issue in the interpretation of the NA50 data is the (non-) presence of a threshold 
behavior in semi-central Pb-Pb collisions, around E^ ~ 30-40 GeV. To address this question experimen- 
tally, the successor experiment NA60 has measured J/ip production for an intermediate-size system, 
i.e., In-In collisions (A/„ = 115). The pertinent data are compared theoretical predictions in the left 
panel of Fig. [36], all 3 of which are in reasonable agreement with the Pb-Pb data. While the comover 
model [233j and the threshold-melting scenario [231] do not reproduce the data, the thermal rate equa- 
tion approach in QGP + hadronic phase [38] I232J shows a fair agreement (except for central collisions 
where the data suggest a reduced suppression). The suppression pattern at SPS energies has been 
revisited within the Hadron String Dynamics (HSD) transport model in Ref. [183j . where charmonium 
suppression was studied within the QGP threshold-melting scenario (middle panel of Fig. [5B]) and a 
comover-interaction model (right panel of Fig. [5B]) . The former is characterized by critical energy den- 
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Figure 36: Left panel: NA60 data for J/ijj production in In(158 AGeV)-In collisions at SPS energies |230] 
compared to theoretical predictions based on the QGP threshold-suppression scenario |231j (upper 
dashed line), QGP and hadronic suppression in a thermal fireball [232] (middle dashed-dotted line) and 
the comover model |233] (lower dotted line); all models are in fair agreement with NASO data in Pb-Pb 
collisions. Middle and right panels: postdictions for Pb-Pb |32] and In-In |234] data computed in the 
HSD transport model [183J using either a QGP threshold-melting scenario (middle panel) or comover 
model (right panel); note that the NA60 data are the preliminary ones from Ref. 



sities above which the corresponding state does not forno while the latter has been implemented via a 
single matrix element for all charmonium states (including phase space and spin effects). The threshold 
melting leads to a rather pronounced drop in J/ip production at small centralities (iVpart<75) together 
with a leveling off at large Npart, whereas the comover scenario produces a more gradual suppression 
pattern which reproduces the trend in the data somewhat better |183j . Note that the level of agree- 
ment of the HSD-|-comover scenario with NA60 data is quite comparable to the thermal rate-equation 
prediction [3H1 1232] in the left panel of Fig. 131] (the latter is sometimes quoted as not reproducing the 



NA60 data [Ml [2201 1235] )• 

Let us now turn to RHIC energies, focusing on the approaches discussed above in the SPS context. 
The hydrodynamic model with gluo-dissociation has been extended to include regeneration reactions, 
c + c-^g + '$, in Ref. [165] . The charm-quark distributions figuring into the gain term, Eq. fl50|) . have 
been approximated by two limiting scenarios: (a) primordial spectra in p-p given by perturbative QCD 
(pQCD), or (b) fully thermalized spectra including the collective flow of the hydrodynamic background. 
The corresponding difference in the inclusive J/tp yields turns out to be rather small, cf. upper left 
panels in Fig. [371 The underlying open-charm production cross section has been adopted from the 
PHENIX measurement in Ref. [12] , (jpp=622±57 /ilo. The experimental uncertainty translates into an 



uncertainty in the regeneration yield reflected by the colored bands in the left panels of Fig. [37] (the cross 
section has been scaled by the number of binary N-N collision at given impact parameter). The overall 
magnitude of preliminary PHENIX J/ijj data |236] is roughly reproduced (the agreement is somewhat 



better with updated PHENIX data [37] displayed in the upper right panel). The contribution from 
regeneration is moderate but significant, reaching ca. 50% of the total yield in central Au-Au. 

The right panels of Fig. [37] shown results of the thermal rate-equation approach [38] , extended to 
include the finite 3- momentum dependence of the inelastic charmonium reaction rates [HZ]. Here, the 
impact of incomplete charm-quark thermalization has been schematically modeled by a "thermal re- 



laxation time" factor, TZ{t) 



-t/Tc 



multiplying the charmonium equilibrium abundance. 



Nl\ 



which figures into the gain term of the rate equation (!22|) . This implementation leads to a stronger 



"'^"it may be thought of as a sudden transition of the quarkonium width from zero to infinity. 

^^An update of the PHENIX open-charm cross section amounts to gpp=567ib57 ^b fl69| : for both measurements the ad- 
ditional systematic uncertainty is about 150-200 /ib; STAR data [170] give significantly larger values, (Tpp=1.3±0.1±0.2mb, 



which is a currently unresolved discrepancy. 
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Figure 37: J/V' production in Au-Au collisions at RHIC computed within kinetic rate equations in 
a thermally evolving background medium, corresponding to the approaches displayed in Fig. |35] (left 
panels: gluo-dissociation with vacuum binding energies in a hydrodynamic simulation |165[ 1227] . right 
panels: quasifree dissociation with in- medium binding energies in a thermal fireball expansion [38|l97]). 
The upper panels show the nuclear modification factor, R^a impart) vs. centrahty for inclusive J/ifj 
yields, while the lower panels display the centrality dependence of the J/ip's average momentum squared 
(vs. Ncoii in the left panels and vs. Npart in the right panel). The left panels contain preliminary 
data [236] ■ the right panels final ones [37] . 



dependence of the regeneration yield on charm-quark equilibration than in Ref. |165] : the original pre- 
dictions of Ref. |38] were done with r^^ = 3 fm/c and app = 570 /ib [92] (consistent with the PHENIX 
cross section), which somewhat overestimate the J/ifj yield in central Au-Au as compared to the most 
recent PHENIX measurement [37] (cf. dash-double-dotted line in the upper right panel of Fig. [37|) . 
Using r^^=5-7fm/c [97] improves the agreement, implying a decomposition into primordial and re- 
generated component rather comparable to the one shown in the upper left panels [165] . A potential 
problem for both calculations are the (two) data points at Npart = 100-200, and has been addressed 
in threshold- melting scenarios [237] . HSD calculations at RHIC based on suppression-only scenarios 
(comover or threshold melting) [184] underpredict existing PHENIX data, thus supporting the presence 
of a regeneration component. 

The J/ip regeneration yield is rather sensitive to the open-charm content of the system. First and 
foremost, this pertains to the number of charm- ant icharm pairs, Ncc, in the system. In the grand- 
canonical limit {Nop 3> 1) the charmonium equilibrium number depends quadratically on N^s- But 



even at SPS energies, where charm-quark production is small and the canonical limit applies (implying 
that N^ depends linearly on Ncc), it has been argued that the (semi-) central NA50 J/ip data can be 
accounted for by statistical production alone |162] using a factor of ~2 increase of the pQCD open-charm 
cross section (amounting to dacc/dy = 5.7 fih; a recent data compilation of open-charm data [172] results 
in a value of alf{Eiab = 158AGeV) ~ 3.6-5.2 /ib for the total cross section). It furthermore matters 
over what range in rapidity the statistical model is applied. In the canonical limit, it follows from 
Eqs. (H71) and (HH!) that A^^ oc l/VpBiNpart), which reproduces the observed centrality dependence for 
semi-/central collisions. The width of the considered rapidity window also affects the total yield due to 
different ^-distributions of Nch and Ncc- Finally, recalling Fig. |30l even for fixed Ncc, the charm-quark 
fugacity, and thus the equilibrium charmoniun numbers, are sensitive to the underlying open-charm 
spectrum (e.g., due to medium effects in charm hadrons and/or quarks). 



Finally, we comment on two recent developments. In Ref. |239j the regeneration processes at RHIC 
have been studied by performing Langevin simulations for both c and c quarks in a hydrodynamically 
evolving medium and treating their mutual interaction after production via a potential taken as the 
heavy-quark internal energy provided by IQCD. The combimation of a small spatial c-quark diffusion 
coefficient with a strong c-c attraction (deep potential) keeps the original cc pair, whose primordial 
production essentially occurs at the same space-time point, close together and results in a sizable prob- 
ability to form a charmonium bound state in the subsequent QGP and mixed-phase evolution. This 
"diagonal" recombination of c and c from the same primordial cc pair is not unlike the "canonical limit" 
of the statistical model, augmented by a "correlation volume" to localize it in space. The correlation 
volume is of course replaced by a dynamical calculation using a microscopic diffusion process. The cen- 
trality dependence of the PHENIX data for the inclusive J/ip yield can be approximately reproduced, 
and the pt spectra exhibit a significant softening caused by the softening of the (partial) thermalization 
of the c-quark spectra. 

In Ref. [238j the thermal rate-equation approach (underlying the right panels in Figs. [35] and 1371) 
has been employed to study the consequences of weak- and strong-binding scenarios for charmonia. 
Pertinent charmonium spectral functions have been constructed for these two scenarios based on the T- 
matrix calculations of Ref. |103j for either the free (F) or internal energy (U) as potential, respectively. 
In quasiparticle approximation, the main quantities characterizing the spectral functions and figuring 
into the rate equation are the temperature-dependent charmonium binding energies, charm-quark mass 
and inelastic width (the former two define the charmonium masses). The resulting spectral functions 
were constrained to yield correlator ratios close to one (in line with IQCD). In addition, updated nu- 
clear absorption cross sections have been implemented into the initial conditions, i.e., increasing from 
(Tabs — 4.5 mb to ~7.5mb at SPS energy in line with recent p-A data by NA60 at 158 GeV projectile 
energy [188] . and from 1.5 mb to ~4mb at RHIC energy (which also accounts for shadowing) [259j. 



Some of the findings are as follows: the strong-binding scenario (with a J/ip dissociation temperature 
of Tdiss — l-8Tc) entails a very small "anomalous suppression" of primordial J/ip^s at SPS; most of the 
observed anomalous suppression (which is small due to the large cxabs — 7.5 mb) is due to Xc and ip' sup- 
pression, and the regeneration contribution is small (not unlike the right panels in Fig. [35]) . In central 
Au-Au at RHIC, the partition of (suppressed) primordial and regenerated components is comparable 
(not unlike the right panels in Fig. [H71) . On the other hand, in the weak-binding scenario (with a J/ip 
dissociation temperature of Tdiss — l-lTc), the regeneration component exceeds the suppressed primor- 
dial one already for central Pb-Pb at SPS, and dominates the yield for central Au-Au at RHIC. While 
both scenarios can be compatible with the centrahty dependence of the total yields, i.e., RAA^Nps^^t), the 
different pt dependencies of the two components offer a more promising observable for discriminating 
weak- and strong-binding scenarios. 



3.3.2 Transverse-Momentum Spectra and Elliptic Flow 

Charinonium py-spectra are believed to shed more light on the suppression mechanisms (in particular, its 
momentum dependence, recall Fig. [18]), as well as to disentangle (suppressed) primordial and regenerated 
production components. Since latter is expected to follow an exponential shape (plus effects from 
collective flow), vs. the power-law spectra of primordial production, it should primarily contribute at low 
momentum (where the open-charm phase space density is largest). The pr-dependence of suppression 
may thus be more directly studied toward higher momenta (unless regeneration is absent (or small) 
altogether, as is expected at lower collision energies). At high momenta, additional effects may lead to a 
reduction in charmonium suppression, e.g., escape from the fireball ("leakage effect") ^51 II 77^11 78il240j . 
delayed formation time due to Lorentz time dilation [177^ 1197^ 1241] (coupled with an inhibited formation 
in the QGP zone not unlike a leakage effect, or using reduced dissociation cross sections during the build- 
up of the charmonium wave package |1811I242] ). or the Cronin effect (a broadening of hadron pt spectra 
in p-A relative to p-p collisions) |243l 12441 1245] . 

In lower panels of Fig. [35] the centrality dependence of the J/tp^s average p^ in the hydro+gluo- 
dissociation (left) and fireball-|-quasifree dissociation (right) calculations is compared to NA50 data. 
Most of the observed increase of {p'^) with centrality is accounted for by the Cronin effect. In both 
calculations the latter is implemented into the initial condition for the pt spectra (at the start of the 
thermal evolution, r = tq) using an impact-parameter dependent average nuclear pathlength, {l){b), of 
the parton prior to fusing into charmonium, {j)'^){h; Tq) = {Pt)pp + CigN ■ (O(^) with ag7v=0.76GeV^/fm 
as extracted from p-A data [246] . The reduction of the final (p^) toward central collisions in the lower 
right panel of Fig. [35] [97] is mostly caused by the 3-momentum dependence of the quasifree dissociation 
rate (cf. right panels in Fig. [TSll . leading to larger suppression for higher momentum charmonia (the 
smallness of the regeneration yield, evaluated in blast- wave approximation at Tc, renders its impact on 
(Pt) practically negligible at SPS). A similar reduction is observed in the lower left panel of Fig. [351 [227] . 
despite the decrease of the underlying dissociation rate with 3-momentum. In Ref. [247] . the J/ip and ip' 
spectra in central Pb(158 AGeV)-Pb collisions have been analyzed using a blast-wave description alone; 
the extracted temperature (T ~ 170 MeV) and collective expansion velocity {v ~ 0.2c) are consistent 
with the charmonia forming at the hadronization transition (quite similar to the regeneration component 
in the lower left panel of Fig. [35]) . However, at high p^ > 3 GeV, and toward more peripheral collisions, 
this description most likely needs to be supplemented by non-thermal (primordial) contributions. 

At RHIC energies, the regeneration component in both the hydro+gluo-dissociation and the fire- 
ball-|-quasifree dissociation calculations predicts a softening of (p^) with centrality, relative to the 
Cronin-enhanced initial production, see the lower panels of Fig. [37] Current PHENIX data support 
this scenario, but one should note that the Cronin effect is not yet well constrained from available 
d-A data. The analysis of J/ip p^-spectra in terms of a blast-wave description does not yet allow for a 
quantitative identification of the temperature and transverse flow associated with the kinetic decoupling 
of charmonia. However, the current PHENIX data suggest that a decoupling at the phase boundary 
results in spectra which are too soft [162] . especially for noncentral collisions. This is illustrated more 
explicitly in Fig. [38] where the J/ip Raa{pt) is displayed in 4 centrality bins and compared to the 
fireball -|-quasifree dissociation approach [97]. The rather soft regeneration spectra (corresponding to 
decoupling at Tc with a flow velocity increasing for more central collisions), together with the much 
harder primordial spectra (including suppression), result in a rather fiat Raa{pt) (suggesting that both 
primordial and regeneration components are relevant) Recent measurements of high-pr J/ip^s in Cu- 
Cu(y'sJviv=200GeV) indicate that RAA{PT>5GeV) is compatible with one [249] . At first sight, and in 
view of the strong suppression of other measured hadrons thus far {RaI — 0.25), this appears to be 
a surprise. However, estimates of bottom feed-down [B — t- J/tp + X) and formation-time effects (via 
a reduced absorption cross section at early times), indicate that this observation may be understood 
in the framework of the rate-equation approach [250] . Formation time effects coupled with threshold 
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Figure 38: J/^ nuclear modification factor as a function of transverse momentum in Au- 
Au(y^S7V7v=200GeV) collisions. Left panels: PHENIX data [37J compared to calculations in a fire- 
ball+quasifree dissociation model (including regeneration) [97] in various centrality classes (dotted 
lines: primordial production including Cronin effect and nuclear suppression, dashed line: primordial 
production after QGP and hadronic suppression, dash-double-dotted line: primordial suppressed pro- 
duction without leakage effect, dash-dotted line: regeneration yield in blast- wave approximation, full 
line: total). Right panels: suppression of primordial charmonium states (second from right) and in- 
clusive J 1'^ (far right; including feed-down as well as the sensitivity to the speed of sound employed 
in the hydrodynamical expansion) within a hydrodynamic model employing a QGP threshold-melting 
scenario supplemented by formation time effects [248j . 
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Figure 39: Theoretical predictions for elliptic flow of Jj^l) mesons in 20-40% central Au- 
Au(-^/s/vAr=200GeV) collisions (corresponding to an average impact parameter of 6 ~ 7.8 fm). Left 
panel: transport calculations in the hydro+gluo-dissociation approach |165] : right panel: weighted av- 
erage of coalescence model and suppression calculations using the pr-dependent weights computed in 
the fireball+quasifree dissociation approach [97 t 1250] . 



suppression scenarios |177[ll78j have indeed predicted i?*^ — )■ 1 at high p^ ~20 years ago, see, e.g., right 
panel of Fig. |3H] [248j. A rather unique signature of this effect is a stronger suppression of J j^l) relative 
to Xc; which is opposite to both conventional suppression and statistical hadronization scenarios. 



A promising observable to discriminate suppressed primordial production from regeneration is be- 
lieved to be the azimuthal asymmetry of J /ip production in the transverse plane, quantified by the 
second Fourier coefficient, V2{pt), in the azimuthal angle, 0, relative to the direction of the impact 
parameter vector (commonly identified with the positive x-axis), 

;t^ = 7^ ^ 1 + 2v2ipT) cos 20 + . . . . 57 

d^PT 2'KpTdpT 

Direct and regenerated charmonia are expected to exhibit a large difference in V2- For the former, 
the azimuthal asymmetry is solely due to different absorption path lengths of the J /if) traversing the 
almond-shaped overlap zone in a noncentral nuclear collision: a shorter path length along the "short" 
(x-) axis of the almond implies less suppression, while there is larger suppression along the "long" {y-) 
axis. Calculations typically predict a corresponding V2 of up to 1-3 % |165[ 1251] , see Fig. |39l On 
the contrary, if charmonia are regenerated from c-c coalescence, their V2 is largely determined by the 
underlying charm-quark elliptic flow, fg''^, and expected to approximately "scale" as V2{pt) — 2f2(pT/2) 
(assuming V2=v^ [252| 1164] 1255)^4 Such an approximate scaling has been observed at RHIC for light and 
strange mesons and baryons, with the inferred light-quark V2 reaching up to ~7-8%. If the charm-quark 
f 2 reaches similar values (as suggested by theoretical predictions |253] for semileptonic electron-decay 
spectra |254] ). the elliptic flow of regenerated J/V^'s may reach up to ~15% |252[ I164[ 1255] (see right 
panel of Fig. [39]) . about an order of magnitude larger than in suppressed primordial production. In a 
two-component picture, the combined V2 is given by the weighted sum, v^^ipr) = /^"™(pt) V2^^^{.Pt) + 
f^'^^ipr) v^2^{pt)i where P^^™{pt) + f^^^{PT)=^ describe the fraction of primordial and regenerated 
J/ip's at each value oi pr- The left panel of Fig. [391 displays the results of the hydro+gluo-dissociation 
transport model (employing thermal c-quark distributions) |165j . While the V2 of the regenerated 
component is large, it carries a small weight (decreasing toward higher px), leading to a total V2 of 
around 1%. This will be very difficult to detect experimentally. In the right panel of Fig. |39l the total 
J/ijj V2 (represented by the band) has been estimated by combining blast- wave results for the regenerated 
component |164t 1255] with the suppression calculations of Refs. |165t 1251] , using the weighting following 
from the fireball-|-quasifree dissociation model [97]. The maximum value for the total V2 of ~3±1% is 
somewhat larger than in Ref. |165] . mostly due to the larger V2 in the coalescence component. Most of the 
difference in the f 2 of the regenerated component in Refs. [165] and |250] is presumably due to the fact 
that the former accounts for continuous regeneration throughout the QGP while the latter approximates 
the production at T^. Another difference could be due to the underlying formation reaction, which is 
c + c — 7- J /ip + g in Ref. |165] compared to c + c+p — )■ J /ip+p {p=q, q or g) in Ref. [97]. The main point, 
however, is the overall smallness of the J ftp elliptic flow in both approaches. Note that for semi-central 
collisions [Npart — 150), for which the V2 is evaluated, both calculations underestimate the inclusive 
yield (cf. upper panels of Fig. ET]), leaving room for extra regeneration and thus larger elliptic flow. 

3.3.3 Rapidity Distributions 

The large mass of charm and bottom quarks implies that their rapidity {y) distributions in hadronic 
(and nuclear) collisions are narrower than those of light and strange quarks/hadrons, even at the LHC 
Since the equilibrium number of charmonia is given by N^ oc N'^^/Vfb (grand-canonical limit), and 
the volume of the flreball at a given temperature (e.g., T^) is determined by the light particles, the 
y-distribution of the equilibrium number of charmonia will be narrower than that from primordial 
production, typically given by N^^™ oc N^c- This is to be contrasted with the suppression of primordial 
production, which is expected to become weaker with decreasing light-particle density, thus generating 



^^More precisely, the Q- value of the reaction should be small, i.e., the meson mass should be close to the quark-antiquark 
threshold j255) : this should be a good approximation for loosely bound J/ip^s. 
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Figure 40: Rapidity dependence of J /if) production in Au-Au(ys/vw=200GeV) collisions at RHIC. 
Left panels: calculations for the rapidity density in central (first left) and semi-central (second left) 
Au-Au within the statistical hadronization model |162] . compared to PHENIX data |37j. The cen- 
tral lines (and shaded error bands) correspond to (the uncertainty in) pQCD charm cross sections, 



^pp—^^^-im 



^'^^^ /ib |258] . Right panels: calculations for RAA{.Npart) within the comover- interaction model 
(including regeneration) at mid- (second from right) and forward (far right) rapidity |186] . compared 
to PHENIX dielectron and dimuon data 



a broader \E' rapidity distribution than the primordial one (more suppression at mid- relative to forward 
rapidity). 

The rapidity distributions of the PHENIX J /ip data in Au-Au exhibit a narrowing relative to p- 
p collisions which is quite consistent with regeneration only, as demonstrated by calculations in the 
statistical hadronization model |162] . cf. the two left panels in Fig. HOI Since the J/ip abundance in 
the statistical model is directly proportional to the (squared) open-charm number in the system, the 
open-charm cross section is the key input to compute the J/ip rapidity distribution. Alternatively, 
the rapidity dependence has been evaluated in the comover-interaction model in Ref. |186j (two right 
panels in Fig. HU]) . Including shadowing and comover dissociation on the primordial component, the 
total suppression of the latter is about equal at mid- and forward rapidity (a larger comover density 
around y=0 is essentially compensated by larger shadowing at y=1.2-2.2). Only upon inclusion of a 
gain term does the J/ip production at mid-rapidity slightly exceed the one at forward y. Very similar 
results are obtained in the thermal rate-equation approach [25& j. Recent d-Au data indicate that CNM 
effects may account for a large part of the observed rapidity dependence [257] . 



3.3.4 Excitation Function 

A promising observable to identify the interplay of (suppressed) primordial and secondary charmonium 
suppression is an excitation function. This is due to the much stronger increase of charm- relative to 
light-quark production over a large range in energy. Original predictions of a 2-component model [H] 
for central A-A collisions envisaged a transition from a suppression-dominated regime at SPS to a 
largely regeneration-dominated one at top RHIC energy, cf. left panel of Fig. Hi] While the more recent 
theoretical analyses (recall upper panels of Fig. EZ]) suggest that the current RHIC data contain a smaller 
fraction of regeneration than originally predicted, it is important to realize the following (experimental) 
uncertainty pertaining to secondary production in the commonly used nuclear modification factor. The 
latter is generically defined as the yield (or spectrum) of hadron h in A-A normalized to the binary- 
collision scaled yield in p-p, R\a = ^aa/ i^couNpp) ■ For (suppressed) primordial J/ip production, the 
dependence on the initial production thus drops out (to the extent that primordial production follows 
binary collision scaling). This is not the case for the (absolutely normalized) regeneration component. 
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Figure 41: Collision-energy dependence of J/ip production, normalized to the number of cc pairs, in 
ultrarelativistic heavy-ion collisions. Left panel: original predictions of the 2-component model [HJ [92] 
for central A-A collisions from SPS to RHIC energies (using 2 thermal fireballs). Right panel: predictions 
within the statistical hadronization model for the centrality dependence at SPS, RHIC and LHC (for 
Ay=l) including corona effects for primordial production jl62] : the data point is obtained from PHENIX 
p-p measurements |169[ 1259] . 



which in Raa plots is therefore beset with additional uncertainty due to the input cross section from 
p-p collisions (at RHIC the pertinent experimental uncertainty currently amounts to ~20-25% [259J). 

The right panel of Fig. UT] shows calculations of the statistical hadronization model (augmented by 
primordial production from a corona region) for the centrality dependence of J/ip over open-charm 
numbers at midrapidity for top SPS, RHIC and LHC energies [162] . The open-charm production 
cross section is taken from pQCD (using iVcozz-scaling for A-A) along with assuming a fixed value 
of Nji^/Ncc='\-% in p-p. The increasing trend toward more peripheral collisions at SPS and RHIC 
is largely induced by the canonical suppression factor (for central collisions, the results are in fair 
agreement with the excitation function in the left panel of Fig. W\\ . At LHC, however, the transition 
to the grand-canonical ensemble is essentially completed at A'paj.t— 100, so that for larger centralities a 
distinct increase of Nj/^/Ncc occurs, signaling the dominance of the recombination of uncorrelated c 
and c quarks into J/ip^s. 



3.3.5 Excited Charmonia 

In all phenomenological applications discussed in this section thus far the number of J/ip's refers to 
inclusive production, i.e., contains the feed-down contributions of excited charmonia (~30(10)% from 
Xc{4'') ill p-p reactions). E.g., a comprehensive interpretation of SPS and RHIC data in terms of a 
schematic "sequential melting" scenario has recently been presented in Ref. [260] . In addition to their 
role in J/ip production, excited charmonia provide complementary information on the produced medium 
as their (in-medium) binding energies, dissociation widths and dissolution temperatures (are expected 
to) differ significantly from the J/ip. Consequently, future experimental programs put a large emphasis 
on direct measurements of Xc and tp', as elaborated in Sec. |H E.g., ratios of different charmonia can be 
used to test the occurrence of relative chemical equilibrium. 

NA50 data for the tp' / {J / ip) ratio in Pb-Pb collisions at SPS exhibit a pronounced decrease with 
centrality up to a factor of ~3, see left panel of Fig. HH |261] . It has been suggested |262] that a 



dropping a-meson mass (as a consequence of (partial) chiral symmetry restoration) accelerates the 
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Figure 42: ip' production in ultrarelativistic heavy-ion collisions. Left and middle panel: ip' /{J /ip) ratio 
as computed in the thermal rate-equation approach with quasifree dissociation in the QGP and hadronic 
dissociation with (sohd lines) and without (dashed hues) in-medium D-meson masses at SPS (left; data 
from Ref. [261] ) [38] and RHIC (middle). Right panel: ip' suppression (nuclear modification) factor 
in the HSD transport approach in the comover-interaction (circles) and QGP threshold-suppression 
(triangles) scenarios [263] ; upon dividing the ip' suppression factor by the one for the J /ip (see Fig. 9 in 



Ref. [263] ). the comover result is rather compatible with the rate-equation curve in the middle panel. 



transition rate for ip' -H- J /%p + a close to T^ and can account for the ip' /{J /ip) data. For semi-/central 
Pb-Pb collisions the data can also be explained within the statistical hadronization model [H], although 
here the underlying mechanism is different, i.e., a QGP converting into a hadron gas in thermal and 
(relative) chemical equilibrium (with conserved charm- and anticharm-quark number). In the thermal 
rate-equation approach of Ref. [38], regeneration at SPS energies is a small effect, and the QGP phase is 
too short to provide sufficient ip' suppression to describe the NA50 data for the ip' /{J /ip) ratio (dashed 
line in the left panel of Fig. H2l) . However, if effects of (partial) chiral restoration are implemented 
via an in-medium reduction of D-mesons masses (associated with their light-quark content), the direct 
%p' — )■ DD channel opens and increases the ip' dissociation rate in the hadronic phase close to T^ (in the 
vacuum, the ip' mass is only ~50MeV below the DD threshold). This allows for a good description 
of the NA50 ip' data (solid line in the left panel of Fig. H2|) . Alternatively, a broadening of the in- 
medium D-meson spectral function can induce a similar effect as reduced masses, since it also opens 
phase space for direct ip' — )■ DD decays. Either mechanism is not mutually exclusive to a relative 
equilibration via J /ip + a -v^ ip' (neither is statistical hadronization, which, however, requires the 
prevalence of regeneration). At full RHIC energy, both the thermal rate-equation approach (middle 
panel in Fig. H2|) and HSD transport calculations using pre-hadronic and comover interactions with 
regeneration (right panel in Fig. H2]) predict a ip' /{J/ip) ratio which is very similar to SPS energies. On 
the other hand, threshold-melting scenarios predict significantly sharper suppression patterns in the 
centrality dependence. The Xc/{J I'lp) ratio is expected to drop with centrality in both rate-equation 
and comover approaches, and to decrease more sharply for threshold-melting scenarios. 



3.3.6 Bottomonia 

Bottomonium observables add at least two aspects, relative to charmonia, to the study of the medium 
in heavy-ion collisions at RHIC and LHC: (i) larger binding energies (parametrically, for Coulombic 
wave functions, e^ oc a^rriQ), implying larger dissociation temperatures and thus different dissociation 
patterns; (ii) a larger mass of the bottom quark, reducing its primordial production (e.g., N^^i ^ 0.1 per 
central Au-Au at RHIC [264J), which suggests that recombination effects are suppressed. The latter 
assertion, however, has to be taken with care since primordial T production is also relatively small 
(typically N-c/Nf^i ~ 0.1 % [264[ 1265] compared to 1 % for charm/onium) , and bottom-quark fugacities 
can be rather large. 

Early studies of bottomonium production in heavy-ion collisions have focused on suppression ef- 
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Figure 43: Predictions for T production at RHIC within the thermal rate-equation approach using 
quasifree dissociation with in-medium binding [21] in an expanding fireball. Left panel: time evolution 
of T abundance in central Au-Au; right panel: centrality dependence of inclusive T yield. 



fects [1541 1180"] . Similarly, in Ref. [266], the threshold- melting scenario, combined with formation-time 
effects, has been utilized to compute the pt dependence of suppression factors, S{pt) = Raa{pt), at 
the LHC Assuming initial temperatures of Tq ~ 0.8-1.2 GeV, a large suppression for all states has been 
predicted at low pt, which, however, gradually disappears in the range of pt — 5-30 GeV (qualitatively 
similar to the right panel in Fig. 155]) . depending on the Debye mass and surface effects in the fireball. 
An extension of the threshold-melting scenario to RHIC [248] suggests that directly produced T's do 



not undergo any suppression, while inclusive T's (which receive about 40-50% feed-down from Xb and 
T') are suppressed by up to 45% at low pt due to an almost complete melting of the excited states. 
As mentioned above, the threshold-melting picture represents a sharp transition between stable and 
dissolved quarkonia, i.e., a sudden "jump" of the inelastic width from zero to infinity. This is a sim- 
plification of a more realistic situation where the width increases with temperature more gradually, as 
found in microscopic calculations discussed in Sec. 12.31 (in the context of Fig. [19]). On the other hand, a 
rather abrupt dissolution of a quarkonium bound-state as a function of temperature may be justified by 
potential-model calculations of spectral functions, recall Sec. 12.21 The reason is that, once the binding 
approaches zero, the QQ state in the QGP moves into the continuum and the structureless potential can 
no longer support resonance states. In addition, a rapid increase of the dissociation width is expected 
once the direct decay channel opens, T— )'6 + 6or\l/— j-c + c, i.e., above the quark-antiquark threshold 
in the QGR 

Quantitative studies of bottomonium suppression and regeneration have been performed in Ref. [99] 
within the thermal rate equation fl23]) . solved in an expanding fireball background. When assuming 
vacuum binding energies for the T throughout the QGP evolution at RHIC (with an initial temperature 
of To=370MeV, dropping to 300 MeV within the first 0.5 fm/c), the results of Ref. [248] have been 



confirmed in that no significant suppression of direct T's (i.e. those not from feed-down) occurs (nor 
regeneration). This is understandable due to its long lifetime of several tens of fm/c for temperatures 
below 300 MeV (recall left panel in Fig. [T9]) . The situation changes appreciably when in- medium 
reductions of the bottomonium binding energies are included (based on the in-medium rates displayed 
in the left panel of Fig. Hn]). In central Au-Au at RHIC, direct T's are suppressed by ~ 40%, with most 
of the suppression occurring in the first Ifni/c, i.e., for temperatures above 250 MeV (cf. left panel of 
Fig. US]). In the centrality dependence of inclusive T production (right panel in Fig. |43]) the suppression 
is further amphfied, due to ~50% feed-down contributions [267] from excited Xb (~30% Xbi, ~10% Xb2) 
and T' (~10%) states; regeneration contributions are (very) small. An important point here is that 
direct T suppression at RHIC is a rather direct measure of the Debye screening in the heavy-quark 
potential, especially since regeneration is absent. The binding energies figuring into the calculations of 
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Figure 44: Predictions for the centrality dependence of inclusive T production at LHC. Left panel: 
thermal rate-equation approach [99] displaying the final yield, A^t, within a thermal fireball of rapidity 
width Ay = 1.8 (the total fireball entropy of S" = 20000 in central Pb-Pb converts into dNch/dy ~ 1600). 
Right panel: statistical hadronization model |162j (based on dNch/dy ^ 1800); the yellow uncertainty 
band reflects a variation of the input pQCD bb cross section by a factor of 1.5 up and down. 



Ref. [99] are somewhat smaller than potential models based on the internal energy (but significantly 
larger than those based on the free energy), as discussed in Sec. 12.21 In either case, one could envisage 
a situation where the production of T is as (or even more) suppressed as for J/ip. This would be a 
rather unique QGP signature [99j, and provide indisputable evidence of charmonium regeneration and 
color-Debye screening. 

At the LHC, the situation becomes more involved. The rate-equation approach (left panel in 
Fig. mi) [99j predicts a suppression of up to a factor of ~10 stronger than the statistical hadronization 
model (right panel in Fig. HI]) |162j in central Pb-Pb collisions. Besides moderate differences in input 
cross sections (e.g., a factor of 0.8 shadowing correction in Ref. [99]), there are several sources of this 
difference: (a) the equilibrium limit in the rate equation is computed in the QGP with (in-medium) 
bottom-quark masses (rather than bottom hadrons), leading to smaller bottom-quark fugacities; the 
resulting equilibrium number is up to a factor of 3-6 smaller, depending on the (not very well known) 
spectrum of bottom hadrons. (b) the equilibrium limit in the rate equation is reduced by a schematic 
thermal relaxation-time factor, TZ{t) < 1, to mimic incomplete 6-quark thermalization; at the end of 
the QGP phase this factor has reached to about 0.3-0.5 (it is smaller earlier in the evolution). However, 
close to Tc, the inelastic reaction rates for the bottomonia are already rather small; (c) as mentioned 
for charmonia at SPS, there is a question over which range in rapidity the thermal ensemble of bottom 
quarks should be defined (e.g., Ay=1.8 and 1 in Refs. [99] and |162j . respectively). More work is needed 
to clarify these questions. 



4 Future experiments 



This chapter is devoted to the presentation of future opportunities of quarkonia measurements in heavy- 
ion collisions. Detectors are currently being designed, built, installed, commissioned or upgraded around 
three accelerators which will deliver heavy-ion beams in three different energy domains: 



• RHIC-II (starting time scheduled for 2013) is an upgrade program of the existing RHIC (Rela- 
tivistic Heavy Ion Colhder) at Brookhaven National Laboratory, USA, which delivers beams up 
to ^/SNN = 200 GeV for Au. The upgrade of the machine consists of an increase of the present lu- 
minosity by an order of magnitude as well as a low-energy program {^/sNN = 5-15 GeV). Whereas 
the low-energy physics program will mainly focus on evidence for the critical point, the high- 
luminosity high-energy physics program aims at completing existing measurements with a strong 
emphasis on the heavy-flavor sector. To achieve this goal, an upgrade program of two detectors 
(PHENIX and STAR) is underway; 

• The LHC (Large Hadron Collider) at CERN will provide, from 2009 on, the largest energy ever 
delivered by an accelerator with proton beams up to ^/s = 14 TeV and Pb beams up to ^/SNN = 
5.5 TeV. Three of the four LHC experiments (ALICE, ATLAS and CMS) will take heavy-ion data 
with an important part of their physics program devoted to quarkonia measurements. Whereas 
ATLAS and CMS are designed for proton-proton physics, ALICE is dedicated to the study of 
heavy- ion collisions; 

• The SIS-300, located at the FAIR (Facility for Antiproton and Ion Research) at GSI, Darmstadt, 
will deliver, from 2015 on, heavy ion beams in fixed target mode with beam energy from 10 AGeV 
to 45 AGeV (^s^ = 4.5-9.3 GeV). The Compressed Baryonic Matter (CBM) detector will al- 
low to explore, with high-statistics measurements of rare signals like quarkonia, the moderate 
temperature and large baryon chemical potential region of the QCD phase diagram. 

In two of these energy ranges (LHC and SIS-300) heavy quark resonances will be measured for the 
first time. This ensures a rich and exciting future physics program. In the following, we present 
the motivations for studying quarkonium production in these three energy ranges as well as a short 
overview of the detector capabilities needed to achieve this goal. All detectors have in common the 
goal of performing quarkonium measurements as a function of centrality and reaction plane, in a broad 
acceptance, with large statistics and low background, to ensure a mass resolution good enough to 
separate all vector resonances and to identify feed-down from excited states. For more details we refer 
to [26H], to [2691 |2IQ1 [2711 M and to [273] for the RHIC-II, the LHC and the SIS-300 experimental 
programs, respectively. 

4.1 Quarkonia in heavy-ion collisions at the RHIC-II 

The RHIC has started delivering heavy-ion beams in 2000. Since then, the four heavy-ion experiments 
(BRAHMS, PHENIX, PHOBOS and STAR) have collected a large set of data for various systems at 
energies a/sjvtv = 20 — 200 GeV. The status of information gathered up to 2004 has been summarized 
both from the theoretical side |274] and from the experimental side |275] . Hard probe precision mea- 
surements have become accessible at RHIC with the high statistics runs performed since 2004. The 
measurements, done by the PHENIX and STAR collaborations, have revealed some of the most exciting 
aspects of the hot and dense matter produced at RHIC. A particularly striking and unexpected obser- 
vation in the quarkonium sector is the J/ip suppression whose magnitude in mid-central and central 
heavy- ion collisions is similar to that measured at SPS energy for the Pb-Pb system (recall Figs. 1531 and 
[57|) . As a consequence, several models which offer a successful quantitative description of J/ip suppres- 
sion at SPS energy (assuming that the suppression is driven by the system temperature or its energy 
density) over-predict the suppression observed at RHIC. Furthermore, the observed J/ip suppression is 
larger in the forward region than in the central region (Fig. HOj) . This observation is again striking as 
most of the suppression models predict an opposite effect, i.e., the higher the energy density the larger 
the suppression. The interpretation of these observations has been discussed in Sees. 13.21 and 13.31 



Further detailed measurements are mandatory in order to disentangle the different theoretical inter- 
pretations and to achieve a comprehensive understanding of quarkonia production in heavy-ion collisions 
at RHIC Indeed, direct evidence for resonance dissociation by color screening requires the centrality 
dependence of J /'ip, Xc and ip' to be measured simultaneously. Regeneration models can be constrained 
from the inspection of J /ip rapidity distributions which are expected to narrow |162[I168] (recall Fig. HOl) . 
as well as from a precise measurement of the (open) charm production cross-section. Trends in the exci- 
tation function of J /ip suppression could also reveal the relative importance of the regeneration process 
as it should increase with the multiplicity of cc pairs [H] . More stringent tests include measurements of 
J /iIj elliptic flow since, as open charm shows non-zero flo"fto, J /i^^ if produced by regeneration, should 
inherit this flow (recall Fig. I39l) . In addition, both J /ip flow and suppression should exhibit a specific 
Pt dependence marked by QGP effects at relatively low pt whereas the high px region should be mostly 
populated by J /ip which could escape the medium unaffected. In this respect it is also important to be 
able to identify secondary J /ip from bottom decay whose yield might be sizable at high pt- Further- 
more, detailed measurements in p{-Y\ke)-A collisions are mandatory in order to disentangle cold and hot 
nuclear effects. We finally note the possibility to identify the presence of the deconfined medium from 
the measurement of J/ip polarization |276] . In addition to the above observables based on charmonia, 
the measurement of bottomonium states should provide additional insights. The temperature of the 
medium at RHIC is expected to be large enough to induce the break-up of T{2S) and T(3S') whereas 
the T(IS') is supposed to melt at higher temperatures (it might still be affected by inelastic reactions 
with partons in the QGP, as discussed in Sees. 12.31 and [3. 3. 6p . 

The previously exposed physics program can partially be covered in the next years with the existing 
RHIC machine and detectors. However, it is clear that the full physics program requires high statistics 
and excellent data quality. This quality will be achieved by means of an upgrade of the machine 
luminosity and of the detectors. 

4.1.1 Experimental conditions 

With the luminosity upgrade of the machine, the integrated luminosity per week is expected to be 
2500 /ib~^ and 33pb~^ for Au-Au and p-p collisions at ^/sNN = 200 GeV, respectively. A typical run at 
RHIC consists of 24 weeks of data taking. This is significantly longer than a LHC run in the heavy-ion 
mode (see below). As a consequence, although the heavy flavor production cross sections are much 
bigger at the LHC, the lower cross sections at RHIC-II are compensated by the integrated luminosity 
so that the heavy-flavor yields for one year of running are expected to be similar at RHIC-II and at the 
LHC. 



4.1.2 The PHENIX experiment 

The PHENIX experiment is made of four spectrometers and a set of specialized detectors to determine 
the collision centrality and to provide triggers (Fig. |35]left). In the central region, two arms allow 
to detect electrons, photons and charged hadrons. They consist of a complex arrangement of drift 
chambers, pad chambers. Time Expansion Chambers (TEC) operated as transition detectors, time-of- 
flight (ToF) detectors, RICH detectors, aerogel detectors and electromagnetic calorimeters (EMCs). A 
magnet supplying a fleld parallel to the beam is placed around the interaction vertex. In the forward and 
backward regions, two spectrometers allow muon measurements in the rapidity range —2.25 < y < — L15 
and 1.15 < y < 2.44. They consist of a muon tracker (three stations of multi-plane drift chambers) 



^•^In contrast to light-flavor hadrons, the heavy-flavor flow is, so far, not measured experimentally through identified 
hadrons, but in an inclusive way via the flow of non-photonic electrons. The latter is obtained from full distributions of 
electrons after subtraction of Dalitz-decay electrons from light hadrons and photon-conversion electrons. 
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Figure 45: Longitudinal view of the PHENIX (left) and STAR (right) detectors. 



placed inside a radial magnetic field, followed by a muon identifier (alternating layers of steel absorbers 
and streamer tubes for tracking) both with full azimuthal coverage. 

The upgrade program of PHENIX [268] consists of a barrel and two endcaps Silicon Vertex Detectors 
(SVDs), two Nose Cone Calorimeters (NCCs) and an upgrade of the muon trigger. The SVD will 
provide inner tracking with full azimuthal coverage for |?7| < 2.5. This will improve the dilepton 
mass resolution and reduce the background both in the electron and in the muon channels. It will 
further allow to measure displaced vertices which is mandatory for identification of secondary J/ip^s 
from 5-meson decay. The NCCs contain electromagnetic and hadronic sections covering the acceptance 
0.9 < |?7| < 3.5. These are key detectors for measuring Xc -^ J/ip + j decays. The muon trigger upgrade 
consists of adding three layers of RPC detectors in each muon arm together with associated front-end 
electronics and trigger logic. This will improve the level-1 trigger selection for both single and di-muons. 

Typical expected yields with the upgraded PHENIX detector for a 12-week physics run with Au 
beams at RHIC-II are 45000 (395000) J/^ and 400 (1040) T in the electron (muon) channel. 



4.1.3 The STAR experiment 

The different sub-systems of the STAR experiment (Fig. HS] right) are placed in a solenoidal magnet 
operating at 0.5 T. The main detectors are (i) a large azimuthally symmetric Time Projection Chamber 
(TPC) providing charged particle tracking within |?7| < 1.8, (ii) the Silicon Vertex Tracker (SVT) and 
Silicon Strip Detector (SSD) for charge particle tracking close to the interaction point and vertexing, 
(iii) the Forward TPCs (FTPCs) covering 2.8 < |?7| < 3.8 to track particles at forward and backward 
rapidities, (iv) the ToF system providing additional particle identification in the acceptance of the 
central detectors and (v) the EMC (both in the central region and at forward angles) which ensures 
neutral particle measurements. Additional specialized detectors complement the setup for centrality 
measurements and trigger. 

The upgrade program of STAR |268j includes a ToF system covering the full acceptance of the central 
barrel, new front-end electronics for the existing TPC, an upgrade of the data acquisition system and 
a new inner tracking system based on two layers of silicon pixel sensors and three layers of silicon strip 
detectors. The ToF will extend the momentum range for hadron identification and, in conjunction 
with the electromagnetic calorimeter, provide a level-2 trigger for J/ijj — )■ e~^e~ measurements. The 



new inner tracking system will make displaced vertex measurements accessible. As for the PHENIX 
experiment, this will be used to tag J /ip from B hadron decays, to reconstruct charm hadron in their 
hadronic decay channels and will allow open heavy flavor measurements from single electrons. Finally, 
the DAQ upgrade will allow to collect rare signals with large statistics. 

Simulations indicate that up to 220 ■ 10'^ J /ip and 11200 T could be reconstructed in STAR during 
a 12 week physics run at RHIC-II in Au-Au collisions. 

4.2 Quarkonia in heavy ion collisions at the LHC 

With a nucleus-nucleus center-of-mass energy nearly 30 times larger than the one reached at RHIC, the 
LHC will provide the biggest step in energy in the history of heavy-ion collisions and open a new era 
for studying the properties of strongly interacting matter under extreme thermodynamical conditions. 
This new energy regime will lead to a much higher energy density, to faster equilibration and to a longer 
lifetime of the deconfined system, resulting in an enhanced role of the QGP over final-state hadronic 
interactions |277j . The high temperature and close to vanishing baryon chemical potential of the system 
will render it close to the conditions of the primordial universe. In addition, heavy-ion collisions at the 
LHC access unprecedented small Bjorken-x values where low-momentum gluons are expected to be 
close to saturation and lead to a significant shadowing effect. As a consequence high-density parton 
distributions could be responsible for a large fraction of particle production |278j . Another exciting 
aspect of this new energy regime is the massive production rate of hard processes. They will provide 
an ideal tool for a detailed characterization of the deconfined medium. The heavy-flavor sector at LHC 
energy is also subject to significant differences with respect to SPS and RHIC energies. First, the large 
production rate offers the possibility to use new and a large variety of observableo. The magnitude 
of most of the in-medium effects is therefore anticipated to be dramatically enhanced. Some of these 
aspects are discussed in the following. 

New observables. The T(IS') state is expected to dissolve significantly above the critical tempera- 
ture, at ~3-4 Tc ^621 EU EH] , which presumably comes into reach only with the LHC. The spectroscopy of 
the T family at LHC should then reveal unique characteristics of the QGP. In particular, the pt depen- 
dence of the T(2S')/T(1S') ratio presents significant sensitivity to the dissociation temperatures |266] 
(see Refs. |99l 1280] for updates). Measuring the T(25') is also particularly interesting in order to un- 
ravel J /ip suppression versus regeneration: the T(25') and J /ip dissociation temperatures are predicted 
to be similar |62] whereas, in contrast to charmonia, bottomonia are expected to be little affected by 
regeneration processes [99]. Electroweak Z° and W^ bosons will be available with large statistics and 
can serve as reference processes for quarkonium suppression studies. 

Large resonance dissociation rate. In addition to comoving hadrons and color screening, quarkonia 
can be destroyed by parton ionization |281] . This mechanism, induced by thermal partons in the QGP, 
starts being effective for temperatures above the critical temperature but below the temperature of 
resonance dissociation by color screening (cf. Sec. 12. 3p . Recall, however, the interplay between the two 
mechanisms on the final quarkonium yields: a large dissociation by color screening implies a low binding 
which facilitates a large break-up rate by parton ionization. Recent estimates |279j (see Ref. [91] for 
an update) of the quarkonium dissociation cross sections show that none of directly produced the J /ip 
survives the deconfined phase at LHC, and that about 20% of the T are destroyed, possibly more if 
color-screening is strong, cf. Fig. HD [99J . 



^''^From RHIC to LHC, it is expected that the number of cc [hh) pairs produced in central nucleus-nucleus collisions will 
increase by a factor 10 (100) |270j : up to 115 cc (5 hb) pairs are expected to be produced in 5 % central Pb-Pb collisions at 
5.5 TeV. This estimate is based on next-to-leading order pQCD calculations and includes shadowing; see Refs. |270[ 1279] 
for more details. 



Large secondary charmonium production. Secondary charmonia yields can arise due to statisti- 
cal hadronization [282J and/or kinetic recombination ^2]. These processes result in an increase of J /ip 
yield with collision centrality roughly proportional to N"^^ (grand-canonical limit). Again, this effect is 
expected to entail dramatic consequences at the LHC. Due to the large number of cc pairs produced in 
heavy-ion collisions at LHC, models predict a qualitatively different centrality dependence of the J /%p 
yield [2791 [283], recall right panel of Fig. gH 

Large charmonium rate from B hadron decay. Another source of non-direct charmonia arises 
from the decay of B mesons. The ratio N{B — )■ J /i(j)/N {direct J/ip) can be determined as follows. 
The number of directly produced J/ip in central (5%) Pb-Pb collisions at 5.5 TeV is 0.5 [279^^ The 
corresponding number of bb pairs (with shadowing) amounts to 4.6 [279] . The B — )■ J/ipX branching 
ratio is 1.16 ± 0.10% [M]. Therefore N{B -^ J/^)/iV(direct J/ip) = 20% in 47r0. The secondary 
J/ip from 5-meson decay must be subtracted from the primary J/ip yield prior to J/tp suppression 
studies. They can further be used in order to measure the 5-meson production cross section in p-p 
collisions |226] , to estimate shadowing in p-A collisions and to probe the medium-induced 6-quark energy 
loss in A-A collisions. Indeed, it has been shown |285j that the pt and rj distributions of those J/ip 
exhibit pronounced sensitivity to 6-quark energy loss. 

4.2.1 Experimental conditions 

The LHC will be operated seven months per year in the p-p mode and one month per year in the 
heavy-ion mode. The corresponding estimated effective running time is 10^ s and 10^ s for p-p collisions 
and A-A collisions, respectively. The expected luminosity for Pb-Pb collisions is about 5 ■ 10^^ cm~^s~^ 
which results in a minimum-bias interaction rate of 4 kHz. As described in Ref. [269], the heavy- (and 
light-) ion runs include, over the first five years of operation, one Pb-Pb run at low luminosity, two 
Pb-Pb runs at high luminosity, one p-A run and one light ion-ion run. In the following years different 
options are considered depending on the first results. 

4.2.2 The ALICE experiment 

ALICE (A Large Ion Collider Experiment) is the only LHC experiment dedicated to the study of nucleus- 
nucleus collisions [2691 1270] . The ALICE physics program also includes the study of p-p collisions which 
will provide the reference for heavy-ion data. The ALICE experiment is designed to perform high- 
precision measurements of numerous observables based on hadrons, leptons and photons, in a broad 
acceptance. 

The detector (Fig. l46l) consists of a central part, a forward muon spectrometer and forward/backward 
small acceptance detectors. The central part of ALICE consists of four layers of detectors placed in 
the solenoidal field {B < 0.5 T) provided by the L3 magnet previously used at the Large Electron 
Positron collider (LEP). From inside out, these detectors are (i) the Inner Tracker System (ITS) con- 
sisting of six layers of silicon detectors, (ii) the TPC, (iii) the Transition Radiation Detector (TRD) 
and (iv) the ToF system based on multi-gap resistive plate chambers. They provide charged particle 
reconstruction and identification in the pseudo-rapidity range \ri\ < 0.9, with full azimuthal coverage 
and a broad px acceptance. The ALICE central barrel will later be equipped with a large acceptance 
{\ri\ < 1.4, A$ = 110°) EMC (not shown in Fig. H6l) . These large area devices are complemented by two 



^^Including shadowing and feed-down from higher states. 

^^This ratio is subject to large uncertainties since (i) the predictions from the color-evaporation model (CEM) on 
J/ip cross sections differ by up to a factor of 2 at LHC energies |279j . and (ii) the total bb production cross section in 
p-p collisions at LHC is predicted within a factor of 2-3 uncertainty resulting from the choice of the quark mass, the 
renormalization and factorization scales and the parton distribution function. In addition, resonance suppression and/or 
regeneration, heavy-quark energy loss and other effects not taken into account here could play a significant role. 
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Figure 46: Longitudinal view of the ALICE detector. 



smaller acceptance detectors: the High Momentum Particle IDentification (HMPID) and the PHOton 
Spectrometer (PHOS). In the forward/backward region, additional detectors (TO, VO and FMD, not 
shown in Fig. \^\i enable a fast characterization and selection of the events, as well as charged-particle 
measurements in the pseudo-rapidity range —3.4 < rj < 5.1. At large rapidities, photon multiplicity 
and spectator nucleons in heavy-ion collisions will be measured by the Photon Multiplicity Detector 
(PMD) and the Zero-Degree Calorimeters (ZDC; not shown in Fig. H6|) . respectively. A forward muon 
spectrometer covering the pseudo-rapidity range —4 < rj < —2.5 complements the central part. It 
consists of a front absorber, a dipole magnet, ten high-granularity tracking chambers, a muon filter and 
four large area trigger chambers. 

The performance of the ALICE detector in the quarkonium sector is summarized in Fig. HTl Quarko- 
nium states will be measured in both electron (|?7| < 0.9) and muon channels (—4 < t] < —2.5). The 
acceptance will allow reconstruction of differential distributions down to very low transverse momentum 
in most cases. The resolution of the apparatus, better than 100 MeV/c^ for invariant masses around 
10 GeV/c^, allows for a separation of all quarkonium vector resonance states. The expected number of 
reconstructed J/i/j and T during one month of data taking is 677 ■ 10^ and 9600, respectively, in the 
muon channel for minimum-bias Pb-Pb collisions. The corresponding numbers in the electron channel 
for central (10%) collisions are 121 ■ 10^ and 1800. 



4.2.3 The CMS experiment 

CMS (Compact Muon Solenoid) is a general-purpose detector (Fig. HH] left) designed to measure muons, 
electrons, photons and jets. Although the detector is optimized for p-p collisions, a strong heavy-ion 
program has been developed |271j . CMS is composed, from the interaction point to the outer side, of 
a tracking system, an EMC, a hadronic calorimeter and muon chambers arranged in a central barrel 
and two endcaps. The central element of CMS is a 13 m long, 3 m diameter magnet which delivers 
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Figure 47: Invariant- mass spectra of low-mass dimuons (left), high-mass dimuons (middle), and dielec- 
trons (right), expected to be measured with the ALICE detector within one month of Pb beam running 
(extracted from Ref. [270] ). In the left and middle plots, the non-correlated background is subtracted 
from the total spectrum assuming a perfect subtraction, i.e., the statistical error of the "full" spectrum 
is assigned to the remaining spectrum of the sum of the correlated sources. 
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Figure 48: Longitudinal view of the CMS (left) and of the ATLAS (right) detectors. 



a i? = 4 T solenoidal field surrounding the tracking and calorimetric systems. The tracker is based 
on several layers of silicon pixel and strip counters and covers the pseudo-rapidity region \ri\ < 2.5. 
The EMC, made of lead-tungstenate crystals, covers the pseudo-rapidity region \ri\ < 1.5 in the central 
barrel. This coverage is extended to \t]\ < 3 with the endcaps. The hadronic calorimeter is made of 
copper plates and plastic scintillator sandwiches. Its acceptance is jr^l < 2 in the central barrel and 
reaches \ri\ < 5.3 with the endcaps. Two additional very forward calorimeters ensure coverage in the 
pseudo-rapidity range 3 < jr^l < 5. The muon system is located outside the central magnet. It consists 
of four layers of detectors (three for tracking and one for trigger) covering the pseudo-rapidity range 
\r]\ < 2.4 {\ri\ < 1.5 in the barrel). Very forward calorimeters, including two ZDCs (3 < |?7| < 5.2) 
and a quartz fiber calorimeter (5.3 < \ri\ < 6.7) allow measurements of the collision centrality and 
electromagnetic energy (these detectors are not shown in Fig. US]). 

Figure HH shows the dimuon invariant- mass spectra expected to be measured for charmonia (left) 
and bottomonia (middle). The muon acceptance ranges in rj from —0.8 to +0.8 in the central barrel 
and from ±0.8 to ±2.4 in the endcaps. The excellent mass resolution of the apparatus will allow a very 
clean separation of quarkonia from the J/V' to the T(3S'). In the acceptance of the central barrel, J/ip 
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Figure 49: Invariant- mass spectra of opposite-sign and like-sign muon pairs in the J/ip (left) and T 
(middle) mass regions expected to be measured with the central barrel of the CMS detector within one 
month of Pb-Pb collisions |271] . Right: corresponding invariant-mass spectrum of high-mass opposite- 
sign muon pairs reconstructed with the ATLAS detector |272] . 



measurements will be limited to pr > 4 GeV/c due to a high pt threshold induced by the calorimeters 
on single muons. This can be extended to lower pt when detecting muons in the endcaps thanks to the 
extra longitudinal Lorentz boost. Due to its higher mass, bottomonium will be measured down to zero 
Pt in the whole r] acceptance. Simulation results indicate that the typical number of reconstructed J/ip 
and T will be 184 ■ 10^ and 37.7 ■ 10^ in minimum-bias Pb-Pb collisions for one month of data taking. 



4.2.4 The ATLAS experiment 

Like CMS, ATLAS (A Toroidal LHC ApparatuS) is designed for p-p physics. The detector capabilities 
for heavy ion physics have been recently investigated }272] . The design of the detector (Fig. HH] right) 
is similar to that of CMS with a tracking system, an EMC, a hadronic calorimeter and muon chambers 
placed in a barrel and two endcaps. The tracking system of ATLAS is composed of silicon pixel detectors, 
the Semiconductor Tracker (SCT) made of silicon strip detectors, and the Transition Radiation Tracker 
(TRT). It is placed inside a -B = 2T solenoidal magnet and covers the pseudo-rapidity range \ri\ < 2.5. 
The electromagnetic calorimeter is a liquid-argon device covering \r]\ < 4.9. It is surrounded by the 
hadronic calorimeter which consists of lead scintillators in the barrel (|?7| < 1.7) and liquid argon in the 
endcaps (1.5 < \ri\ < 3.2). Two additional (electromagnetic and hadronic) calorimeters cover the very 
forward region (3.2 < |?7| < 4.9). The muon spectrometer consists of a toroidal magnet providing a 
B = AT field and several muon chambers with different technologies in the barrel and in the endcaps. 
The acceptance of the muon spectrometer covers \ri\ < 1.0 in the barrel and extends to rj = ±2.5 in the 
endcaps. 

As an illustration of the ATLAS performance for measuring quarkonium signals, the right panel 
of Fig. Sni shows the invariant-mass spectrum of high-mass muon pairs. In the full rj acceptance, the 
mass resolution is 145MeV/c^ at around lOGeV/c^. The expected statistics for T for one month of 
Pb-Pb minimum-bias collisions is 10-15-10^ (with a mnon-pT threshold of 3GeV/c). The number of 
reconstructed J/ijj ranges from 8000 to 216 ■ 10^ depending on the trigger threshold. 



4.3 Quarkonia in heavy-ion collisions at the SIS-300 

While the future heavy-ion experiments at RHIC-II and at the LHC will focus on the study of the QCD 
phase diagram at large temperatures and small chemical potential, the SIS-300 accelerator at FAIR aims 
at exploring the region of moderate temperatures and large baryon chemical potentials. In this region 



of the phase diagram, IQCD calculations predict a critical endpoint whose location is not precisely 
known. Beyond this critical endpoint, for higher baryon chemical potentials and lower temperatures, 
one expects a first order phase transition from hadronic to partonic matter. This regime has been very 
little explored experimentally so far. The corresponding beam energy range (from ~10-45AGeV) has 
only been partially covered in the past by pioneering experiments at the AGS in Brookhaven and, more 
recently, by experiments at CERN-SPS. A non-monotonous behavior was observed in the excitation 
function of several observables, such as the {K~^) / (n^) ratio, at a beam energy of around 30 AGeV [286] . 
An interpretation of this intriguing finding as a signal for the onset of deconfinement is still under debate. 
A large low-mass dielectron enhancement was observed at 40 AGeV by the CERES experiment |287j . 
corroborating the trend in the {K^) / {-k^) excitation function. In fact, experimental limitations, both 
from the accelerators (beam intensities) and from the detectors (capabilities to measure rare signals) 
prevented comprehensive measurements to be performed at such low beam energies. Consequently, 
charm hadrons (and charmonia) have never been measured in this energy range. Such measurements 
could be especially interesting as this beam-energy range is close to the production threshold for cc pairs. 
Therefore, the characteristics of the produced charm hadrons are expected to be particularly sensitive 
to the properties of the medium in the early stage of the collision. We further note that this region of 
the phase diagram is well suited for investigating properties of hadrons in a dense baryonic medium. 
This is based on the increasing gap between the putative phase boundary and the chemical freeze-out 
line when decreasing the beam energy from RHIC and top-SPS to FAIR. It furthermore suggests that 
at FAIR energies, contrary to RHIC and top-SPS, the produced system spends a relatively long time 
in a dense baryonic phasi 
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4.3.1 Experimental conditions 

The SIS-300 synchrotron at the FAIR facility will deliver heavy-ion beams from 10-35 AGeV for uranium 
and up to 45 AGeV for lighter ions with Z/A = 0.5. Proton beams will be available up to 90 GeV. 
Thanks to the unprecedented high beam intensities of 2 ■ 10^/s, rare probes such as J/ip, whose typical 
production rate is 2 ■ 10~^ per central Au-Au collision at 25 AGeV, will become available. Detection 
of such a small signal in such a high-intensity environment, and among about 1000 charged tracks per 
event, is highly challenging. It requires excellent detector performance in terms of radiation hardness, 
read-out speed, online-event selection, particle identification and data processing. The CBM detector 
is currently being designed to meet these requirements. 

4.3.2 The CBM experiment 

The CBM experiment is a next generation fixed-target universal detector capable of measuring hadronic, 
leptonic and photonic probes in a large geometrical acceptance with good vertexing. Its modular 
configuration (Fig. [50]) will allow quarkonia measurements both in the electron and muon channels. 
Hadron-track reconstruction and momentum measurements will be performed by the Silicon Tracking 
System (STS) which consists of several stations of radiation hard silicon strips located inside the aperture 
of a 1 Tm bending dipole magnet. Displaced vertices will be measured with high precision by means of a 
Micro Vertex Detector (MVD) based on monolithic active pixel sensors operated in vacuum close to the 
target. Electrons will be identified by a combination of RICH and TRD detectors placed downstream. 
For muon measurements, the RICH detector will be replaced by an active absorber system (MUCH) 
made of sandwiches of detection planes and iron layers. The setup is completed by a RPC-based ToF 
system for charged-hadron measurements and by an EMC for neutral particle measurements. 

Figure [51] shows the J/ip signal extracted from 4 ■ 10^° central Au-Au collisions at 25 AGeV. The 
corresponding expected data taking time is 11 hours in the muon channel and 55 hours in the electron 
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Note, however, that the uncertainty from IQCD calculations on the phase boundary is still large. 
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Figure 50: Longitudinal view of the CBM detector. Left: configuration for electron and liadron mea- 
surements. Right: configuration for muon measurements. 



channel. The simulation assumes a vr^ suppression of 10'* from the RICH, TRD and ToF systems in 
the electron channel and a very good suppression of muons from vr^, K^ decay in the muon channel 
thanks to kink detection with the STS. The obtained reconstruction performance is comparable in both 
channels. 
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Figure 51: Invariant- mass spectrum showing J /ip reconstruction in central Au-Au collisions at 25 AGeV 
with the CBM detector in the electron (left) and muon (right) channels |288j . 



5 Summary 

The description of quarkonium production in ultrarelativistic heavy-ion collisions remains a challenging 
task, but new insights keep emerging on the three major frontiers, i.e., theoretical, phenomenological 
and experimental. In this review we have mostly focused on recent developments and have attempted to 
combine the current knowledge toward developing a consistent framework which can be checked against 
lattice QCD (IQCD) "data" and realistically applied to experimental data. In each of the following 



paragraphs, we briefly summarize the main points for the three frontiers and subsequently identify 
directions of future work. 

First principle IQCD computations of quarkonium correlation functions at finite temperature, now 
also available including dynamical light quarks, have consolidated the finding that ground-state char- 
monia (bottomonia) are little affected up to temperatures of ~2(3) T^. This translates into the presence 
of bound (or resonance) states in the pertinent spectral functions when extracted from the correlators 
using probabilistic methods. On the other hand, IQCD computations of the free energy of a heavy-quark 
(HQ) pair in a Quark-Gluon Plasma (QGP) clearly exhibit the effect of color screening, increasing with 
temperature. When these free energies are injected as a potential into a Schrodinger equation, the J/ip 
melts at significantly smaller temperatures than suggested by the correlators, slightly above T^. Higher 
melting temperatures are obtained from potential models when employing the IQCD-based internal 
energies, which apparently agrees better with the spectral function results. However, in both cases, i.e. 
using free or internal energies, potential models are able to reproduce the approximate constancy of the 
correlation functions. This has recast some doubt on whether ground-state charmonia really survive well 
above Tc. One of the main reasons for this redundancy is that a reduced binding can be compensated 
by a reduced open-charm threshold, corresponding to a reduced thermal charm-quark mass - an effect 
which is well established for bound states in electromagnetic plasmas. Furthermore, nonperturbative 
rescattering effects typically lead to a large enhancement in the continuum part of the quarkonium 
spectral functions, especially close to the quark- ant iquark threshold. Thus a reliable determination of 
the in-medium charm-quark mass, as well as of finite-width effects (for both charmonia and open-charm 
states), is needed to improve the evaluation of quarkonium spectral functions. In addition, it will be 
important to develop better criteria for the applicability of potential models at finite temperature, e.g., 
by setting up effective field theories if suitable scale separations can be applied. 

In-medium binding energies turn out to play an important role in quarkonium dissociation reactions. 
For strong binding, leading-order gluo-dissociation is the most important process. For weakly bound 
quarkonia, next-to-leading order processes ("quasifree" destruction) become dominant, especially close 
to the "melting" temperatures; both thermal quarks and gluons contribute to these reactions. The 
pertinent inelastic widths of charmonia can reach up to several hundred MeV for temperatures as 
low as 1.5 Tc (depending on the value of ag), translating into lifetimes (well) below Ifm/c. For phe- 
nomenology in heavy-ion collisions this implies that quarkonium dissociation will be operative well 
before the pertinent bound states dissolve due to color screening. It also reiterates that a consistent 
treatment of screening and finite-width effects is mandatory. Furthermore, the spectral properties of 
open heavy-flavor states are expected to change appreciably in the medium, with direct consequences 
for quarkonia: while effective HQ (or heavy-meson) masses determine the onset of the open heavy-flavor 
continuum, flnite widths can open additional quarkonium decay channels even below threshold. Taking 
guidance from well-developed concepts in electromagnetic plasma physics, flnite temperature T-matrix 
approaches seem to be a suitable tool to incorporate all of these aspects. 

Equilibrium properties of quarkonia in medium (masses and widths) are key ingredients to transport 
theoretical evaluations of quarkonium observables in heavy-ion collisions, especially if carried out in a 
thermally evolving bulk medium (which appears to be a good approximation at SPS and RHIC). A 
quantitative understanding of cold-nuclear-matter effects on initial quarkonium spectra is mandatory 
to extract "anomalous" suppression in the subsequently formed hot and dense medium. Cold-nuclear- 
matter effects include modiflcations of the nuclear parton distribution functions ("shadowing"), initial- 
state energy loss, Cronin effect and nuclear absorption. Experimental evidence is mounting that the 
relative importance of these effects changes signiflcantly with collision energy and rapidity, and thus 
empirical constraints from p-A (or d-A) data are indispensable. In addition, secondary charmonium pro- 
duction ("regeneration") via c-c coalescence (in the QGP or at hadronization) is expected to become 
increasingly relevant with increasing collision energy. This implies a coupling of charmonium yields 
and spectra to the abundance and spectra of open charm, which themselves are subject to medium 



modifications. This further adds to the complexity of the calculations reinforcing the need of compre- 
hensive analyses of open- and hidden-flavor observables in collision centrality, energy, system size and 
3-momentum. The current state of phenomenology suggests that (a) the observed J /if) suppression at 
the SPS mostly occurs in a hot/dense medium at energy densities above the critical one found in IQCD, 
with initial temperatures of Tq ^220-250 MeV; (b) at RHIC energies, a significant regeneration compo- 
nent develops, accounting for roughly 50% of the measured J /ip yield in central Au-Au; this assertion 
is corroborated by the observed narrowing in the rapidity distributions and an approximately constant 
average transverse momentum with centrality, with regeneration contributions prevalent at low px- If 
the regeneration yield becomes dominant at the LHC, a qualitatively new (increasing) centrality de- 
pendence of the nuclear modification factor may emerge (depending on the open-charm cross section). 
A more direct probe of color screening is provided by the T: if it retains its vacuum binding energy 
up to ~2Tc, it will be rather inert under RHIC conditions; however, an in-medium reduction in its 
binding increases the dissociation rates noticeably which might suppress T production as much as J /ip 
production. Such a signal benefits from the absence of bottomonium regeneration at RHIC. A similar 
feature could be even more pronounced at the LHC, especially if regeneration of the J /ip is large, while 
it is still expected to be small for T. Transverse-momentum spectra will further illuminate the different 
mechanisms, and/or reveal new ones. E.g., first RHIC data for J/ip^s at high pt > 5GeV indicate a 
nuclear modification factor close to one, very different from light-hadron spectra. Elliptic flow could 
be another good discriminator between direct and secondary production, but only if the regeneration 
component in semi-central collision persists to sufficiently large px, above at least 2GeV. Finally, ex- 
cited \E' and T states are of great interest to test suppression patterns as regeneration contributions are 
usually suppressed due to their larger mass. 

The quantitative realization of the above observables requires an advanced experimentation which 
will become available at RHIC-II and the LHC. While the former appears to operate at an energy 
suited to study the properties of a strongly coupled, liquid-like QGP, the latter may, for the first 
time, reach into a weakly coupled gas-like regime. At lower energies, the CBM experiment at FAIR 
will explore strongly interacting matter at maximum net baryon density; little (nothing) is known 
theoretically (experimentally) about charmonium properties under these conditions. In view of these 
developments, and with an adequate theoretical support, quarkonia promise to remain a rich, and 
eventually quantitative, probe of high-energy density matter in heavy-ion collisions. 
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